Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_bme_pgi 4,235,408 Δ1 bp intergenic (‑2175/‑499) lysC ← / → yjbE lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,235,4070A.98.4% 282.0 / ‑1.2 64intergenic (‑2174/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (27/35);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCA  >  NC_000913_3_bme_pgi/4235259‑4235545
                                                                                                                                                     |                                                                                                                                          
tCAAGGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGa                                                                                                                                            <  2:11979/149‑1 (MQ=255)
  aaggCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGGGCCGGAACTTTTGCCGGTTGTAGCGGGCAAGGAAGGGGCGCGTGGTCGCGACGCTTCCACGCAAGGGGTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGatc                                                                                                                                          <  2:867500/146‑3 (MQ=255)
  aaTGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGatc                                                                                                                                          >  1:549621/1‑147 (MQ=255)
        gtCGATCACTGGGGCGTGGAACTGGGAAAGGAACTGGGGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCAtcgtc                                                                                                                                   <  2:158930/148‑1 (MQ=255)
            aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACCCTCCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGata                                                                                                                               <  1:340023/149‑1 (MQ=255)
                aTGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTa                                                                                                                           <  2:978532/149‑1 (MQ=255)
                 tGGGGGGTGGAACTGGGTAAGGAACTGGGGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTCCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAg                                                                                                                          <  2:877866/149‑1 (MQ=255)
                  ggggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAgg                                                                                                                         >  1:713846/1‑149 (MQ=255)
                         ggAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATa                                                                                                                  >  2:477382/1‑149 (MQ=255)
                           aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                <  1:758741/149‑1 (MQ=255)
                            aCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAgg                                                                                                               <  1:290585/149‑1 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCTAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                           >  1:664470/1‑149 (MQ=255)
                                  tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGttaa                                                                                                         >  2:366128/1‑147 (MQ=255)
                                     ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                      <  2:549066/149‑1 (MQ=255)
                                     ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                      <  1:284477/149‑1 (MQ=255)
                                            gCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                               >  1:554594/1‑149 (MQ=255)
                                               aCGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                            >  2:734696/1‑149 (MQ=255)
                                               aCGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                            >  2:735419/1‑149 (MQ=255)
                                                 ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGAGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGAGAAGGCGTTCCCGCCGCATCCGGCaa                                                                                          >  1:983451/1‑149 (MQ=255)
                                                  gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                         >  1:438957/1‑149 (MQ=255)
                                                   aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                        <  1:965212/149‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                     <  1:507671/149‑1 (MQ=255)
                                                          gCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                 <  2:371586/149‑1 (MQ=255)
                                                            cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                               >  2:456530/1‑149 (MQ=255)
                                                            cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCAACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                               >  1:156603/1‑149 (MQ=255)
                                                             ggTTGTATCGGGCAAGGAAGGAACGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCAAGCGCGA‑AGGAAAGCTTG‑ATAATCGTCGATATGTCGGCCGGATACGGCGTTCCCGCCGCATCCGTGCACCGATGCCtcag                                                                              >  2:436698/1‑146 (MQ=255)
                                                               ttGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGc                                                                            >  1:262157/1‑149 (MQ=255)
                                                                  tATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCACGCATCTGCATGCGCGA‑GGGAAAGGTTG‑AGCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCACCCGATGCCTGCTGAGAc                                                                         >  2:192341/1‑149 (MQ=38)
                                                                      gggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCgg                                                                     >  2:671755/1‑149 (MQ=255)
                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGatc                                                                                                                                          >  1:864844/1‑78 (MQ=255)
                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGatc                                                                                                                                          <  2:864844/80‑3 (MQ=255)
                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAGCGATGCCTGGTGCGACGCggg                                                                    >  1:121479/1‑148 (MQ=255)
                                                                       ggCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGt                                                                    >  1:706915/1‑149 (MQ=255)
                                                                          aaggaaggGTCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                 <  1:708714/149‑1 (MQ=255)
                                                                           aggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                <  1:714385/149‑1 (MQ=255)
                                                                           aggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                <  1:714430/149‑1 (MQ=255)
                                                                            ggaaggGGCGAGTGGTCGCGACGCTTCCACGCATGCGCTTTTATTACTTTTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                               <  2:464073/149‑1 (MQ=25)
                                                                            ggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGt                                                               <  2:868186/149‑1 (MQ=255)
                                                                              aaggGGCGAGTGGTCCCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAGGCCTGATGCGACGCGGTCGCGTCt                                                             <  2:619049/149‑1 (MQ=255)
                                                                               aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                            >  2:340226/1‑149 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                >  1:33823/1‑141 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgcg                                                                <  2:33823/141‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                        <  2:706915/149‑1 (MQ=255)
                                                                                    cggtgggTGGCGACGCTTCTACGCAAGGGCTTGTCGCGCATTTGCATGCGCGA‑AGGAAAGCTTG‑CTCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGGATGCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                       <  2:408041/144‑1 (MQ=25)
                                                                                    cgggTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                       <  1:68482/146‑1 (MQ=255)
                                                                                    cGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAg                                                       <  1:340226/149‑1 (MQ=255)
                                                                                      aGTGGTCGCCGCGCTTCCACGCCAGGGAAGGTCGCGCAGTTGCATGCGCGA‑AGGAAAGCGTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                     <  2:664470/149‑1 (MQ=25)
                                                                                         tgTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCCTCAGGGC‑AGG‑TTGCCTATCTCAGCGTGTATATGTAGGCCGGATAATGCGTTCACTCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatc                                                  <  2:960493/148‑4 (MQ=12)
                                                                                         ggTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTa                                                  <  1:477382/149‑1 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                 >  2:496667/1‑149 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATAAGACGCGACCGCGTCGCATCAGGcatcga                                                >  1:960493/1‑144 (MQ=38)
                                                                                            cgcgACGCCTCCACGCAAGGGCTTGTCGCGCATCTGAATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAg                                               <  2:334649/149‑1 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCAGCATGCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                              >  1:443991/1‑149 (MQ=255)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGGCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAgg                                              <  2:438957/149‑1 (MQ=255)
                                                                                              ggACGCATCCACGCAAGGGCTGGTCGAGGACCTGCATCCGCAACCGGAACGCT‑‑ATTTCTCGTCGATATGTAGGCCGGACAAGGCGTTCGCGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                             <  1:282686/148‑1 (MQ=17)
                                                                                                    tgccACCCAATGTCTTGTCGTGCATGTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCTTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGc                                       <  2:198637/147‑1 (MQ=25)
                                                                                                      ccACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                     >  1:310590/1‑85 (MQ=255)
                                                                                                      ccACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                     <  2:310590/85‑1 (MQ=255)
                                                                                                             aGGGCTTGTCGCGCATCTGAATGCGCGA‑AGGAATGCTTG‑ATCTTCGTCGATATGTAGGCCGGATAAGGCGTTCACGTCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATACCGATATGTa                              <  1:422512/149‑1 (MQ=25)
                                                                                                                gCTTGCCGCGCATTTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACAt                           <  1:291334/149‑1 (MQ=255)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTTGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                  <  2:401376/140‑1 (MQ=255)
                                                                                                                  ttGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTTGCGTCTTATCAGGCCTACAGGTCGATGCCGata                                  >  1:401376/1‑140 (MQ=255)
                                                                                                                    gTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑atcatc                                                                                                                                      <  2:609625/38‑1 (MQ=255)
                                                                                                                    gTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTG‑atcatc                                                                                                                                      >  1:609625/1‑38 (MQ=255)
                                                                                                                      cgcgCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                              <  1:34245/108‑1 (MQ=255)
                                                                                                                      cgcgCATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTc                                                              >  2:34245/1‑108 (MQ=255)
                                                                                                                       gcgcATCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAttg                    >  2:607864/1‑148 (MQ=255)
                                                                                                                           aTCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATACGGCATGCGCTCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                <  1:735419/149‑1 (MQ=37)
                                                                                                                           aTCTGCATGCGCGA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGCCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCa                <  1:734696/149‑1 (MQ=255)
                                                                                                                                       gA‑AGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAg    >  1:817836/1‑149 (MQ=255)
                                                                                                                                          aGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa  <  1:456530/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                          
TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGA‑AGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCA  >  NC_000913_3_bme_pgi/4235259‑4235545

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑