Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,296,0610CT33.3% 109.2 / 73.5 96intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (37/27);  new base T (19/13);  total (56/40)
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.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913_3_bme_pgi/4295913‑4296209
                                                                                                                                                    |                                                                                                                                                    
gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                                                                                                     >  1:1044624/1‑150 (MQ=32)
  cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                   <  2:807867/150‑1 (MQ=17)
   ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                  <  1:909416/150‑1 (MQ=17)
    gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                 >  1:1094071/1‑150 (MQ=17)
     aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTgc                                                                                                                                                >  2:1131330/1‑149 (MQ=17)
     aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATATAGCGTCATACCGTAGGTCGGATAATGCGTTTACGCCGCATCCGACATCAATCCCTGa                                                                                                                                                >  2:1338993/1‑150 (MQ=11)
     aTAAGGCGTTTATGCCCCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                <  1:923753/150‑1 (MQ=17)
      tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGGTCGGGTAAGGCGTTTACGCCGCATCCGCCATCAACGCCCGac                                                                                                                                               >  1:849650/1‑149 (MQ=255)
         ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTCATTTTGTAGGTCGGATAAGGCCTTTACGCCCCCTCCGACATCAATTCATGATGCg                                                                                                                                            >  2:565750/1‑150 (MQ=11)
         ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCg                                                                                                                                            >  2:745056/1‑150 (MQ=17)
          gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGCCTGATCTAGCGTGATTTTGTAGCTCGGCTAAGCCGTTTCCGCCGCATTCGACATCAATGCCTGAcgcgg                                                                                                                                           >  2:277495/1‑149 (MQ=11)
          gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCCGATGCgc                                                                                                                                           >  2:140876/1‑149 (MQ=12)
            gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                         >  1:94057/1‑150 (MQ=17)
            gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                         >  1:947465/1‑150 (MQ=17)
             tttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGc                                                                                                                                        >  2:765636/1‑150 (MQ=32)
              ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                                                                                                       >  2:508785/1‑150 (MQ=17)
               tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCtt                                                                                                                                      <  2:1044624/150‑1 (MQ=32)
               cacgccgccTCCGAAACCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCtt                                                                                                                                      <  1:403167/147‑1 (MQ=11)
                aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                     >  2:1181305/1‑150 (MQ=32)
                  gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGTCAGACAGCGCAGTAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcg                                                                                                                                   >  1:829185/1‑149 (MQ=255)
                   ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTATGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcg                                                                                                                                  >  1:923377/1‑150 (MQ=25)
                    cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                 >  1:1253412/1‑150 (MQ=25)
                      tatcgagcATCAACGCCTGAGGGGACGAATGCCGCGTCTTATCAGGCCTACGTCAGACAGCGAAATAGGCTGATTTAGCGTGAATTTGTAGGTCGGATAAGGCGTTTACGCCGACACCGTTATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                               <  2:304191/143‑1 (MQ=11)
                       aTCCGACATCAACGCCTGATGCGACGCTTAACGCTTCTTAACAGACATACGCTACACAGAGAAATACACAAATGCAGCGTGACTTTGTCTACCGGATAAGGCGCCAACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  2:1317343/150‑1 (MQ=2)
                           gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              >  2:1068958/1‑146 (MQ=38)
                           gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  1:1068958/146‑1 (MQ=39)
                           gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  1:89981/1‑150 (MQ=37)
                              aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAGGCTTGCCGCGTCTTATCAg                                                                                                                       >  1:880550/1‑150 (MQ=37)
                              aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGACTTAGCGTGGTTTTGTAGGTCGGCTAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCCTCTTATAAg                                                                                                                       >  2:315940/1‑150 (MQ=11)
                               tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCCTCCGACATCAATGCCTGACGCGACGCTTGCCGCGTCCTCTCAgg                                                                                                                      >  2:739037/1‑150 (MQ=11)
                                cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAACGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                     <  1:810253/150‑1 (MQ=255)
                                 aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:245850/1‑150 (MQ=17)
                                  aagCCTGATGAGACGCTGAACGGGTGTAAAAAGGCGTACTCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGCGTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                   <  2:644580/148‑1 (MQ=11)
                                   cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       >  2:510742/1‑145 (MQ=255)
                                   cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  1:510742/145‑1 (MQ=255)
                                   cGCCTGAGGCGACTCTTAACGCGTCTTATCAGGCCTACGCCAGACCGAAAATCTGGCGGATTTAGCGTGATTTTTTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                  <  1:265184/150‑1 (MQ=11)
                                      cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                               >  2:547014/1‑147 (MQ=17)
                                       tGATGCGACGCTTAACGCGTGTGATGAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                              <  2:1306028/150‑1 (MQ=255)
                                       tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       >  1:1082718/1‑141 (MQ=38)
                                       tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  2:1082718/141‑1 (MQ=39)
                                         aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatcttac                                                                                                            >  1:69228/1‑144 (MQ=9)
                                                  ttAACGCGTCTTATCAGGCCTACGCCAGAGAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGTTTGCCGCGTCTTATCCGGCCGACCCCCGACAACGCaa                                                                                                   >  2:1262227/1‑150 (MQ=11)
                                                  ttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATTAGGCCTACGTCAGCCAGCGCaa                                                                                                   >  2:1017771/1‑150 (MQ=255)
                                                        ggtggtATCAGGCCAGCTCCAGACAGCGCAAGAGCCTGATTTAGCGTGAGTTTGTAGGCCGGATAAGGGGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                             <  2:93555/145‑1 (MQ=11)
                                                            ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                   >  2:1362966/1‑92 (MQ=21)
                                                            ttATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGcc                                                                                                                                                   <  1:1362966/92‑1 (MQ=21)
                                                             tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGCCAGCGCACTAGCCTGAttt                                                                                        >  2:1141354/1‑150 (MQ=255)
                                                                cAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCCCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                     <  1:301605/150‑1 (MQ=38)
                                                                 aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGGTGCGACGCTTGCCGCGTCTTATCAGGCCTAGGCCAGACAGCGCATTAGCCTGAATTAGCg                                                                                    >  1:132245/1‑150 (MQ=255)
                                                                   gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                               >  2:294652/1‑121 (MQ=18)
                                                                   gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                             <  1:652320/91‑1 (MQ=33)
                                                                   gCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                               <  1:294652/121‑1 (MQ=40)
                                                                    ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                    <  1:991146/83‑1 (MQ=25)
                                                                     cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                             >  2:115094/1‑89 (MQ=33)
                                                                     cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                             <  1:115094/89‑1 (MQ=35)
                                                                      tACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGACCGCGCAATAGCTTTGTTTAGCGTCAtt                                                                               >  2:621978/1‑150 (MQ=11)
                                                                         gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                >  2:1158889/1‑82 (MQ=18)
                                                                         gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                <  1:1158889/82‑1 (MQ=25)
                                                                           cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                  <  1:1343259/78‑1 (MQ=255)
                                                                                aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                           >  1:909130/1‑112 (MQ=18)
                                                                                aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                           <  2:909130/112‑1 (MQ=255)
                                                                                aGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCg                                                                     <  2:88680/150‑1 (MQ=38)
                                                                                 gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCgg                                                                    >  1:88680/1‑150 (MQ=39)
                                                                                    cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                     <  2:1231560/130‑1 (MQ=21)
                                                                                    cAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                     >  1:1231560/1‑130 (MQ=18)
                                                                                       tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATAATCCTGCTTTAGCGTGCTTTTGTATGTCGGATAAgg                                                              >  2:1297489/1‑150 (MQ=14)
                                                                                          ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACGGCGCAATAGCCTGCTTTCGCGTGATTTTGTAGGTCGGACAATGCGt                                                           >  2:579338/1‑150 (MQ=11)
                                                                                                    cGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                         <  1:1048245/126‑1 (MQ=40)
                                                                                                      tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                                <  2:1375845/53‑1 (MQ=255)
                                                                                                      tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                >  1:1375845/1‑53 (MQ=38)
                                                                                                         ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                           >  2:743561/1‑87 (MQ=255)
                                                                                                         ttttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGa                                            <  1:765636/150‑1 (MQ=39)
                                                                                                            tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACAt                                         <  2:89981/150‑1 (MQ=39)
                                                                                                                 gTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGa                                                                                                                                           <  1:1234917/47‑1 (MQ=255)
                                                                                                                  tCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTCCGCCAGCCAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTATGCCCCATCCGACATCAACGc                                   >  2:581841/1‑150 (MQ=255)
                                                                                                                    ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTGGTTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCt                                 >  1:34133/1‑150 (MQ=255)
                                                                                                                       tAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                      <  1:923263/62‑1 (MQ=255)
                                                                                                                       tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGCTGCGGCGATTGGAGGGTATTTTACGATCGACCCCAGACAGCGCCATAGTCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                              <  2:1340107/150‑1 (MQ=11)
                                                                                                                       tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                      >  2:923263/1‑62 (MQ=255)
                                                                                                                       tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGAt                              <  1:1181305/150‑1 (MQ=39)
                                                                                                                        aaGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCATTTATCCCGCATCCGGCGTTACCGCCTAATg                             >  1:1360550/1‑150 (MQ=11)
                                                                                                                         aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGc                            >  1:1058904/1‑150 (MQ=39)
                                                                                                                          ggCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCg                           >  2:909416/1‑150 (MQ=39)
                                                                                                                          ggCGTTTACGCCGCATCCGAACACAACGCCTGATGCGACGTTCGCCGCGTCTTATCAGGCCTACACCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGTTCGGGTCAGGCGTTTATGCCGCATCCGCCATGAAGGCCTGATGcc                           >  1:631695/1‑149 (MQ=11)
                                                                                                                            gctttCCGCCGCATCCGCCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGAc                         <  1:1072541/148‑1 (MQ=255)
                                                                                                                             gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACg                        >  1:746079/1‑150 (MQ=39)
                                                                                                                                     cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGCTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc                >  1:1317343/1‑150 (MQ=25)
                                                                                                                                        aTCCGACATCAGCGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGCTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTc             <  1:1262227/150‑1 (MQ=255)
                                                                                                                                        aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTc             >  2:845593/1‑150 (MQ=35)
                                                                                                                                        aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCAGACAGCGCAATAGCCTGATATAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGCTGCGACGCCTGAcgcgcc             >  2:845644/1‑148 (MQ=11)
                                                                                                                                         tCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                       <  1:210508/75‑1 (MQ=255)
                                                                                                                                           cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTa          <  1:883798/150‑1 (MQ=34)
                                                                                                                                              cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTCTGCCGCATCCGACATCAACGCCTGATGCGACGCTTTACGCGTCTTATCa       >  2:348934/1‑150 (MQ=25)
                                                                                                                                               aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAg      <  1:307625/150‑1 (MQ=34)
                                                                                                                                                 cAACGCCTGATGCGACCCTTGCCGCGTCTTATCAGGCTTTCGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc    <  2:245850/150‑1 (MQ=255)
                                                                                                                                                   aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGCCGCTTAACGCGTCTTATTAGGCCt  >  2:176876/1‑150 (MQ=14)
                                                                                                                                                    |                                                                                                                                                    
GTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913_3_bme_pgi/4295913‑4296209

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: