Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi4,296,0610CT33.3% 66.2 / 54.8 69intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (18/28);  new base T (11/12);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_bme_pgi/4295914‑4296207
                                                                                                                                                    |                                                                                                                                                  
tCGGATAAGACGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTATTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTCGTAGATCCGATAAGTCGTTTACGCCGCCTCTGACCTCAAcc                                                                                                                                                   >  2:387012/1‑148 (MQ=11)
 cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                  <  1:466151/149‑1 (MQ=17)
 cGGATAAGGCGTTTATGCCCCATCCGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                  <  1:303798/149‑1 (MQ=12)
   gttAAGGCGTTTATGCCGGTTCCGGCATCAATGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                <  2:556100/147‑1 (MQ=255)
   gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGGGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCTCACTAGCCTGGTTTAGC‑GTGTTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGccc                                                                                                                                                >  2:276885/1‑148 (MQ=11)
     tAATGCGTTTATGCCGCATCCTACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                              <  2:533346/149‑1 (MQ=17)
         gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCACCGCCTGATGCg                                                                                                                                          >  1:128477/1‑149 (MQ=17)
         gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                          >  1:808965/1‑149 (MQ=32)
           gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                            <  1:733143/145‑1 (MQ=17)
           gTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATg                                                                                                                                            >  2:733143/1‑145 (MQ=17)
                  ccGCATCCGACATCCACGCCTGATGCGCCGCTTAACGCGTCTTATCAGGCCTCCTCCAGACAGCTCAATAGCCTGTTTCGGC‑GTGATATTCTACATCGGACAAGCCGTTAACGCCGCATCAGACATAAATGCCTGGTGCGACCGTTGcc                                                                                                                                 >  2:118359/1‑149 (MQ=2)
                  ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGGCGGtttag                                                                                                                                 >  2:621372/1‑146 (MQ=12)
                  ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                 <  1:363611/149‑1 (MQ=34)
                   cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                >  2:94582/1‑149 (MQ=35)
                    gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                               <  2:808965/149‑1 (MQ=25)
                    ccATCCGCTATAAACGCCTGATGCGACGCTTGACGTGTCTTATCAGGCCTAAGCCAGACAGCGCAATATCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                               <  2:510704/148‑1 (MQ=11)
                       tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGCCAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATCCCTGATCCGACGCTTGCCGCGTc                                                                                                                            >  2:366532/1‑149 (MQ=255)
                          ggtattaaCGCCTGATGCGACGCTTGACGTGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                         <  2:628311/143‑1 (MQ=255)
                          gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                           >  1:131594/1‑131 (MQ=18)
                          gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCACAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                           <  2:131594/131‑1 (MQ=16)
                           aCATCAACGCCTGATGCGACGCTTAACGCGTTTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                        <  1:217973/149‑1 (MQ=25)
                             aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                      >  1:471843/1‑149 (MQ=37)
                             aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGAACTGTGCCAGTCATCGCAATCGCCTGATTGAGCTATGGTTTTG‑AGTTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTTCCGCGTCTTATct                                                                                                                      >  2:596555/1‑148 (MQ=11)
                              tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATccg                                                                                                                     >  2:206494/1‑147 (MQ=25)
                               gaaCGCCTGAGGCGGCGCTTAGCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAGAGGCGGATTTAGC‑GTGATTTTGTAGGTCGGATAGGGCGTTTACGCCGGAACAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                    <  2:96331/148‑1 (MQ=11)
                                 aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                  >  2:189590/1‑149 (MQ=25)
                                    ccTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTATAGGTCGGATAAGGCGTTTACGCCGCCTCCGACATCAACGCCTGATGGGACGCTTGGCGCGTCTTCTCAGGCCtgc                                                                                                               >  1:90794/1‑147 (MQ=11)
                                     cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                                                                                                                                >  1:249826/1‑115 (MQ=32)
                                     cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCtccc                                                                                                              >  2:363135/1‑146 (MQ=37)
                                       gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGAGNCNGCGCTTCCCGCGTCTTATCACGCCtccctt                                                                                                            >  2:270694/1‑144 (MQ=11)
                                        aTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           <  1:115247/149‑1 (MQ=39)
                                             tcgtttAACGCGTCTTATCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGataa                                                                                                      <  2:128477/145‑4 (MQ=255)
                                                   aaCGCGTCTTATCAGGCCTACGCCAGACAGCTCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                 <  2:670604/116‑1 (MQ=255)
                                                   aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                 >  1:254752/1‑116 (MQ=37)
                                                   aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                 <  2:254752/116‑1 (MQ=25)
                                                   aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                 >  1:670604/1‑116 (MQ=25)
                                                    ccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                             >  1:668872/2‑119 (MQ=17)
                                                    ccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGt                                                                                                                             <  2:668872/118‑1 (MQ=40)
                                                       cgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGcc                                                                                            <  1:477149/149‑1 (MQ=39)
                                                         tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                            >  2:157413/1‑99 (MQ=34)
                                                         tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                            <  1:157413/99‑1 (MQ=34)
                                                           ttATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGTTCGTATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAt                                                                                        <  1:443710/149‑1 (MQ=255)
                                                            tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGT‑GTGACTTCGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTGCGCCAGACAGCGCAATAGCCTGAtt                                                                                       <  1:635433/149‑1 (MQ=11)
                                                                aGGCCTACGCCGGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   <  1:94582/149‑1 (MQ=255)
                                                                    cTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   >  2:400355/1‑113 (MQ=12)
                                                                    cTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   <  1:400355/113‑1 (MQ=25)
                                                                           aGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTGCGCTGCAGAAGTCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                        <  1:367471/149‑1 (MQ=11)
                                                                           aGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   >  2:384603/1‑106 (MQ=17)
                                                                           aGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   <  1:384603/106‑1 (MQ=255)
                                                                           aGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                        >  1:595912/1‑149 (MQ=38)
                                                                                   cAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCAGAATGTTAGCAGGCCTACGACAGATAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGAt                                                                <  1:747087/149‑1 (MQ=11)
                                                                                         ccTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                  <  1:578618/61‑1 (MQ=255)
                                                                                         ccTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCATTAGCCTGATTTAGCGTGATTTTGTAGGTCGGCTAAGGCg                                                          >  1:314006/1‑149 (MQ=255)
                                                                                          cTGATTTAGC‑GTGATTTTGTAGGGCGGGTAGGGCGTTTACGCCGCATCCGATATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCTGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                                         <  2:729795/149‑1 (MQ=11)
                                                                                               ttggC‑GTGATTTTGTGGGTGGGATAAGGGGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCTTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATg                                                    <  1:734001/146‑1 (MQ=11)
                                                                                                 aGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGTCAACCAGCGGTATAGCTTGCTTTAGCTTGACTTCGTAGATCAGATACGGCGATTCTGct                                                  >  2:579171/1‑148 (MQ=11)
                                                                                                  gC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                   <  1:495939/83‑1 (MQ=255)
                                                                                                   c‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGAGACGCtt                                                                                                                                    <  2:518648/65‑1 (MQ=11)
                                                                                                     gTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACGGTGCATTAGCTTGAATTAGCGTGATTTTGTGGGTCGGACAAGTCGCTTTTGCCgta                                               >  2:53121/1‑147 (MQ=11)
                                                                                                              tAGGTTGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATc                                      <  1:595280/149‑1 (MQ=38)
                                                                                                               aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCa                                     <  1:276885/149‑1 (MQ=38)
                                                                                                                ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCaa                                    <  1:575600/149‑1 (MQ=37)
                                                                                                                 gTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGACAAGCCGTTTATGCCGCAACCGACATCAAc                                   >  2:78136/1‑149 (MQ=255)
                                                                                                                 gTCGGATAAGGCGTTTACGCCGCATACGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                   <  2:560241/149‑1 (MQ=255)
                                                                                                                           gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           <  1:23727/67‑1 (MQ=255)
                                                                                                                           gCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGt                                                                         <  1:77972/101‑1 (MQ=255)
                                                                                                                                     cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                   <  1:94258/81‑1 (MQ=40)
                                                                                                                                               aTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCTTCCGACATCAACGCCTGATGCGACGCTTTACGCGTCTTATCa     >  1:507418/1‑149 (MQ=21)
                                                                                                                                                  aaCGCCTTATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc  <  1:415701/149‑1 (MQ=17)
                                                                                                                                                  aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc  <  1:415680/149‑1 (MQ=32)
                                                                                                                                                    |                                                                                                                                                  
TCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGC‑GTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913_3_bme_pgi/4295914‑4296207

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: ‑