Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA42.2% 87.2 / 25.8 65pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (32/5);  new base A (0/27);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.76e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC  >  NC_000913/1222083‑1222373
                                                                                                                                                       |                                                                                                                                             
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                               >  2:252417/1‑151 (MQ=255)
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                               <  1:805088/151‑1 (MQ=255)
 tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                               >  2:264789/1‑150 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                           >  2:976070/1‑150 (MQ=255)
      aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat                                                                                                                                         <  2:804012/151‑1 (MQ=255)
          cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGt                                                                                                                                    >  2:669550/1‑151 (MQ=255)
             aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCt                                                                                                                                 >  2:201796/1‑151 (MQ=255)
               aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTg                                                                                                                                >  2:228720/1‑150 (MQ=255)
               aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                               <  1:817275/151‑1 (MQ=255)
                    ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccg                                                                                                                          <  2:1116940/151‑1 (MQ=255)
                     gggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCccgc                                                                                                                         >  1:569899/1‑151 (MQ=255)
                        cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGt                                                                                                                      >  2:756790/1‑151 (MQ=255)
                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                  >  1:790216/1‑151 (MQ=255)
                             aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCt                                                                                                                 >  1:967030/1‑151 (MQ=255)
                               ttgggggaAAATTAAACTTGGCGTAAAGGGCACCCTGAACCCAGCGCCCAGCGTCGGGGGAAGGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGCCCCCGGAGTGTTGGGGGGCTTAAAATAAAAGTTCTGACCCCGCCGTTAGCt                                                                                                                 <  2:532148/146‑1 (MQ=255)
                                 tGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaa                                                                                                              >  1:241641/1‑150 (MQ=255)
                                   ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTaaaaa                                                                                                            >  2:504840/1‑150 (MQ=255)
                                           tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                   >  2:745157/1‑151 (MQ=255)
                                            aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                  >  1:435827/1‑151 (MQ=255)
                                            aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                  >  2:785581/1‑151 (MQ=255)
                                               cTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                               >  2:1066435/1‑151 (MQ=255)
                                                 tGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                              >  2:1079649/1‑150 (MQ=255)
                                                      tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                        >  1:1056713/1‑151 (MQ=255)
                                                       aaaTGGCACCCTAAACCCAGCGCCCAGCGTCTGGGGAAGGGGGGCGTGCGGGGGGGTAATAATGGCTACAATGACCCCGCGGTGATGGTGGGCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                           <  2:701855/147‑1 (MQ=255)
                                                             cacccTGAACCCAGCGGCCGGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTAATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                  <  1:205806/147‑1 (MQ=255)
                                                             cAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAGTATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                  >  2:52802/1‑150 (MQ=255)
                                                               aCCTGAATCCAGCAGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGCTGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGCACat                                                                               >  2:1191474/1‑151 (MQ=255)
                                                                ccTGACTCCGGCGCCCGGCGTCGGGGGAAGGGGGGGGTGCAGGGGGGTGATAATGGCTAAAATCCCCCCGAGGTGATGGGGGCCCTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                    <  2:299886/145‑1 (MQ=255)
                                                                ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                               >  1:907669/1‑150 (MQ=255)
                                                                ccTGAACCCAGCGCCCGGGTTCGGGGGAATGTGGGCGTGCAGCTGGGTAATAATGGCAAAAATGCCCCCGCGGTGATGGGGGCCTTAAAATATAATTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                   <  2:1091268/146‑1 (MQ=255)
                                                                  ggagcccggCGGCCAGCGTCTGGGAAAGGTGGGGGTGCGGCGGGGGAATATTGGCTACAAGGCCGCCGGAGTGTGGGGGGGTTTAAAAATTAGTTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                             <  1:696718/142‑1 (MQ=255)
                                                                    aaTCCAGCGGCCGGCGTCGGGGGAAGGTGGGGGTGCAGCGGGGTGATAATGGCTAAAATGACCCCGCGGTGATGGGGGGCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                          <  1:740455/151‑1 (MQ=255)
                                                                    aaTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                           >  2:37208/1‑150 (MQ=255)
                                                                     acccggcgcCCGGCGTCGGGGAAGTGGGGGCGTGCGGCGGGGTAATATTGCCCACATTCACTCGGCGGTGGTGGTGGCCCTAAAATAAACGTTATGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                           <  1:951269/144‑1 (MQ=255)
                                                                       ccAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                        >  2:129648/1‑150 (MQ=255)
                                                                           cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                   >  2:81533/1‑151 (MQ=255)
                                                                            ggCCGGGGTCGGGGGAAGGGGGGCGTGCAGGTGGGTAAAAATGGCTAAAAGGACCCCGCAGTGAGGGGGGGCTAAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                    <  2:522177/149‑1 (MQ=255)
                                                                            ggCCAGGGTTGGGGGAAGGGGGGCTTGCAGCGGGGTAATAATGGCTAAAATGACCCCGCAGTATTGGTGGCCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                    <  2:362647/149‑1 (MQ=255)
                                                                            ggCCAGCGTGTGGGGAAGGGGGGGGTGCAGGTGGGTAATAATGGCTAAAATGCCCCCGCGGTGAGGGGGGGCCAAAAATATAATTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                    <  2:352444/149‑1 (MQ=255)
                                                                               gcgggttggggaAATGGGGGCGTGGAGCGGGGAAAAATTGGCTACAGTGCCACGGCGGTGGTGGTGGGCCTAAAATAAAAGTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                <  1:1060316/142‑1 (MQ=255)
                                                                               cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                               >  2:197028/1‑151 (MQ=255)
                                                                                 gCGTCGGGGGAAGGGGGGCGTGCAGGGGGGTGATAATGGCAACAACGACCCCGCAGTGATGGGGGGCCTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                               <  1:126156/149‑1 (MQ=255)
                                                                                 gCGTCGGGGAAGTGTGGGGGTCAGGTGGGTTGATAAGGAATACAATACCCCCGGAGGGATGGTGGCCTAAAAAAAAAATTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                              <  2:798838/150‑1 (MQ=255)
                                                                                   gtggggggAATGTGGGCTTGCAGTTGGGTAATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                            <  2:50243/146‑1 (MQ=255)
                                                                                   gTCGGGGGAATGGGGGGGGGCAGCGGGGTAATAATGGCTAAAATGCCCCCGAAGTGATGGTGGGCTTAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                               <  1:1022902/147‑1 (MQ=255)
                                                                                     cTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                          >  1:191396/1‑150 (MQ=255)
                                                                                       ggggAAGGGGGGGGTGCAGGTGGGTAATAATGGCTACAATGACCCCGAAGGGTTGGGGGCCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                            <  2:113637/146‑1 (MQ=255)
                                                                                       ggggAAGGGGGGCGTGCGGCGGGGTAATAAGGTTAACAAGGCCCCCGTGGTGATGGTGGGCAAAAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                       <  2:1102615/151‑1 (MQ=255)
                                                                                        gggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                        >  1:1072787/1‑149 (MQ=255)
                                                                                           aaTGTGGGCGTGCAGGTGGGTAATAATGGCTACAATGACACCGCAGTGAGGGGGGCCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                    <  2:292767/150‑1 (MQ=255)
                                                                                               tgGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                               >  2:1115340/1‑151 (MQ=255)
                                                                                                gggCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg                                              >  1:713305/1‑151 (MQ=255)
                                                                                                gggCGTGCAGCGGGGTAATAATGACTACAATGCCCCCGCGGTGAGGGTGGCCCAAAAATATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                               <  2:1111980/150‑1 (MQ=255)
                                                                                                 ggcgtgggggggggTAATAATGGCAACAATGCCCCCGGAGGGAGGGGGGGCCAAAAAAATATGTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                             <  1:206226/140‑1 (MQ=255)
                                                                                                 ggCGTGCAGCGGGGTAATAATGGCTAAAAGAACACCGCGGTGATGGGGGGCTTAAAATATAGGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                             <  1:680320/151‑1 (MQ=255)
                                                                                                    gTGCGGCTGGGTGATAATGGCAACAATGACACCGCGGTGATGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                           <  1:71788/150‑1 (MQ=255)
                                                                                                            gagTGATAAGGGCTAAAAGGACGCCGGAGTGGGGGGGGGCAAAAAATATAAGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca                                        <  2:78170/143‑1 (MQ=255)
                                                                                                                    aagggCTAAAATGACACCGGGGTGATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                           <  1:505521/148‑1 (MQ=255)
                                                                                                                    aaTGCCTACAAGCACGCCGCAGTGAGGGTGGGCGTGAAATATAGGTTTTCGTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                            <  2:349132/150‑1 (MQ=255)
                                                                                                                           aCAATGACACCGCAGTGATGGTGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                    <  1:431531/151‑1 (MQ=255)
                                                                                                                            cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                   >  2:79428/1‑151 (MQ=255)
                                                                                                                             aaTGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                   >  2:414485/1‑150 (MQ=255)
                                                                                                                               tcacacCGCAGTGATGGGGGGCCTAAAATATAAGTTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                 <  2:398478/148‑1 (MQ=255)
                                                                                                                                  ccccGCAGTGATGGTGGGCGTAAAATATAAGTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg              <  1:641273/148‑1 (MQ=255)
                                                                                                                                    ccGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg           >  2:1164999/1‑151 (MQ=255)
                                                                                                                                    ccGCAGTGAGGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg            <  1:1096828/150‑1 (MQ=255)
                                                                                                                                        aGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt       >  2:990104/1‑151 (MQ=255)
                                                                                                                                             tggtggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc  >  2:70564/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC  >  NC_000913/1222083‑1222373

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