Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,222,2330GA42.9% 111.8 / 50.1 91pseudogene (60/75 nt)ycgHpseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (40/12);  new base A (0/39);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 1.96e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCT  >  NC_000913/1222083‑1222381
                                                                                                                                                         |                                                                                                                                                      
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                                        <  2:992385/151‑1 (MQ=255)
gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                                        <  2:185910/151‑1 (MQ=255)
 tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACGTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                                        <  1:470059/150‑1 (MQ=255)
  gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaa                                                                                                                                                      <  2:387385/151‑1 (MQ=255)
  gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGa                                                                                                                                                       >  2:952141/1‑150 (MQ=255)
  gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGa                                                                                                                                                       <  1:1316805/150‑1 (MQ=255)
  gtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGa                                                                                                                                                       <  1:1371857/150‑1 (MQ=255)
    gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGaaa                                                                                                                                                     <  1:626978/150‑1 (MQ=255)
    gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    >  2:905312/1‑151 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    >  1:1062560/1‑150 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    >  2:676151/1‑150 (MQ=255)
     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAat                                                                                                                                                    <  1:116755/150‑1 (MQ=255)
           gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑tt                                                                                                                                           <  2:1474773/151‑1 (MQ=255)
            aaCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTc                                                                                                                                          <  1:1071869/151‑1 (MQ=255)
                aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGAt                                                                                                                                      <  2:429713/151‑1 (MQ=255)
                  aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATcc                                                                                                                                    <  2:741189/151‑1 (MQ=255)
                        cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCccgccg                                                                                                                               >  2:225041/1‑150 (MQ=255)
                           gTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGtt                                                                                                                             >  1:1188660/1‑149 (MQ=255)
                           gTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGtt                                                                                                                             >  2:573404/1‑149 (MQ=255)
                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGc                                                                                                                          >  2:1445373/1‑151 (MQ=255)
                             aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCt                                                                                                                         >  1:1043872/1‑151 (MQ=255)
                             aaTTTGGGGCAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCt                                                                                                                         >  1:933178/1‑151 (MQ=255)
                                  gggggAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGAGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTaaaa                                                                                                                     >  2:792063/1‑150 (MQ=255)
                                  gggaaaaaTAAAATTTGGAGAAAATGCAAA‑CGGAAACCCAGCGGCCAGGTGGTGGGGAA‑TGGGGGCGGGCAC‑CTGGGGGATAATGCCTCCAACAACACCGAAGTGAAGGGGGGCATAAAAT‑ATAGG‑TTCGTATCCCGCCGTTAGCTaaaaa                                                                                                                    <  1:975114/146‑1 (MQ=255)
                                    gggAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAAc                                                                                                                  >  2:727243/1‑151 (MQ=255)
                                       aaaTTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACcgc                                                                                                               >  2:1174510/1‑151 (MQ=255)
                                        aaTAAAACTGGGCGAAAAGGGCAC‑CCTAAACCCAGCGGCCGGCGTCGGGGGAA‑TGTGGGCGTGCGG‑TTGGGTGATAATGGCTAAAATGCCACCGCGGTGATGGGGGCCCTAAAAT‑ATAAG‑TTCTGACCCCGCCGTTAGCTAAAAAACcgcg                                                                                                              <  2:288848/151‑1 (MQ=255)
                                          ttAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGcgtc                                                                                                            >  2:2091358/1‑151 (MQ=255)
                                           tAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                           >  1:48417/1‑151 (MQ=255)
                                             aaCTTGGCGTAAAGGGCAC‑CCTGAACCCAGCGGCGTGCGTGGGGGGAG‑TGGGGGCGTGCGG‑GGGGGTAAAAAGGCCAACAAGAACCCCGCAGTGATGGGGGGCATAAAAA‑AAAAT‑TTTGAAGCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                         <  1:1764717/151‑1 (MQ=255)
                                                 gggCGTAAATGGCAC‑CCTGAATCCAGCGGCCAGCGTCGGGGGAG‑TGGGGGCGTGCTG‑TGGGGTGATAATGGCAAAAATGCCCCCCCAGTGTTGGTGGGCCTAAAAT‑ATATG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                       <  2:793460/148‑1 (MQ=255)
                                                   gggTAAATGGCAC‑CCTAAACCCGGCGCCCAGCGTCGGGGGAA‑TGGGGGCGTGCAG‑CTGGGTGATAATGGTTACAATGCCCCCGCAGTGATGGGGGGCCGAAAAT‑ATATG‑TTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                   <  2:830040/149‑1 (MQ=255)
                                                   ggcgaaaaGGCAACCCTAAACCCGGCGGCCAGCTTCTGGGGAT‑TGGGGGGGTGCGG‑TTGGGTAAAAAGGGCTAAAATGACCCCGCAGTGAGGGTGGGCCTAAAAT‑ATAAT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                                    <  2:1737283/147‑1 (MQ=255)
                                                   gCGAAAAGGCCAC‑CCTAAACCCAGGGGCCAGCGTCGGGGGAA‑TGGGGGCGTGCAG‑CTGGGTGATAAGGGCTACAATGCCCCCGCAGTGAGGGGGGGCCTAAAAA‑ATAAT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                      <  2:981264/148‑1 (MQ=255)
                                                     gTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                  >  2:631925/1‑150 (MQ=255)
                                                     gTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                  >  1:988557/1‑150 (MQ=255)
                                                       aaaGGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTAATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTAAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                  <  1:198876/148‑1 (MQ=255)
                                                        aaTGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                               >  1:1141587/1‑150 (MQ=255)
                                                         agggCAA‑CCTGAATCCGGGGGCCAGTTTCGGGGGAA‑GGGGGGCGTGCGG‑CTGGTTTATATTGCCTCAACTGACCCCGGGGTGATGGGGGCCATAAAAA‑GTAAG‑TCCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgcg                                                                                             <  1:632761/149‑1 (MQ=255)
                                                         aatgacacccgaaacccAGCGGCCGGCGTGGGGGGAT‑GGGGGGCGTGCAG‑GGGGGTAAAAAGGGCTACAAGGCCCCGGAGGGGGGGGGGGGCCTAAAAT‑ATAAG‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                              <  1:2078184/136‑1 (MQ=255)
                                                               a‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                        >  1:903434/1‑150 (MQ=255)
                                                                 cctaaacccgggggCCGGGGTCGGGGGAG‑GGGGGGCGGGCGG‑GGGGGTGATAAGGGCTAAAATGACCCCGCAGTGTGGGGGGGCTTAAAAT‑AAAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                        <  2:1157705/137‑1 (MQ=255)
                                                                 ccgaaacccAGCGCCCGGCGTTGGGGGAA‑GGGGGGCGTGCAG‑GGGGGTGATAAGGGCTACAATCCCCCCGAGGTGAGGGGGGGCCTAAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                       <  1:331124/143‑1 (MQ=255)
                                                                 ccTGAACCCGGCGGCCAGGTTCGGGGAAA‑TGGGGGCGTGCAG‑GGGGGTGATAATGGCTACAATAACACCGCAGTGAGGGTGGGCTTAAAAT‑ATAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                      <  2:556329/151‑1 (MQ=255)
                                                                     accccAGCGGCCAGCGTCTGGGGAT‑TGGGGGCGTGCGG‑CTGGGTAATAATGGCTACAATGACCCCGCGGTGGTGGGGGCCCTAAAAT‑AAAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                                   <  2:52405/147‑1 (MQ=255)
                                                                     aagcgaggggCCGGCGT‑TGGGGAAGTGTGGGGGGGCAG‑GTGGGAGATAATGGCTAAAATGACCCCGCAGTGTGGGGGCGCAAAAAAA‑ATAAT‑TTTTGTGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                     <  2:1530797/140‑1 (MQ=255)
                                                                        cgaGCGGCCGGCGTCGGGGGAA‑TGTGGGCGGACAG‑GGGGGTAAAAATGCTTAAAACCACCCCGGAGTGATGGTGGGCTTAAAAT‑TTAAG‑TTTTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                     <  2:33118/143‑1 (MQ=255)
                                                                        ccAGGGGCCGGCGTCGGGGGAA‑TGGGGGGGGGCGG‑TGGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGCCCTAAAAT‑ATAAG‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                               <  2:1163905/151‑1 (MQ=255)
                                                                        ccAGCGCCCAGCGTGGGGGGAA‑TGTGGGCGGGCAG‑GTGGGTAAAAATTGCTACAATGCCCCCGAAGTTA‑GGGGGGTCAAAAATAATATT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                               <  1:135556/151‑1 (MQ=255)
                                                                         cggCGGCCAGCGTGGGGGGAG‑TGGGGGGGGGCAG‑GTGGGTGATAGGGGCTACAATGCCCCCGCGGTGATGGGGGGCTAAAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                <  1:612388/147‑1 (MQ=255)
                                                                          aGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                             >  1:377718/1‑151 (MQ=255)
                                                                             ggCCGGCGTGGGGGGAA‑TGGGGGGGTGCAG‑GTGGGTAATAATGGCTAAAAGGACCCCGCGGTGATGGGGGGCCTAAAAT‑ATAAT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                           <  2:322697/150‑1 (MQ=255)
                                                                             ggCCAGCGTCTGGGAAT‑TGGGGGCGTGCAG‑CTGGGTAAAAATGGCTACAATGACCCCGCGGTGAGGGTGGGCAAAAAAT‑ATAAT‑TTCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                           <  2:926258/150‑1 (MQ=255)
                                                                              gCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAACtattat                                                                           >  1:820936/1‑149 (MQ=255)
                                                                                 ggcgtggggggAA‑TGGGGTGGCGCAG‑CTGGTTGAAGAGGGCTACAATGACCCCGTGGTGAGGGTGGGGGTGAATT‑ATAAT‑TTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                      <  1:1485042/144‑1 (MQ=255)
                                                                                   cGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                    >  2:1792745/1‑151 (MQ=255)
                                                                                    gTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                   >  1:1139483/1‑151 (MQ=255)
                                                                                     tCTGGGGAA‑TGTGGG‑GGGCGGTGTGGGTGATAATGGCAACAATGACCCCGCGGTGTTGGTGGGCCTAAAAT‑ATAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                    <  2:2007430/149‑1 (MQ=255)
                                                                                      tcggggAATTGGGGGCGTGCCG‑CTGGGTGATAATGCCTACAATACCCCCGTGTTGAGGGTGGCCTTAAAAT‑ATAAG‑TTCTCCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                    <  2:1022323/147‑1 (MQ=255)
                                                                                      cTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                                  >  2:1802645/1‑150 (MQ=255)
                                                                                        ggggAT‑TGGGGGGGTGCAG‑GTGGGTAATAATGCCTACAATGACCCCGCGTTGAGGGTGGGCATAAAAT‑ATAAT‑TTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                      <  2:880086/144‑1 (MQ=255)
                                                                                        ggggAA‑GGGGGGCGTGCAG‑GGGGGTGATAATGGCTACAATGACCCCGCAGTGAGGGTGGGCCTAAAAT‑ATAAG‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGa                                                                 <  2:828605/149‑1 (MQ=255)
                                                                                         gggaattggggggGTGCAG‑CTGGGTGATAAGGGCTAAAATGACCCCGCAGTGATGGGGGGCCTAAAAT‑AAAAT‑TTCTTACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                    <  2:987091/137‑1 (MQ=255)
                                                                                               tgtgGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                         >  1:659678/1‑151 (MQ=255)
                                                                                               tgggggCGTGCAG‑CTGGGTGATAATGGCTACAATGACCCCGGAGTGTGGGGGGGCTAAAAAT‑ATAGT‑TTCTGGTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                          <  2:1360833/147‑1 (MQ=255)
                                                                                                 tgGGCGTGCAG‑CGGGGTAATAATGGCAACAATGACCCCGCAGTGTGGGTGGGCCTAAAAT‑ATAAT‑TTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCa                                                        <  2:1119860/150‑1 (MQ=255)
                                                                                                      gTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                   >  1:1940087/1‑150 (MQ=255)
                                                                                                       tGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtc                                                 >  1:536594/1‑151 (MQ=255)
                                                                                                       tGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGGGtc                                                 >  1:691446/1‑151 (MQ=255)
                                                                                                             cTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca                                             >  1:1473730/1‑150 (MQ=255)
                                                                                                              tGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                           >  2:1950946/1‑151 (MQ=255)
                                                                                                               gggTGATAATGGCTACAAGCACACGGCAGTGAGGGGGGGCCTAAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                           <  2:706415/150‑1 (MQ=255)
                                                                                                               gggTGAAAATGGCTAAAATGCCCCCGCGGTCATGGTGGGCCTAAAAT‑ATAGG‑TTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAtt                                          <  2:1167435/151‑1 (MQ=255)
                                                                                                               gggTGAAAAGGCCAAAAATGCCCCCGAAGTGATGGTGGGCCTAAAAT‑ATAGT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa                                            <  1:895976/149‑1 (MQ=255)
                                                                                                                    aaaaggccaaaaaTGACCCCGCGGTGGTGGGGGGCCAAAAAT‑AAAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                         <  1:192443/136‑1 (MQ=255)
                                                                                                                    aTAATGCCTACAATGCCCCCGGGGGGGGGGTGGGCCAAAAAT‑ATAGT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg                                      <  1:866264/150‑1 (MQ=255)
                                                                                                                        tGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATc                                 >  2:568729/1‑151 (MQ=255)
                                                                                                                          gCTACAATGACCCCGAAGAGAGGGGGGGCTTAAAAT‑ATAAT‑TTCTTATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                                   <  2:1000249/147‑1 (MQ=255)
                                                                                                                               aaTGACACCGCGGTGAGGGGGGGCTTAAAAA‑ATAAT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                        <  2:1086199/137‑1 (MQ=255)
                                                                                                                               aaTGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                           >  2:1026747/1‑150 (MQ=255)
                                                                                                                                aggACACCGCAGATATGGTGTGCC‑GAAAT‑AAAGGTTTCTGTTACCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGTCGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                            <  1:913501/146‑1 (MQ=255)
                                                                                                                                   accccGCAGTGATGGGGGCACAAAAAT‑ATATT‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt                        <  1:634603/146‑1 (MQ=255)
                                                                                                                                    cacCGCGGTGATGGTGGCCCTAAAAA‑ATAAT‑TTTTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt                        <  1:1635472/148‑1 (MQ=255)
                                                                                                                                      ccGGAGTGAGGGTGGGCTTAAAAT‑ATAAT‑TTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                     <  1:44795/149‑1 (MQ=255)
                                                                                                                                         cAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc                >  2:474318/1‑151 (MQ=255)
                                                                                                                                          agcgATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCt               >  1:1489672/4‑151 (MQ=255)
                                                                                                                                           gTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt              >  1:701111/1‑151 (MQ=255)
                                                                                                                                               tggtggGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc          >  1:1308379/1‑151 (MQ=255)
                                                                                                                                                ggtggGCGAGAAAT‑ATAAG‑TTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc          <  2:352853/150‑1 (MQ=255)
                                                                                                                                                ggtggGCCTAAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg           <  2:1374373/149‑1 (MQ=255)
                                                                                                                                                 gtggGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc        >  2:181020/1‑151 (MQ=255)
                                                                                                                                                  tggGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc        >  1:935780/1‑150 (MQ=255)
                                                                                                                                                  tggGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc        >  2:1601344/1‑150 (MQ=255)
                                                                                                                                                       cTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCt  >  1:707962/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                      
GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAA‑CCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAG‑CTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAAT‑ATAAG‑TTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCT  >  NC_000913/1222083‑1222381

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