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)
      ccTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCGTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcg                                                                                                                                           >  1:775969/1‑149 (MQ=255)
        ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCAtcgtc                                                                                                                                         >  2:406042/1‑149 (MQ=255)
            aTCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGata                                                                                                                                     <  1:115373/149‑1 (MQ=255)
                tacgggCGTGGAAGTGGGTAGGGAACTGGTGATGGAAATTTTGTCGGTTGTATCGGGCAAGGAAGGGGTGGTGGGTCGAGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTa                                                                                                                                 <  1:208726/146‑1 (MQ=255)
                        tGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCggct                                                                                                                         >  2:842670/1‑147 (MQ=255)
                          gAACTGGGTAAGGAACTGGTGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGGCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATaa                                                                                                                       <  2:883221/149‑1 (MQ=255)
                           aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGCGCAAGGAAGGGGCGAGCGGGCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAACCTTGATCATCGCCGAGATATAGACCCGATaca                                                                                                                      >  2:207590/1‑147 (MQ=255)
                                gggaaGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGtt                                                                                                                 <  2:117821/146‑1 (MQ=255)
                                ggTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGACGGAAAGCTTGATCATCGTCGATATGTAGGCCGGACAAGGCGtt                                                                                                                 >  2:403101/1‑149 (MQ=255)
                                  tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAAGCGTTCa                                                                                                               >  2:821779/1‑149 (MQ=255)
                                  tAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAAAGGCGAAGGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGCTCTTCATCGATATGTAAGCCGAATAAGGCGTTCa                                                                                                               >  2:371805/1‑149 (MQ=255)
                                     ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGTAAGGGTTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                            <  1:848278/149‑1 (MQ=255)
                                      aaaCTGGGGACGGAACTTGTGCCTGTTGTATCGGGCAAGTAAGGGGCGAGTGGTCGCTACTCTTCCACGCAAGGGCTTGTCGTGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                           <  1:533302/148‑1 (MQ=255)
                                            gCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAGGGCGTTCACGCCGCATcc                                                                                                     <  1:406042/149‑1 (MQ=255)
                                             ggacggacgTTTGGCCGGTTGTATGGGGGAAGGACGGGGAGCGTGGCCGAGGAGTGTGCAGGCTAGTGGTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  2:638751/148‑1 (MQ=19)
                                             cGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCTAGGGTTTGTTGCGTATCTGCATGCGCGAAGGAAACCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCg                                                                                                    <  1:802829/149‑1 (MQ=39)
                                                 ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGTATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                >  1:296533/1‑149 (MQ=255)
                                                  gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                               >  2:152935/1‑149 (MQ=255)
                                                   aaCTTTTGCCGGTTGTATCGGGCAAGGAAGGGGTGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAcc                                                                                              <  1:726631/149‑1 (MQ=255)
                                                      ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGCTATGTAGGCCAGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                           >  1:49957/1‑149 (MQ=255)
                                                        ttGCCGGTTGTATCGGGCAAGGAAGGGGGGAGTGGTCGCGACGATTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGTGAAGGAAAGCTTGATCCTCGTCTTTATGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCAACCGATGc                                                                                         <  2:48905/149‑1 (MQ=38)
                                                          gCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCCCGCCTCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATAATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCt                                                                                       <  1:165159/149‑1 (MQ=39)
                                                             ggTTGTATCGGGCAAGGAAGGTGCGAGTGGTCCCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                    >  1:354457/1‑149 (MQ=255)
                                                             ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTTGCGACGTTTCCATTCTAGGGCTTTTTGTGCTTCTGTATGCGTGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                    <  1:295024/149‑1 (MQ=37)
                                                              gTTGTATGGGGTAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGAGAAGGAAAGCTTGATCATCTTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                   <  2:145962/149‑1 (MQ=39)
                                                              gTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATg                                                                                   <  2:372439/149‑1 (MQ=255)
                                                              gTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGACGATATGTAAGCCGGATAAAGCGTTCACGCCGCATCCGGCAACCGATACCTGATg                                                                                   >  1:468461/1‑149 (MQ=255)
                                                                   cgcGGGCAAGGAAGGGGCTAGTGGTCGACACCGTGCCACTCAAGGCCTTGTCGCGCATCTGCATGCGCGAAGGAACGCTTGCTCCTTGTCGATATGTAGGCCGGCGAACGCGTTTATCCCGCATCCGGCAATTGATGCCTGATGCGTcg                                                                              <  2:184267/147‑1 (MQ=255)
                                                                          aaggaaggGGTGAGTGGTCGCGTCGCTTCCACGCAAGGGCTTGTCGCGTATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                       <  1:870016/149‑1 (MQ=255)
                                                                          aaggaaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCACATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                       >  1:531336/1‑149 (MQ=255)
                                                                           aggaaggGGCGAGTGGTCGCGACGCTTCCGCGCAAGGGATTGTCGCGCATATGCGTGAGCGTAGGAAAGCTTGATCATCGTCGATAAGTAGGCCGGATAAGGCATTCAAGCCGCATCAGGCAACCGATGCCTGATGCGATGCGGAcgcg                                                                      >  1:760461/1‑149 (MQ=255)
                                                                              agggggCGGGGGGTGGCGACGCTTCCACGGAAGGGCTAGTCGGGCATTTGTAGGTGCGAAGGGAAGCTTGATCGTGGCCGATCTGTAGGCCGGTTAAGGTGTTCACGCCGCTTCCGGCAACGGATGCCTGATGCGACGCGGTCGCGTCt                                                                   <  2:176600/147‑1 (MQ=255)
                                                                              aaggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCACATCTGCATGAGCGAAGGAAAACATGAGCATCGTCGATATGTAGGCCGGATAAGCGGTTCACGACGCATCAGGCAACCGATGCCTGCTGAGACGCGGTCGCGTCt                                                                   >  2:433674/1‑149 (MQ=255)
                                                                                   gcggggggTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  2:537388/143‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGTAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  2:876794/149‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  1:633215/149‑1 (MQ=255)
                                                                                   gCGAGTGGTCGCGACCCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGTCTGATGCGACGCGGTCGCGTCTTATCa                                                              <  1:823255/149‑1 (MQ=255)
                                                                                      agtgtacgcgACGCTTCCACGCAAGTGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGGCGGAGAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGc                                                           <  1:767644/143‑1 (MQ=38)
                                                                                        tGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCt                                                         >  1:696974/1‑149 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTAc                                                       >  1:446598/1‑149 (MQ=255)
                                                                                          gTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCAAGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGCCGCGGTCGCGTCTTATCAGGCCCAc                                                       >  2:565209/1‑149 (MQ=255)
                                                                                           tCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGCCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGATGCCCGATGCGACGCGGTCGCATTTTATCAGCGCTACa                                                      >  1:790857/1‑149 (MQ=38)
                                                                                             gcgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCCACGCGGTCGCGTCTTATCAGGCCTACAgg                                                    >  2:233198/1‑149 (MQ=255)
                                                                                              cgAGGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                   <  2:280967/149‑1 (MQ=255)
                                                                                              cgACGCTTCCCCGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCTATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGt                                                   <  1:57651/149‑1 (MQ=255)
                                                                                              cgACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCCTATCAGGCCTACAgtc                                                   >  2:850241/1‑147 (MQ=255)
                                                                                                 cGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGACAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGa                                                <  2:863498/149‑1 (MQ=255)
                                                                                                  gCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGAt                                               <  2:296533/149‑1 (MQ=255)
                                                                                                   cgtCCCCGCACGGGCTTGCCGCGCATCTGCATCCGCGAAGGAAAGCTTGATCATCGTTGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCCTATCAGGCCTACAGGTCGATg                                              <  2:616407/147‑1 (MQ=25)
                                                                                                          gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATCTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATACCTGATGCGACGCGGTCGCGTCTTACCAGGCCTACAGGTGGAGGCAGaaaa                                       >  1:61027/1‑146 (MQ=25)
                                                                                                             aGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGAAATGTa                                    >  1:529123/1‑149 (MQ=255)
                                                                                                                   tGTCGCGCATCTGTGTTGTCGAAGGAAAGCTTTATCATCTTCTATTTGTATGCCGGATAAGGCTTTCACGCCGCATCCGGCAACCGATGCCTCATGCGACGCGGTCGCGTCTTATCAGGTCGACAGGTCGATGCCGATATGTACATCGt                              <  2:734739/149‑1 (MQ=255)
                                                                                                                      cgccaaTCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGAAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTAtt                           <  1:293062/144‑1 (MQ=255)
                                                                                                                       gcgcAGCTGCATGCGGGAGGGAAAGCTTGAACATCGTCGATTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAACGATGCCTGATGCGACGCGGGCTCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTc                          <  2:326993/149‑1 (MQ=25)
                                                                                                                          cATCTGCATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGACCTACAGGTCGATGCCGATATGTACATCGTATTCGGc                       >  1:687481/1‑149 (MQ=255)
                                                                                                                            tgtGCATACGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCaa                     <  2:790857/147‑1 (MQ=255)
                                                                                                                                cATGCGCGAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTaa                 <  1:152935/149‑1 (MQ=255)
                                                                                                                                  tGCGCGAAGGAAAGCTGGAGCATCGTCGATAGGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTTATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATTCCGATATGTACCTCGTATGCGGCAAGTAATa               <  2:329849/149‑1 (MQ=25)
                                                                                                                                       gAAGGAAAGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGACTTATCGGGCCTACAGGTCGATGCCGATATGTACATCGCATTCGCCAATTAATACAtca          >  2:775747/1‑147 (MQ=37)
                                                                                                                                              aGCTTGATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTCATACATAACACgatt   >  2:128451/1‑149 (MQ=255)
                                                                                                                                               gCTTGATCCTCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGAAAGGTACATCGTAATCGGCAACTAATGCATAGCCCgattg  >  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: