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

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC  >  NC_000913/1222083‑1222373
                                                                                                                                                       |                                                                                                                                             
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                               >  2:252417/1‑151 (MQ=255)
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                               <  1:805088/151‑1 (MQ=255)
 tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                               >  2:264789/1‑150 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                           >  2:976070/1‑150 (MQ=255)
      aaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAatat                                                                                                                                         <  2:804012/151‑1 (MQ=255)
          cGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑t                                                                                                                                    >  2:669550/1‑151 (MQ=255)
             aCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCt                                                                                                                                 >  2:201796/1‑151 (MQ=255)
               aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTg                                                                                                                                >  2:228720/1‑150 (MQ=255)
               aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGa                                                                                                                               <  1:817275/151‑1 (MQ=255)
                    ggggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCccg                                                                                                                          <  2:1116940/151‑1 (MQ=255)
                     gggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCccgc                                                                                                                         >  1:569899/1‑151 (MQ=255)
                        cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGt                                                                                                                      >  2:756790/1‑151 (MQ=255)
                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGc                                                                                                                  >  1:790216/1‑151 (MQ=255)
                             aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCt                                                                                                                 >  1:967030/1‑151 (MQ=255)
                               ttgggggaAAATTAAACTTGGCGTAAAGGGCACCCTGAACCCAGCGCCCAGCGTCGGGGGAAGGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGCCCCCGGAGTGTTGGGGGGCTTAAAATAAAAG‑TTCTGACCCCGCCGTTAGCt                                                                                                                 <  2:532148/146‑1 (MQ=255)
                                 tGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTaaa                                                                                                              >  1:241641/1‑150 (MQ=255)
                                   ggggAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTaaaaa                                                                                                            >  2:504840/1‑150 (MQ=255)
                                           tAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                   >  2:745157/1‑151 (MQ=255)
                                            aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                  >  1:435827/1‑151 (MQ=255)
                                            aaaCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                  >  2:785581/1‑151 (MQ=255)
                                               cTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                               >  2:1066435/1‑151 (MQ=255)
                                                 tGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                              >  2:1079649/1‑150 (MQ=255)
                                                      tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                        >  1:1056713/1‑151 (MQ=255)
                                                       aaaTGGCACCCTAAACCCAGCGCCCAGCGTCTGGGGAAGGGGGGCGTGCGGGGGGG‑TAATAATGGCTACAATGACCCCGCGGTGATGGTGGGCCTAAAATATAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                           <  2:701855/147‑1 (MQ=255)
                                                             cacccTGAACCCAGCGGCCGGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TAATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATAAAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                  <  1:205806/147‑1 (MQ=255)
                                                             cAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAGTATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                  >  2:52802/1‑150 (MQ=255)
                                                               aCCTGAATCCAGCAGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGCTGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGCACat                                                                               >  2:1191474/1‑151 (MQ=255)
                                                                ccTGACTCCGGCGCCCGGCGTCGGGGGAAGGGGGGGGTGCAGGGGGG‑TGATAATGGCTAAAATCCCCCCGAGGTGATGGGGGCCCTAAAATATAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                    <  2:299886/145‑1 (MQ=255)
                                                                ccTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                               >  1:907669/1‑150 (MQ=255)
                                                                ccTGAACCCAGCGCCCGGGTTCGGGGGAATGTGGGCGTGCAGCTGGG‑TAATAATGGCAAAAATGCCCCCGCGGTGATGGGGGCCTTAAAATATAAT‑TTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                   <  2:1091268/146‑1 (MQ=255)
                                                                  ggagcccggCGGCCAGCGTCTGGGAAAGGTGGGGGTGCGGCGGGG‑GAATATTGGCTACAAGGCCGCCGGAGTG‑TGGGGGGGTTTAAAAATTAGTTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                             <  1:696718/142‑1 (MQ=255)
                                                                    aaTCCAGCGGCCGGCGTCGGGGGAAGGTGGGGGTGCAGCGGGG‑TGATAATGGCTAAAATGACCCCGCGGTGATGGGGGGCCTAAAATAAAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                          <  1:740455/151‑1 (MQ=255)
                                                                    aaTCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                           >  2:37208/1‑150 (MQ=255)
                                                                     acccggcgcCCGGCGTCGGGGAAGTGGGGGCGTGCGGCGGGG‑TAATATTGCCCACATTCACTCGGCGGTGGTGGTGGCCCTAAAATAAACG‑TTATGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                           <  1:951269/144‑1 (MQ=255)
                                                                       ccAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                        >  2:129648/1‑150 (MQ=255)
                                                                           cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                   >  2:81533/1‑151 (MQ=255)
                                                                            ggCCGGGGTCGGGGGAAGGGGGGCGTGCAGGTGGG‑TAAAAATGGCTAAAAGGACCCCGCAGTGAGGGGGGGCTAAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                    <  2:522177/149‑1 (MQ=255)
                                                                            ggCCAGGGTTGGGGGAAGGGGGGCTTGCAGCGGGG‑TAATAATGGCTAAAATGACCCCGCAGTATTGGTGGCCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                    <  2:362647/149‑1 (MQ=255)
                                                                            ggCCAGCGTGTGGGGAAGGGGGGGGTGCAGGTGGG‑TAATAATGGCTAAAATGCCCCCGCGGTGAGGGGGGGCCAAAAATATAAT‑TTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatta                                                                    <  2:352444/149‑1 (MQ=255)
                                                                               gcgggttggggaAATGGGGGCGTGGAGCGGGG‑AAAAATTGGCTACAGTGCCACGGCGGTGGTGGTGGGCCTAAAATAAAAG‑TTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                <  1:1060316/142‑1 (MQ=255)
                                                                               cAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                               >  2:197028/1‑151 (MQ=255)
                                                                                 gCGTCGGGGGAAGGGGGGCGTGCAGGGGGG‑TGATAATGGCAACAACGACCCCGCAGTGATGGGGGGCCTAAAATAAAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                               <  1:126156/149‑1 (MQ=255)
                                                                                 gCGTCGGGGAAGTGTGGGGGT‑CAGGTGGGTTGATAAGGAATACAATACCCCCGGAGGGATGGTGGCCTAAAAAAAAAAT‑TTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                              <  2:798838/150‑1 (MQ=255)
                                                                                   gtggggggAATGTGGGCTTGCAGTTGGG‑TAATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                            <  2:50243/146‑1 (MQ=255)
                                                                                   gTCGGGGGAATGGGGGGGGGCAGCGGGG‑TAATAATGGCTAAAATGCCCCCGAAGTGATGGTGGGCTTAAAATAAAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                               <  1:1022902/147‑1 (MQ=255)
                                                                                     cTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                          >  1:191396/1‑150 (MQ=255)
                                                                                       ggggAAGGGGGGGGTGCAGGTGGG‑TAATAATGGCTACAATGACCCCGAAGGGTTGGGGGCCCTAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                            <  2:113637/146‑1 (MQ=255)
                                                                                       ggggAAGGGGGGCGTGCGGCGGGG‑TAATAAGGTTAACAAGGCCCCCGTGGTGATGGTGGGCAAAAAATAAAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                       <  2:1102615/151‑1 (MQ=255)
                                                                                        gggAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                        >  1:1072787/1‑149 (MQ=255)
                                                                                           aaTGTGGGCGTGCAGGTGGG‑TAATAATGGCTACAATGACACCGCAGTGAGGGGGGCCTTAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                    <  2:292767/150‑1 (MQ=255)
                                                                                               tgGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                               >  2:1115340/1‑151 (MQ=255)
                                                                                                gggCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATg                                              >  1:713305/1‑151 (MQ=255)
                                                                                                gggCGTGCAGCGGGG‑TAATAATGACTACAATGCCCCCGCGGTGAGGGTGGCCCAAAAATATATG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                               <  2:1111980/150‑1 (MQ=255)
                                                                                                 ggcgtggggggggg‑TAATAATGGCAACAATGCCCCCGGAGGGAGGGGGGGCCAAAAAAATATG‑TTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                             <  1:206226/140‑1 (MQ=255)
                                                                                                 ggCGTGCAGCGGGG‑TAATAATGGCTAAAAGAACACCGCGGTGATGGGGGGCTTAAAATATAGG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                             <  1:680320/151‑1 (MQ=255)
                                                                                                    gTGCGGCTGGG‑TGATAATGGCAACAATGACACCGCGGTGATGGTGGGCCTAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                           <  1:71788/150‑1 (MQ=255)
                                                                                                            gag‑TGATAAGGGCTAAAAGGACGCCGGAGTGGGGGGGGGCAAAAAATATAAG‑TTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtca                                        <  2:78170/143‑1 (MQ=255)
                                                                                                                    aagggCTAAAATGACACCGGGGTGATGGGGGGCCTAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                           <  1:505521/148‑1 (MQ=255)
                                                                                                                    aaTGCCTACAAGCACGCCGCAGTGAGGGTGGGCGTGAAATATAGG‑TTTTCGTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                            <  2:349132/150‑1 (MQ=255)
                                                                                                                           aCAATGACACCGCAGTGATGGTGGGCCTAAAATATAAG‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                    <  1:431531/151‑1 (MQ=255)
                                                                                                                            cAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                   >  2:79428/1‑151 (MQ=255)
                                                                                                                             aaTGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                   >  2:414485/1‑150 (MQ=255)
                                                                                                                               tcacacCGCAGTGATGGGGGGCCTAAAATATAAG‑TTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                 <  2:398478/148‑1 (MQ=255)
                                                                                                                                  ccccGCAGTGATGGTGGGCGTAAAATATAAG‑TTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg              <  1:641273/148‑1 (MQ=255)
                                                                                                                                    ccGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg           >  2:1164999/1‑151 (MQ=255)
                                                                                                                                    ccGCAGTGAGGGTGGGCCTAAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg            <  1:1096828/150‑1 (MQ=255)
                                                                                                                                        aGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt       >  2:990104/1‑151 (MQ=255)
                                                                                                                                             tggtggGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc  >  2:70564/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                             
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGG‑TGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGC  >  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: ‑