Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 845,227 A→C 19.3% V85G (GTG→GGG)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913845,2270AC19.3% 106.3 / 2.9 66V85G (GTG→GGG) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base A (19/29);  new base C (10/5);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 8.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGATACAGCGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGA  >  NC_000913/844989‑845464
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
tcgaaccAGCGGCAGTGCGCACAGGCGAAACAGCCAGTAAAAACCAAACACTAATCCCGCGAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCa                                                                                                                                                                                                                                                   <  2:298441/233‑1 (MQ=255)
       aGCGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTTCACCCCGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCCC‑TAtttt                                                                                                                                                                                                                                           >  1:1308864/1‑239 (MQ=255)
         cGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTTCCTTCTGTGGCCTTTTGCAGGCGGGTTTGCTCTTC‑GACCCGCGCCGGC‑GGCCC‑TAtttttt                                                                                                                                                                                                                                         >  1:1069814/1‑238 (MQ=255)
           gCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGCGCCGGGGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCCTCCTCCGGGGACCCTTTGCGGCCCGGTTTTTCCTTC‑GCCCCGAGTCGGCGGGCCC‑TTTTTTTgt                                                                                                                                                                                                                                        >  2:1337251/1‑238 (MQ=255)
              gTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCCACGTCGGC‑GGGCC‑TATTTTTGGTAcc                                                                                                                                                                                                                                    >  1:633899/1‑239 (MQ=255)
                gcgcACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTAc                                                                                                                                                                                                                                     <  1:708891/236‑1 (MQ=255)
                gcgcAAAGGGAAACCAGCCAGTAAAAACCAAACACTAATACCGGTAACATGGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGGGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTAc                                                                                                                                                                                                                                     <  2:745123/236‑1 (MQ=255)
                  gaacaGGGGAATCAGCCAGTAAAAACCAAACCATAATGCCGCTAACATTGAAAAATGGGTCAGCGGATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGt                                                                                                                                                                                                                                 <  2:191379/236‑1 (MQ=255)
                          aaaTCAGCCAGTAAAAACCAAACCCTAATACCGCTAACATTGAAAAATGGGGCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCCC‑TATTTTTGGTACCGGTTCAGCAgg                                                                                                                                                                                                                         <  2:1303240/237‑1 (MQ=255)
                                  gagaaaaagCCAAACAGTAATAACACTCAACTGGAAGAAGGGGTCCGTGCATTGGTAAAG‑CTGTGTGGGTTAAACGGATTAGGGGTGGAGCCCGTGATATTGCGATAAAGTTGCCCGAAGCGGGCTGAAAGGGCTTCACCAGAATTGCGGCTGACTTCGTTGACCTTTTGCACCACGGTTTACGCTTC‑GACCAGCGTCGGC‑GGACC‑TATTTTTGGTAACGGTTCCGCAGGGGGCGt                                                                                                                                                                                                                   <  2:1083236/232‑1 (MQ=255)
                                       aaaaCCAAACCCTAACTCCGCAAACATGGAAAATGGGGTACGCGCATTGCTGAAG‑GTTTGTGGTATAAACGGATTATGCGGGGAGCCGGGTATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCaa                                                                                                                                                                                                          <  2:869740/239‑1 (MQ=255)
                                       aaaaCCAAACACTAATTCCGCTAACATTGAAAAATGGGTTAGCGCATTGCTGAAG‑GTTTGTGGATTAAACCGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCaa                                                                                                                                                                                                          <  2:422755/239‑1 (MQ=255)
                                            cAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GGTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTg                                                                                                                                                                                                     <  2:1192135/239‑1 (MQ=255)
                                               acacTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCGTTTGCAGCACGGGTTTGTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTTGTACCGGTTTCGCCAGGGGCGGAgg‑gggcaagaggacg                                                                                                                                                                                                  >  2:139031/1‑228 (MQ=255)
                                                     aTACCGCTAACATGGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAACCGGCTTATGCGGGGAGCCGGTGATAGTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTgg                                                                                                                                                                                            <  1:389702/239‑1 (MQ=255)
                                                           cTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGGCCTTTTGGAGGACGGTTTGCTCTTC‑GGCCCCCGGCGGC‑GGCAC‑TATTTTTGGTACCGGTTCACCAGGGGGCGTAGC‑GGCAAAGGGGCGCaaacgggcggtc                                                                                                                                                                                      >  2:18377/1‑229 (MQ=255)
                                                           cTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCCCTTC‑GACCAGCGTCGGC‑GGCCC‑TATTTTTGGTACCGGTTCCGCAGGGGGGCTAGC‑GGCCACGGTTGGCCACTGGGCGGTc                                                                                                                                                                                      >  1:362541/1‑239 (MQ=255)
                                                           cTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGGTTGCTCTTC‑GACCAGCGTCGGC‑GGCCC‑TATTTTTGGTACCGGTTCCGCAGGGGGAGTAGC‑GGCAACGGTGCGCAACTggtcggtc                                                                                                                                                                                      >  1:34381/1‑239 (MQ=255)
                                                            aaactgtGAAAAATGGGTGGGGCGATTGCTGAAG‑GTTTGGGGATTAAACGGGTTAGGGGGGGAGCCGGGGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGGGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCa                                                                                                                                                                                     <  2:1080112/233‑1 (MQ=255)
                                                             aaaagtGAAAAAGGGGTCAGGGCATTTCTGAAG‑GTTGGTGGATTAAACGGCTTGTGGGGGGAGCCGGGGATATGTCGATAAAGGTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGAACTTTTGCAGCACGGGTTGCTCTTC‑GACCAGCGGCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCaa                                                                                                                                                                                    <  2:1352694/234‑1 (MQ=255)
                                                             aaCATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCCTTTGCAGCACGGTTTGCTCTTCGGACCA‑CGCTGGC‑GGCCCTTATTTTTGACGCGGGTT‑AGGCCGGGGCGTAGC‑GGGAAAGGTGCGACCCTGGTCgcggaa                                                                                                                                                                                    >  2:110333/1‑234 (MQ=255)
                                                             aaCATGTAAAAAGTGGGCAGCGCAGTGGTGTAG‑GGTTGTGGATTAAACGGGTTTGGGGGGGGGGCGGTGAAATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCaa                                                                                                                                                                                    <  2:990284/239‑1 (MQ=255)
                                                                 ttGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCCCGGTTTGCTCTTC‑GACCAGGGTCTGC‑GGCCC‑TATTTTTGGTACCGGTTCAGCAGGGGGGGCAGC‑GGCCACGGGGGGCAACCGGTGGGTCaaacct                                                                                                                                                                                >  2:420347/1‑235 (MQ=255)
                                                                 ttGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAAC‑GGCAACCGTGGGCAACTGGTCCATCAActct                                                                                                                                                                                >  1:526452/1‑239 (MQ=255)
                                                                      aattggggtcgcgcATTGCTGAAG‑TTTTGTGGATTAAACGGCTTAGGGGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAc‑g                                                                                                                                                                          <  2:673671/229‑1 (MQ=255)
                                                                      aaaTGGGTCGGGGCTTTGTTGAAG‑GTTTGGGGAGTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGCCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTGTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAc‑g                                                                                                                                                                          <  2:564139/239‑1 (MQ=255)
                                                                          gggTCAGGCCATTGGTGGAGTGTTTG‑GGATTAAACTGCTTTGGGGGGGGGCGGGTGATATTGCGAAAAAGGTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAg                                                                                                                                                                      <  2:891517/239‑1 (MQ=255)
                                                                           ggTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAg                                                                                                                                                                      <  1:1337251/238‑1 (MQ=255)
                                                                                           aaG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGTCTTCGGTGACCTTTTTAAGCCCGGTTTGCTCTTC‑GACCGCCGCCGGG‑GGCGC‑TTTTTTTTGTTCCGGGTCATCAGGGGGCGTGGCGGGCAAC‑GTGCACACCTGGGCGGTCAACTCTTTAC‑GCGAGGGGGCAATATCCaagc                                                                                                                                                      >  2:565952/1‑235 (MQ=255)
                                                                                           aaG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTAACGGTTCAGCAGGGGGCGTAGG‑GGCAACGGTGCGCAACTGGTCGATCAACCCTTTAC‑GCGAGGGGGCATTTTCCAGcac                                                                                                                                                     >  1:1360511/1‑239 (MQ=255)
                                                                                                 ttGTGGATTAAACGGGTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGGGGGCAACTGGTCGATCAACTCTTTAC‑CCCAGGGGTCATTATCCAGCAAATcc                                                                                                                                                 >  1:801711/1‑238 (MQ=255)
                                                                                                  ttttgATTAAACCGCTTAGGGGGGGAGCCGGTTAATTTGGGATAAATGTGCCCGAAGCGGGCGGAAAGGGCTTCCCCATAATGGCGGCTGCCTTCGGTGACCTTTGTCAGCCCGGTTTGCTCTTC‑GACCAGCGTCGGC‑TGCCC‑TATTTGTGGTCCCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑c                                                                                                                                              <  2:1376391/235‑1 (MQ=255)
                                                                                                    tgGATTAAACGGTGTAGGGGGGGAGCCGGTGATATTCCGATAAAGTTGCCCGAAGCGGGGGGAAAGGGCTTCACCATAATGGCGGCTGACTCCGGTGACTTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCa                                                                                                                                            <  2:1069814/239‑1 (MQ=255)
                                                                                                         tAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCa                                                                                                                                                                                                                                                   <  2:917256/135‑1 (MQ=255)
                                                                                                         tAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCa                                                                                                                                                                                                                                                   >  1:917256/1‑135 (MQ=255)
                                                                                                         tAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGGAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCCGCACATCCG‑CCAGGgc                                                                                                                                        >  1:260634/1‑238 (MQ=255)
                                                                                                          aaaCGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCa                                                                                                                                                         <  1:28198/221‑1 (MQ=255)
                                                                                                                          gAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGGCTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GGCCAGCGGCCGC‑GGCAC‑TATTTTTGGGACCGGGTCAGCAGGGGGCGTAGC‑GGCGACGGTGGGCCACCGGTCGGTCCACCCTTTAC‑CCGAGGGGGCATTATCCAACACATCCGCCCCGTGCCCCACAGCCCCCttttt                                                                                                                       >  2:1206264/1‑239 (MQ=255)
                                                                                                                                     atTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGGTCAACTCTTTAC‑GCGAGGTGCCATTATCCAGCACATCCG‑CCAGTGCGCCTTAGGCCGCTTTTTTTTGTTCGata                                                                                                            >  1:655023/1‑238 (MQ=255)
                                                                                                                                         cGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCaa                                                                                                                                                                                    <  1:703717/164‑1 (MQ=255)
                                                                                                                                           aTAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGGTCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGt                                                                                                      >  1:859858/1‑238 (MQ=255)
                                                                                                                                                    cccGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGCTCAGCAAACCCTTAC‑CCAGGGGGGCCTTATCCAGCCCATCCG‑CCCAGGCGCCCTAAGCCCCTTTTTTTTTTTTGATATCCGGTTCCCGGGGCg                                                                                            >  2:486545/1‑239 (MQ=255)
                                                                                                                                                       gaggcgggcggAAAGGGCTTCACCATATTGGCGGGTGACTTGGGTGCCCTTTGGCAGACGGGTTTGCTCTTC‑GACCGGCGTGGGC‑GGCAC‑TATTTTGGGTACCGGTTAAGCAGGGGGCGTAGG‑GGCAACGGTGCGCAACTGGTCGATCACCTCTTTAC‑GGGGGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGtt                                                                                         <  2:1280161/236‑1 (MQ=255)
                                                                                                                                                          ggggggggAGAGGTTTTCCCCATAATGGGGGGTGGTTTGGGTGACCTTTTGCAGCAGGGTTGCCTCTGC‑GACCAGCGTCGGG‑GGCTA‑TTTTTGTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACGGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGtt                                                                                         <  2:260634/234‑1 (MQ=255)
                                                                                                                                                                tgAAAGGGTTTCCCCATAAGGGGGGGGGACCTGGGGGACCTTTTGCAGCGGGGTTTGCTCTCT‑GACCCGCGTGGGG‑GGCAC‑TATTTTTGGTACCGGTTCAGCGGGGGGCGTAGG‑GGCAACGGTGCGCCACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGtt                                                                                 <  2:1015845/237‑1 (MQ=255)
                                                                                                                                                                ggAAGGGGTTTCCCCTTATGGGGGGCTGTCTTGGGTGACCTTTTGCCGCCCGTTTTGCTTTTG‑GACCAGCGTCGGC‑GGCAC‑TTTTTTTGGTCCCGTTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCACCTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTg                                                                                <  2:34381/239‑1 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCCTAGGCCGCTTTTTTTTGTTCGGTATCCGGTTCCGGGGCCgggt                                                                                         >  1:776050/1‑220 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAg                                                                                      >  1:274492/1‑225 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCcgt                                                                         <  1:893906/238‑1 (MQ=255)
                                                                                                                                                                        cTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGGGTCATTATCCAGCACATCCG‑CCAGAGCGCCATAGGCCGCTTTTTTTTGTTCGGTATCCGGTTCCGGGGCCGGTTCAGTCGTTGGGTCcgtt                                                                        >  1:1127342/1‑239 (MQ=255)
                                                                                                                                                                          tcccaaaaaTGGGGGCTGACTTCGGTGCCCTTTTGCAGCCCGGTTTGCTCTTC‑GACCAGGGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTCGAGTCcgttgt                                                                      <  1:1052195/232‑1 (MQ=255)
                                                                                                                                                                                tAATGGCGGCTGACTTCGGTGACCTTTTGGAGGCCGGTTTGCCCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCcgttgtcgttgt                                                                <  1:18377/239‑1 (MQ=255)
                                                                                                                                                                                tAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCCACGGTGCGCCACTGGTCGGTCAACTCTTTAC‑GCGAGGGGGCATTTTCCAGCACCTCCG‑CCCGGGCGCCCTAGGCCGCTTTTTTTTTTTCGCTATCCGGTTCCGGGGCCGGTTCGGCCGTTGGGTCcgttttggttgt                                                                >  2:893906/1‑239 (MQ=255)
                                                                                                                                                                                tAATGGCGGCCGGCTTCGGTGAGCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGGCCGC‑GGCCC‑TATTTTTGGTACCGGTTCAGCTGGGGGCGTAGC‑GGCAACGGCGCGCAAATGGTCGATACACTTCTTAC‑GCGAGGGGTCATGATCCAGCACACACG‑CCCATGCGCCCTAGGGCGCGTTTTTTTGTTCGCAATCCGGTTCCGGGGGCGGGTCAGTGGTTGGGTGcgttgtcgttgg                                                                >  2:1303833/1‑238 (MQ=255)
                                                                                                                                                                                        gcTGACTTCGGTGACCTTTTGCAGCACGGTTTTCTCTTC‑GACCAGCGTCGGC‑GGCCC‑TATTTTTTGTACCGGTTCAGCCGGGGGCGTAGC‑CGCCACGGTGCGCACCTGGTCGTTCCACTCTTTACTGCG‑GGGGTCCTTATCCAGCCCCTCCG‑CCCGTCCGCCCTGGGCCCCTTTTTTTTTTTCCATATTCGGTTCCGGGGCCGGTTTCGGCGGTTAGTCCGTTTTTGTTGTgggttgca                                                        >  2:653988/1‑233 (MQ=255)
                                                                                                                                                                                        gcTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCAGTTCAGCAGGGGGCGGAGC‑GGCAACGGGGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGGGTCATTATCCAGCACATCCG‑CCAGGGCGCCCTAGGCCGCCTTTTTTTGTTCGTATTCCGGTTCCTGGGCCGGTTACGTGTTTGAGTCCCTTGTGTTggtggtgggta                                                        >  2:553922/1‑234 (MQ=255)
                                                                                                                                                                                           tgtttCGGTGACTTTTTGCTGCACGTTTTGCTCGTC‑GGCCAGGGTCGGC‑GGCCC‑TTTTTTTGGTACCGGTTCAGCGGGGGGCGTGGC‑GGCAACGGTGCGCAACTGGTCGATCAACTTTTTAC‑GCGAGGTGTCATTATCCAGCCCATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGc                                                     <  2:1308864/236‑1 (MQ=255)
                                                                                                                                                                                             cTTCGGTTACCTTTTGCAGCACGGTTTGCTCTTC‑GGCCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCg                                                   <  1:486545/239‑1 (MQ=255)
                                                                                                                                                                                                       ctttgtcagccCGGTTTGCTTTTC‑GTCCAGGGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGAt                                          <  1:17455/232‑1 (MQ=255)
                                                                                                                                                                                                           tGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACcgc                                     <  1:420347/239‑1 (MQ=255)
                                                                                                                                                                                                                cccttttttCTCTTC‑GCCCGGCGGCGGC‑GGCAC‑TATTGTGGGTCCCGTTTCGGCAGGCGGTGGGGC‑GGCAACGGTGCGAAAAGGGCCGACCAACCCTGTCC‑GCGAGGGTCCATTTCCCAGCCCACCCG‑CCCGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTGCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTggg                                <  2:691641/234‑1 (MQ=255)
                                                                                                                                                                                                                 aCGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAA‑TATTGTTGGGACCGGGTCAGCAGGGGGCGTAGC‑GGCAACGGGGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGGGTCATTGTCCAGCACATCCG‑CCAGTGAGCCACAGGCCGCATTTTTTTGTTCGATATCCGGTTCCGGGGCCGGGTCAGTCGTTATGTCCGTTGTGGTTGTGGGTGCAACCCCGGGTTGGGATACGGGGTgggg                               >  2:982364/1‑238 (MQ=255)
                                                                                                                                                                                                                     tttGCTTTCC‑GGCCCGCGTGGGG‑GGCCG‑TATTTTTGGTCCCGTTTCGGGAGGGGGCGTGGC‑GGAACCGGTGCGCACCTGGTGGATCAACTCTTTCC‑GCGGGGTGTTATTACCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGt                           <  2:1360511/239‑1 (MQ=255)
                                                                                                                                                                                                                             tC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTTTCTTTGTGGTTTTACAGCCGGGTTTGGTAACCCCTGGGGCAGGGGCGccccg                   >  2:800123/1‑237 (MQ=255)
                                                                                                                                                                                                                             tC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTTCCGTTGTGGTTTTCAAGCCGGGTATGGGTACCGCGTGGGCAGGGGGGccccc                   >  2:748295/1‑237 (MQ=255)
                                                                                                                                                                                                                                gACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGAGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGCGTCCGTTGTCGTTGTGGTTTTAACCCCGGGTGTGGATACCGCGCGGGCTGGTGCGCCcagca                 >  1:869151/1‑239 (MQ=255)
                                                                                                                                                                                                                                      cGTGGGC‑GCCAT‑TTTTTTTGGTCCCGTTTACGCGGGGGGCGTAGC‑GGCAACGGTGCGCACCTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCaga           <  2:859858/239‑1 (MQ=255)
                                                                                                                                                                                                                                       gTCGGG‑GGCAC‑TATTTTTGGTCCCGGTTCAGCGGGGGGAGTAGG‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GGGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCaga           <  1:553922/238‑1 (MQ=255)
                                                                                                                                                                                                                                         cggc‑ggCTT‑TTTTTTTGGTCCCTGTTCGGCAGGGGGCGTAGG‑GGCAACGGTGCGAAACTGGTCAATCACCTCTTTAC‑GCGAGGTGTCATTATCCAGCACTGCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCagaag         <  2:1127342/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 aC‑TTTTTTGGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACCCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGa  <  1:863415/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 aC‑TGTTTTTGGTGCGGTTTCAGCGGGGGGGGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTAGCCAGCCCATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGa  <  1:653988/238‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
GCGATACAGCGGCAGTGCGCACAGGCGAATCAGCCAGTAAAAACCAAACACTAATACCGCTAACATTGAAAAATGGGTCAGCGCATTGCTGAAG‑GTTTGTGGATTAAACGGCTTATGCGGGGAGCCGGTGATATTGCGATAAAGTTGCCCGAAGCGGGCGGAAAGGGCTTCACCATAATGGCGGCTGACTTCGGTGACCTTTTGCAGCACGGTTTGCTCTTC‑GACCAGCGTCGGC‑GGCAC‑TATTTTTGGTACCGGTTCAGCAGGGGGCGTAGC‑GGCAACGGTGCGCAACTGGTCGATCAACTCTTTAC‑GCGAGGTGTCATTATCCAGCACATCCG‑CCAGTGCGCCATAGGCCGCTTTTTTTTGTTCGATATCCGGTTCCGGGGCCGGTTCAGTCGTTGAGTCCGTTGTCGTTGTGGTTGTAACGCCGGGTATGGATACCGCGTGGGCAGGTGCGCCCAGCAGGCAGAAGAGGATGA  >  NC_000913/844989‑845464

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