Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_bme_pgi 4,235,407 Δ1 bp 100% intergenic (‑2174/‑500) 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.95.1% 240.1 / ‑0.2 61intergenic (‑2174/‑500)lysC/yjbElysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  major base . (27/31);  minor base T (0/1);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_bme_pgi/4235259‑4235551
                                                                                                                                                    |                                                                                                                                                
tCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTACAGGTTGTGTCGGGCAAGGAAGGGGCGAGTGGGCGCGACGCTTCCACGAACGGGATTCTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                                  >  1:208843/1‑149 (MQ=255)
      ccTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCGTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcg                                                                                                                                           >  1:775969/1‑149 (MQ=255)
        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                         >  2:406042/1‑149 (MQ=255)
            aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGata                                                                                                                                     <  1:115373/149‑1 (MQ=255)
                tacgggCGTGGAAGTGGGTAGGGAACTGGTGATGGAAATTTTGTCGGTTGTATCGGGCAAGGAAGGGGTGGTGGGTCGAGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTa                                                                                                                                 <  1:208726/146‑1 (MQ=255)
                        tGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCggct                                                                                                                         >  2:842670/1‑147 (MQ=255)
                          gAACTGGGTAAGGAACTGGTGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATaa                                                                                                                       <  2:883221/149‑1 (MQ=255)
                           aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGCGCAAGGAAGGGGCGAGCGGGCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAACCTTG‑ATCATCGCCGAGATATAGACCCGATaca                                                                                                                      >  2:207590/1‑147 (MQ=255)
                                gggaaGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                 <  2:117821/146‑1 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGACGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGACAAGGCGtt                                                                                                                 >  2:403101/1‑149 (MQ=255)
                                  tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAAGCGTTCa                                                                                                               >  2:821779/1‑149 (MQ=255)
                                  tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAAAGGCGAAGGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑CTCTTCATCGATATGTAAGCCGAATAAGGCGTTCa                                                                                                               >  2:371805/1‑149 (MQ=255)
                                     ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGTAAGGGTTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                            <  1:848278/149‑1 (MQ=255)
                                      aaaCTGGGGACGGAACTTGTGCCTGTTGTATCGGGCAAGTAAGGGGCGAGTGGTCGCTACTCTTCCACGCAAGGGCTTGTCGTGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                           <  1:533302/148‑1 (MQ=255)
                                            gCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAGGGCGTTCACGCCGCATcc                                                                                                     <  1:406042/149‑1 (MQ=255)
                                             ggacggacgTTTGGCCGGTTGTATGGGGGAAGGACGGGGAGCGTGGCCGAGGAGTGTGCAGGCTAGTGGTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  2:638751/148‑1 (MQ=19)
                                             cGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCTAGGGTTTGTTGCGTATCTGCATGCGCGAAGGAAACCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  1:802829/149‑1 (MQ=39)
                                                 ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                >  1:296533/1‑149 (MQ=255)
                                                  gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               >  2:152935/1‑149 (MQ=255)
                                                   aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGTGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                              <  1:726631/149‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGCTATGTAGGCCAGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                           >  1:49957/1‑149 (MQ=255)
                                                        ttGCCGGTTGTATCGGGCAAGGAAGGGGGGAGTGGTCGCGACGATTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGTGAAGGAAAGCTTG‑ATCCTCGTCTTTATGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                         <  2:48905/149‑1 (MQ=38)
                                                          gCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCCCGCCTCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATAATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                       <  1:165159/149‑1 (MQ=39)
                                                             ggTTGTATCGGGCAAGGAAGGTGCGAGTGGTCCCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                    >  1:354457/1‑149 (MQ=255)
                                                             ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTTGCGACGTTTCCATTCTAGGGCTTTTTGTGCTTCTGTATGCGTGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                    <  1:295024/149‑1 (MQ=37)
                                                              gTTGTATGGGGTAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGAGAAGGAAAGCTTG‑ATCATCTTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                   <  2:145962/149‑1 (MQ=39)
                                                              gTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                   <  2:372439/149‑1 (MQ=255)
                                                              gTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGACGATATGTAAGCCGGATAAAGCGTTCACGCCGCATCCGGCAACCGATACCTGATg                                                                                   >  1:468461/1‑149 (MQ=255)
                                                                   cgcGGGCAAGGAAGGGGCTAGTGGTCGACACCGTGCCACTCAAGGCCTTGTCGCGCATCTGCATGCGCGAAGGAACGCTTG‑CTCCTTGTCGATATGTAGGCCGGCGAACGCGTTTATCCCGCATCCGGCAATTGATGCCTGATGCGTcg                                                                              <  2:184267/147‑1 (MQ=255)
                                                                          aaggaaggGGTGAGTGGTCGCGTCGCTTCCACGCAAGGGCTTGTCGCGTATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                       <  1:870016/149‑1 (MQ=255)
                                                                          aaggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCACATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                       >  1:531336/1‑149 (MQ=255)
                                                                           aggaaggGGCGAGTGGTCGCGACGCTTCCGCGCAAGGGATTGTCGCGCATATGCGTGAGCGTAGGAAAGCTTG‑ATCATCGTCGATAAGTAGGCCGGATAAGGCATTCAAGCCGCATCAGGCAACCGATGCCTGATGCGATGCGGAcgcg                                                                      >  1:760461/1‑149 (MQ=255)
                                                                              agggggCGGGGGGTGGCGACGCTTCCACGGAAGGGCTAGTCGGGCATTTGTAGGTGCGAAGGGAAGCTTG‑ATCGTGGCCGATCTGTAGGCCGGTTAAGGTGTTCACGCCGCTTCCGGCAACGGATGCCTGATGCGACGCGGTCGCGTCt                                                                   <  2:176600/147‑1 (MQ=255)
                                                                              aaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCACATCTGCATGAGCGAAGGAAAACATG‑AGCATCGTCGATATGTAGGCCGGATAAGCGGTTCACGACGCATCAGGCAACCGATGCCTGCTGAGACGCGGTCGCGTCt                                                                   >  2:433674/1‑149 (MQ=255)
                                                                                   gcggggggTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  2:537388/143‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGTAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  2:876794/149‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  1:633215/149‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACCCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGTCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  1:823255/149‑1 (MQ=255)
                                                                                      agtgtacgcgACGCTTCCACGCAAGTGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGGCGGAGAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                           <  1:767644/143‑1 (MQ=38)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                         >  1:696974/1‑149 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                       >  1:446598/1‑149 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCAAGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCTTATCAGGCCCAc                                                       >  2:565209/1‑149 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGCCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGATGCCCGATGCGACGCGGTCGCATTTTATCAGCGCTACa                                                      >  1:790857/1‑149 (MQ=38)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCCACGCGGTCGCGTCTTATCAGGCCTACAgg                                                    >  2:233198/1‑149 (MQ=255)
                                                                                              cgAGGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                   <  2:280967/149‑1 (MQ=255)
                                                                                              cgACGCTTCCCCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCTATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                   <  1:57651/149‑1 (MQ=255)
                                                                                              cgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCCTATCAGGCCTACAgtc                                                   >  2:850241/1‑147 (MQ=255)
                                                                                                 cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGACAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                                <  2:863498/149‑1 (MQ=255)
                                                                                                  gCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                               <  2:296533/149‑1 (MQ=255)
                                                                                                   cgtCCCCGCACGGGCTTGCCGCGCATCTGCATCCGCGAAGGAAAGCTTG‑ATCATCGTTGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCCTATCAGGCCTACAGGTCGATg                                              <  2:616407/147‑1 (MQ=25)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATACCTGATGCGACGCGGTCGCGTCTTACCAGGCCTACAGGTGGAGGCAGaaaa                                       >  1:61027/1‑146 (MQ=25)
                                                                                                             aGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGAAATGTa                                    >  1:529123/1‑149 (MQ=255)
                                                                                                                   tGTCGCGCATCTGTGTTGTCGAAGGAAAGCTT‑TATCATCTTCTATTTGTATGCCGGATAAGGCTTTCACGCCGCATCCGGCAACCGATGCCTCATGCGACGCGGTCGCGTCTTATCAGGTCGACAGGTCGATGCCGATATGTACATCGt                              <  2:734739/149‑1 (MQ=255)
                                                                                                                      cgccaaTCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGAAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                           <  1:293062/144‑1 (MQ=255)
                                                                                                                       gcgcAGCTGCATGCGGGAGGGAAAGCTTGAA‑CATCGTCGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGATGCCTGATGCGACGCGGGCTCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                          <  2:326993/149‑1 (MQ=25)
                                                                                                                          cATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGACCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                       >  1:687481/1‑149 (MQ=255)
                                                                                                                            tgtGCATACGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                     <  2:790857/147‑1 (MQ=255)
                                                                                                                                cATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                 <  1:152935/149‑1 (MQ=255)
                                                                                                                                  tGCGCGAAGGAAAGC‑TGGAGCATCGTCGATAGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTTATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATTCCGATATGTACCTCGTATGCGGCAAGTAATa               <  2:329849/149‑1 (MQ=25)
                                                                                                                                       gAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGACTTATCGGGCCTACAGGTCGATGCCGATATGTACATCGCATTCGCCAATTAATACAtca          >  2:775747/1‑147 (MQ=37)
                                                                                                                                              aGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTCATACATAACACgatt   >  2:128451/1‑149 (MQ=255)
                                                                                                                                               gCTTG‑ATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGAAAGGTACATCGTAATCGGCAACTAATGCATAGCCCgattg  >  2:391440/1‑149 (MQ=21)
                                                                                                                                                    |                                                                                                                                                
TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGATTG  >  NC_000913_3_bme_pgi/4235259‑4235551

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