Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT43.5% 30.1 / 104.5 87intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base T (20/17);  total (49/38)
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.14e-03
Rejected as consensus: Frequency below/above cutoff threshold.

CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913/4295914‑4296206
                                                                                                                                                  |                                                                                                                                                  
cGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAAc                                                                                                                                                    >  2:576795/1‑147 (MQ=17)
    tAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATATTGTAGGTCGGATAAGTCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                >  1:908914/1‑147 (MQ=12)
     caGGCGTTGATGTCGGATGCGGCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                               >  1:556593/2‑147 (MQ=11)
     atggCGTTTATGCCGCTTAAGTCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                               <  2:745956/145‑1 (MQ=255)
     aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                               >  2:588052/1‑147 (MQ=17)
      aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                              <  1:1127699/147‑1 (MQ=17)
         cGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCAGCATCCGACATCAACGCCTGATGc                                                                                                                                           >  1:721027/1‑147 (MQ=17)
            ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGAc                                                                                                                                        >  1:703052/1‑147 (MQ=17)
            ttATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCAAGACATCGCAATAGCCTGATTAAGCGTGATTTTGTAGGTCGGATAAGGCTGTTACGACGCATCCGACATCAAAGCCTGATGCGAc                                                                                                                                        >  2:253808/1‑147 (MQ=11)
             tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACg                                                                                                                                       >  1:720510/1‑147 (MQ=17)
                 ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                                   >  1:579512/1‑147 (MQ=17)
                   gCATCAGTCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                 <  2:1091380/147‑1 (MQ=21)
                    cATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCAAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAGATCAATGACTGATGCGACGCTTggcc                                                                                                                                >  2:167005/1‑144 (MQ=11)
                       ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGt                                                                                                                             >  1:863008/1‑147 (MQ=25)
                          aCATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGAGCTATAGCCTTATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCTTCCGTCATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          <  2:303792/147‑1 (MQ=11)
                          aCATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGGGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          <  1:28512/147‑1 (MQ=255)
                          aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          >  1:406486/1‑147 (MQ=25)
                          aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                          >  1:1162513/1‑147 (MQ=37)
                           cATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                         <  1:1292965/147‑1 (MQ=255)
                           cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCttt                                                                                                                         >  2:41467/1‑146 (MQ=25)
                           cATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGACAGACAGCGCAATAGCGTGAATTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATACGACATAAAATCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                         >  2:793352/1‑147 (MQ=11)
                            tgcatcGCCTGATCTGACGCTTGACGCGTCTTATGAGGCGTATGCCTGAGAGAGAAATCGCATGATTTTGCTTGATTTATTATTTCTGTTTTGGGGTTTACTCCGCTTCCGTCATGAACGGCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                        <  2:192572/142‑1 (MQ=2)
                            ataaaCGCCTGATGCGACGATTTACGCGTCTTATCAGGCCTACGGGAGTCAGCTCAATTTGCTGATTTAGCGTGATTTTGTTTTTCGGTTAAGGCTTTTACTCCGCATCAGTCTTCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                        <  2:192571/144‑1 (MQ=2)
                             tCAACGCCTGATGCGACGCTTAACTCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGTATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       <  2:455233/147‑1 (MQ=255)
                             tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGAGAGCTAAATAGCCTGATTTAGCTTGATTTTGTTGGTCGGATAAGGCGTTTACGCCGCTTCCGACTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       <  1:534065/147‑1 (MQ=11)
                             tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       <  2:1091042/147‑1 (MQ=37)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                      >  1:597419/1‑147 (MQ=25)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                      >  2:591470/1‑147 (MQ=25)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCa                                                                                                                      <  1:321478/147‑1 (MQ=25)
                                 cGCCTGATGCGACGCTTAACGCGTCTTATCCGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCTTCTTATCAgtc                                                                                                                   >  1:1167815/1‑145 (MQ=255)
                                 cGCCTGATGCGACGCTTAACGCGTCTTATAAGGCCTACGCCAGACAGCGCAATCGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGAGATTACTCCTCATACGACATCAATTCCTGATTCGACGCTTGAAGCGTCTTATCAGGc                                                                                                                   >  2:708965/1‑147 (MQ=11)
                                  gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATACGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                  >  1:404911/1‑147 (MQ=25)
                                  gCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGACGCCTCTTAACAGGcc                                                                                                                  >  2:955894/1‑147 (MQ=255)
                                   ccTGATGCGACGCTTCACGCGTCTTATCAGGCCTACTCAATACAGCTCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTTCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                 <  2:716668/147‑1 (MQ=11)
                                    cTGATGCTACGATTAAGGCGTCTTATCAGGCGTACGCCTGTCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  1:1238444/147‑1 (MQ=11)
                                    cTGATGCGACGCTTACCGCGTCTTAGCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  1:793352/147‑1 (MQ=255)
                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  2:326129/147‑1 (MQ=25)
                                    cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                <  2:568091/147‑1 (MQ=25)
                                      gataagACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTGTTATCAGGCCTACg                                                                                                              >  1:1215671/6‑147 (MQ=11)
                                      gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                              >  1:903729/1‑145 (MQ=18)
                                      gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       <  2:597755/138‑1 (MQ=37)
                                      gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                       >  1:597755/1‑138 (MQ=37)
                                         gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatctta                                                                                                           <  1:548736/147‑6 (MQ=9)
                                               cTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGACAGACAgc                                                                                                     >  1:1179425/1‑147 (MQ=37)
                                                         cTTATCAGGCCTACGCCAGACAGATCAATTTAATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCTTTTACGCCGCATAAGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCt                                                                                           <  2:7831/147‑1 (MQ=11)
                                                           tATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGa                                                                                         <  1:507971/147‑1 (MQ=37)
                                                              cAGGCCTACGCCAGAATGCTCTATATCCTTATTTCGCGTGATTATGTATTGCTGATATGGGTTTTACGCCGGATAAGATTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGAttt                                                                                      <  2:730795/147‑1 (MQ=2)
                                                               aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                  >  1:820544/1‑86 (MQ=18)
                                                               aGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGc                                                                                                                                                  <  2:820544/86‑1 (MQ=18)
                                                                   cTACGCCAGTCAGCTCAATATCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGt                                                                                 <  1:394241/147‑1 (MQ=255)
                                                                             cAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                           <  1:668266/111‑1 (MQ=39)
                                                                                   aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                                                                                                     <  2:610731/79‑1 (MQ=21)
                                                                                   aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCt                                                                                                                                     >  1:610731/1‑79 (MQ=34)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                           >  2:819856/1‑66 (MQ=25)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGc                                                                                                                                           <  1:819856/66‑1 (MQ=25)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                                   <  1:1058632/74‑1 (MQ=25)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                                   <  1:1274151/74‑1 (MQ=25)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                   >  2:1274151/1‑74 (MQ=35)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTg                                                                                                                                   >  2:1058632/1‑74 (MQ=35)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGTCg                                                          >  1:1279793/1‑147 (MQ=37)
                                                                                          tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                          >  2:1098373/1‑147 (MQ=37)
                                                                                            aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCAAGACATAGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGtt                                                        >  1:943462/1‑147 (MQ=255)
                                                                                            aTTTAGCGTGATTTTGTAGGTCGGATAAGGAGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGAGTGATTTTGTAGGTCGGATAAGGAGtt                                                        >  2:438771/1‑147 (MQ=255)
                                                                                              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                >  2:892720/1‑57 (MQ=255)
                                                                                              ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                <  1:892720/57‑1 (MQ=255)
                                                                                               tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                              >  2:454335/1‑58 (MQ=255)
                                                                                               tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGa                                                                                                                                              <  1:454335/58‑1 (MQ=255)
                                                                                                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAtt                                                                             <  1:32125/117‑1 (MQ=39)
                                                                                                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGAGCAATAGCCTGAATTAGCGTGATATTGTAGGTCGGATAAGGCGTTTATGACGCa                                               >  2:548736/1‑147 (MQ=255)
                                                                                                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                    <  2:619134/70‑1 (MQ=255)
                                                                                                             aGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                    >  1:619134/1‑70 (MQ=255)
                                                                                                               gTCGGATGCGGCGTAAACGCCTTATCCGACCTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCTATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATACggaagac                                     >  1:800066/1‑141 (MQ=11)
                                                                                                                  ggATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACg                                  >  1:449365/1‑147 (MQ=37)
                                                                                                                     tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCt                               >  1:904545/1‑147 (MQ=37)
                                                                                                                     tAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCAGACATCAACGCCt                               >  2:1170750/1‑147 (MQ=255)
                                                                                                                       aGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                                                                                                              <  1:265625/50‑1 (MQ=255)
                                                                                                                       aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcgcg                                                                                                                              >  2:265625/1‑50 (MQ=255)
                                                                                                                          cGTTTACGCCGCATCCGACATCAACGCCTGATGCTATGCTTTCCTCGTCTTTTCAGGCCTACGCCATACAGCGCAATAGCCTGATTTAGCTTGATTTTGTATGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGc                          <  2:721027/147‑1 (MQ=11)
                                                                                                                                  ccgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGcgtt                  <  2:1162513/147‑4 (MQ=255)
                                                                                                                                   cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                           <  2:1037200/105‑1 (MQ=35)
                                                                                                                                   cgcATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                           >  1:1037200/1‑105 (MQ=35)
                                                                                                                                   cgcATCAGACATCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                 <  2:1031990/147‑1 (MQ=12)
                                                                                                                                           aCATCAACGCCTGATGCGACGCTTTCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATaa                                                              >  2:150781/1‑94 (MQ=18)
                                                                                                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGACTGATGCGACGCTTAACGCGTCTTa        >  2:959942/1‑147 (MQ=21)
                                                                                                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACTCGTCTTa        >  2:970501/1‑147 (MQ=21)
                                                                                                                                            cATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTAGTAGGTCGGATAAGGCGTTTATGACGCATCCGACATCAAAGCCTGATGCGACGCTAAACCCGTCTTa        >  1:1270149/1‑147 (MQ=11)
                                                                                                                                              tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGAAAGACAGAGCAATATAATGATATATCGTGATTTTGTAGGTAGGATAATGCGTTTAAGCCGCATCCGACATCAAATCCTGATGCGACGCTTAACGCGTCTTAta      >  2:167753/1‑146 (MQ=11)
                                                                                                                                               cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCa     >  1:320573/1‑147 (MQ=32)
                                                                                                                                                  cGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc  >  1:716668/1‑147 (MQ=11)
                                                                                                                                                  |                                                                                                                                                  
CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGC  >  NC_000913/4295914‑4296206

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: