Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,486 T→C 19.2% S135P (TCG→CCG)  torC → trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,058,4860TC19.2% 120.7 / 6.5 70S135P (TCG→CCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base T (38/14);  new base C (2/12);  total (41/29)
Fisher's exact test for biased strand distribution p-value = 1.03e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTA  >  NC_000913/1058252‑1058708
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                              
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGCAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGTTAAACCTGAAAAATTCCAACCCAAACGC‑GCGGAACTTGCCGAGCCTGATTGGGCGCGAAT‑ACAAGAAAACA‑ACTCGGc                                                                                                                                                                                                                            >  1:538806/1‑239 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑CAGGACATTGCCGAACGTGACTGGGCGCGAAT‑GAAAGAAAACA‑ACTCgcg                                                                                                                                                                                                                            >  1:871094/1‑237 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACACGGCACCCGGCCGc                                                                                                                                                                                                                  >  1:203115/1‑238 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAAAA‑AACCGGCAACCGGCCGc                                                                                                                                                                                                                  >  2:726531/1‑238 (MQ=255)
            tataAACGTTCGGTGC‑TTTTCAGAACCGCCTCGGGCGGGCGGGCGAAATGCCATGACTGTCATACCCCGCCGGATA‑TTCCAGGCATGGGGAGGCCAAACCTGA‑AAGGAAGAAATAATATCTACCAGCCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCAAACCGC‑GGGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCACCCTGCCGCTc                                                                                                                                                                                                                <  2:4340/239‑1 (MQ=255)
             ataAACGGGGGGTTCATTTCCAGAC‑CCCCTCCGGCGTGCGAGCTAAATGCAATGCCTGTAATATCCCGCCGGATA‑TTCCAGGCATGGGGAAGCGCAAACGGA‑AAGGGGACAATGATAACTCCCAGCCCTTTATTGCTCACCCCATTGATACACCGGAAAAATTCGAAGCCAAACGC‑GGGGAACTTGCCGGGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACCCGGCACCCTGCCGcccc                                                                                                                                                                                                               <  2:256085/239‑1 (MQ=255)
                 aCAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGcc                                                                                                                                                                                                           >  1:800379/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACCCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                        >  1:154785/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCCCGAAT‑GAAAGAAAACA‑ACTCGGCAACCCGCCGCCCCTGCCATa                                                                                                                                                                                                        >  2:122203/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACCC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGACT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCAtc                                                                                                                                                                                                        >  1:174589/1‑238 (MQ=255)
                       ggTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGATTGGGCCGGAAT‑GAAAAAAAACA‑ACTCGGCAACCTGACGCACCCGCCATaaaa                                                                                                                                                                                                     >  2:294726/1‑237 (MQ=255)
                           cATTTCCAAAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATACCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  <  1:122203/238‑1 (MQ=255)
                                          ccccggggtgggtgtgaagtcccATGCCTGTAATATCCCGCCGGATT‑TCCCGGGCATGGGAAACCCAAAACTGG‑AAGCGAGAAATGATATCTACCAGACCTTTATTGCGCACTCCATTGATACCCCTGAAAAATTCGAAGCAAAACGC‑GCGGAATTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATg                                                                                                                                                                                   <  2:392546/218‑1 (MQ=255)
                                            ccGGCGGGCGGTTTGAATGCCCTGCCTGACATATCCGCCCGGAAA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GAGGAACTTGCCGAGCGTGAAGGGGCGCGAAT‑GAAAGAAAACA‑ACCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCg                                                                                                                                                                                 <  1:294726/239‑1 (MQ=255)
                                                  tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                           >  1:392546/1‑239 (MQ=255)
                                                  tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAACCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAACCA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                           >  1:662563/1‑239 (MQ=255)
                                                    cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGa                                                                                                                                                                                        >  1:765410/1‑224 (MQ=255)
                                                    cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcagc                                                                                                                                                                          >  1:824582/1‑238 (MQ=255)
                                                       ggttaattgccctgacgtgcaatCCCGCCGGGTATTTCCAGCCGTGGGAAAGCAAAAACGGA‑AAGGGTACAAGAATATTTACCAACCCTTTATTGCCCACCCCATTGATACCCCTGAAAAAGTCGAAGCAAAACGC‑GCGGAACTTGCCGGGCGTGAAGGGGCGCGAAG‑GAAAGAAAACA‑ACCCGGCAGCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATc                                                                                                                                                                       <  2:869292/218‑1 (MQ=255)
                                                            cttccATTACTTTCATCTCCCGCCGATTT‑TCCCAGGAAGGGTGAAGCAAAAACTGT‑AAGCGAGCAATAATATCCACCAGACTTATATGGCTCACTCCATCGATCACCCTGAAAAATTAGAACCCAAACGG‑GCGCACCTTGCCGAGCGTGAAGGGGAGCAAAA‑GAACGAAAACA‑ACTCGCCACCCCCCTGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAag                                                                                                                                                                 <  2:487258/236‑1 (MQ=255)
                                                                   gtttGTCCTACCCCGCGGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGAAATGATATCCACCAGACCTTTATGGCTCACTCAATTGATACACCTGAAAAATTCGAAGCAAAACGC‑GCGAAACTTGCCGGGCGTGAATGGGCGAGAAT‑GAAAGAAAACA‑ACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACg                                                                                                                                                          <  2:897240/236‑1 (MQ=255)
                                                                     cTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GATAGAAAACA‑ATTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGTACCTCGGAAGCCGCACCCTGAAGAGGCACGTc                                                                                                                                                        >  1:31988/1‑239 (MQ=255)
                                                                         cATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGCTGGATCATGCGAACCAGCATCCTGAAGCAGCACGTCAGAt                                                                                                                                                    >  1:450797/1‑239 (MQ=255)
                                                                                  ccgGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAAATAAAGggggg                                                                                                                                           >  2:871632/1‑238 (MQ=255)
                                                                                                gCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GCAAGAAAACA‑ACTCGCCAACCTGCCGCTCCTGCCATAACTACGATGCGAGGGAGCATGCGAACCAGAACACTCAAGCACCACGCAGGAGGAAGGGGGGAGGGAAAGaaaaa                                                                                                                              >  2:465372/1‑238 (MQ=255)
                                                                                                gCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaa                                                                                                                               >  1:842668/1‑238 (MQ=255)
                                                                                                gCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTCCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATACTTAAGCAGCACGTCAGGTGTAGGGAGCAGCGAAAGATaat                                                                                                                              >  1:842674/1‑239 (MQ=255)
                                                                                                             cAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCg                                                                                                                 >  1:456536/1‑239 (MQ=255)
                                                                                                                aCTGGAAAGCGAGAAATAATATCTCCCAGCCCTTTTTTTCTCATCCCATTGATACACCTGAAAAATTCAAACCAAAACGC‑GCGGACGTGGCGGGGCGTGAAGGGGCGCGAAA‑AAAAGAAAACA‑ACTCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                               <  2:203115/239‑1 (MQ=255)
                                                                                                                           aGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAACCA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGCTGCGATGGACCATCCGAAGCAGCATCCTGAACCAGCCCGTCAGATGCAAGGGGCAGCGAAAGGAAATCAAACCGACATCggcgttaataaggg                                                                                                    >  2:518299/1‑226 (MQ=255)
                                                                                                                             cAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  >  1:624729/1‑143 (MQ=255)
                                                                                                                             cAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGTGGCGATGGATCATGCGAAGCAACATCCTGCAGCAGCCCGCCGGCTGAAGGTGGCAGCGAAAGATAATCAAAccctcaccgcacgaccaaaaggaa                                                                                                  >  2:797006/1‑217 (MQ=255)
                                                                                                                                  tataagaaccagcttttATTGCTCACTCCATTGAACCCCCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGcc                                                                                             <  1:784899/225‑1 (MQ=255)
                                                                                                                                   aatCCACCAGGCCTTTATTGCTCACTCCATTGAAACCCCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCTGGAATGGGCGCGAAT‑GAAAGAAAACA‑ACGCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                              <  1:726531/236‑1 (MQ=255)
                                                                                                                                        aCCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCAGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGCTGGCAGCGAAAGATAATCAAGCCTGCATCAACTGTCATAAAGGGATTGCcaacccg                                                                                       >  1:854315/1‑237 (MQ=255)
                                                                                                                                            gccctttattccctcttCCTTTGATACACTTAAAAAATTCGAACACAAAAGC‑GGGGAACTCGCGGAGCTGGAATGGGCGAAAAT‑GAAAAAAAACA‑CCCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                                   <  2:456536/223‑1 (MQ=255)
                                                                                                                                             aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTAcc                                                                                  >  1:534168/1‑239 (MQ=255)
                                                                                                                                                       caaccgcaatgaaacaacCTGAAAAATGAGAGGCAAAACGC‑GTGAAACTTGCCGAGCGATAATGGGC‑AGAATGGAAAGAAACAA‑TCTCGGAAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCGCCAGTTACCGGATATGAGt                                                                        <  1:137999/223‑1 (MQ=255)
                                                                                                                                                        tCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCCGATGAAGGTGCCAGCGAAAGATAATCAATCCTCCATCGCCTGTCAAAAAGGAATTGCCCCCCAGTTACCGGAAATGAgaa                                                                       >  2:900990/1‑237 (MQ=255)
                                                                                                                                                         cACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAACAGCATCCTGAAGCAGCACGTAAGATGAAGGGGGCCGCGAAAGATAATCAATCCTGCATCGACTGCCATAAAGGAATTGCCAACCAGTTAccccgtcccttata                                                                      >  2:784899/1‑227 (MQ=255)
                                                                                                                                                            tctagtagtACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAAAAAACA‑ACGCGTCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                   <  1:115657/231‑1 (MQ=255)
                                                                                                                                                            tCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATCACTACGATGCGATGGATCATGCGAACCAGCATCCTGACGCAGAACGTCCGTGGAAGGTGGCAGGGAAAGATAATCAATCCTGCATCACCTGTCATAAAGGAATTTCCAACCAGTTACCGGATATCAGAAACgg                                                                   >  2:723751/1‑239 (MQ=255)
                                                                                                                                                                   aTACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GACAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATCTGAGTACGGGCTTCCGt                                                            >  1:345353/1‑239 (MQ=255)
                                                                                                                                                                    tACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAgg                                                                                                    >  1:138324/1‑198 (MQ=255)
                                                                                                                                                                    tACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCCCCAGTTACCGGATATGAGTAGCGGCTTCCGTa                                                           >  1:413012/1‑239 (MQ=255)
                                                                                                                                                                      ccccTAAAAAATTCAAAG‑CAAACGCGGCGAAACTTGCCGACGGTGAATGGGCCGGAAT‑GAAAGAAACCA‑ACCCGGAAACCTGCCGTCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAg                                                         <  1:358177/237‑1 (MQ=255)
                                                                                                                                                                      cacCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGCGGCAGCGAAAGATAATCAATCCTTCATCGACTGTCATAAAGGTATAGCCCACACGTTACAGGATATGAGTAGCGGCGTCCGTaa                                                          >  1:703837/1‑238 (MQ=255)
                                                                                                                                                                      cacCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCgg                                                                                >  1:357050/1‑216 (MQ=255)
                                                                                                                                                                      cacCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCCGTTACCGGATATGAGTAGCGGCTTCCGTaac                                                         >  1:255173/1‑238 (MQ=255)
                                                                                                                                                                        ccTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGACTGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAGGCAGCACGTCAGATGAAGGTGGCAGCGAACGAAAATCAATCCTGCATCGACTGTCATAAGGGTATTGCCCACCAGTTACCGAAAATGAGCAGAggacaccgaaagga                                                       >  2:323687/1‑227 (MQ=255)
                                                                                                                                                                        ccTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCCTCCGTaacaa                                                       >  1:916655/1‑236 (MQ=255)
                                                                                                                                                                        ccTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCCCCAGTTACCGGATATGATAAGCGGCTCCCGTAAGCa                                                       >  2:14760/1‑239 (MQ=255)
                                                                                                                                                                                 ttCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCCCCAGTTACCGGAAATGAGCAGCGGCTTCCGTAAAAAGTTCGCTGa                                              >  2:58496/1‑239 (MQ=255)
                                                                                                                                                                                   cGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGCTGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGTGGCAGCGAACGATAATCAATCCTGCATCGACTGTCATAAAGGACTTGCCCACCAGTTACCGGATATGAGCAACGGCATCCGAAAGAAGTGCGAAGAGc                                            >  2:52195/1‑239 (MQ=255)
                                                                                                                                                                                           aaaccggcgGAAATTGCCGGGCGTAAATGGGCGCGAAT‑GAAAGAAAACA‑ACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcgcgc                                     <  1:58496/233‑1 (MQ=255)
                                                                                                                                                                                                   gaacctgcccggggggaaaTGGGCACGAAG‑GAAAAAAAACA‑CCTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCt                              <  1:465372/221‑1 (MQ=255)
                                                                                                                                                                                                      aCTTGCCGAGCTGGAAGGGGCGCGAAT‑GAAAGAAACCA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                           <  1:14760/238‑1 (MQ=255)
                                                                                                                                                                                                       cTTGCCGGGCGTGAAGGGGCGCGAAT‑GAAAGAAAACA‑ACTTGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                         <  1:75687/239‑1 (MQ=255)
                                                                                                                                                                                                               gCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCTTGAAGCAGCACGTCAGATGAAGGTgg                                                                                                                                            >  1:858452/1‑116 (MQ=255)
                                                                                                                                                                                                                 tttatgggggggaaaT‑AAAAGAAACCC‑CCTCGGAACCCCGCCGCTCCTGACATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac               <  2:154785/226‑1 (MQ=255)
                                                                                                                                                                                                                  gaaaTGGTCGCAAAT‑CAGAAAAATCA‑CCTGGGCAACCGGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTTGTGaca              <  1:344165/237‑1 (MQ=255)
                                                                                                                                                                                                                   taaaggggggccatgaaaaaaaaaCA‑ACTCGGCAACCGCCCGCTCCTGCCATAACCACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac               <  2:916655/222‑1 (MQ=255)
                                                                                                                                                                                                                     aTGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAACAGGTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGc           >  2:640545/1‑239 (MQ=255)
                                                                                                                                                                                                                     aTGGGCCAGAAT‑GAAAGAAAAAA‑AGCCGGCACCCTGCCTGCCCTGGCATAACTACGATGCGATGGATCAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg            <  1:555004/238‑1 (MQ=255)
                                                                                                                                                                                                                      agggggcagaaggaagaaaaaCA‑ACCCGCCACCCTGCCGGTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTGCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg            <  1:723751/221‑1 (MQ=255)
                                                                                                                                                                                                                       gggaggaaaT‑GAAAGAAAACC‑ACCTGGACACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTg         <  1:871632/232‑1 (MQ=255)
                                                                                                                                                                                                                           ggaaaaaaaaaaaaaaAC‑ACGCGGAAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa       <  2:413012/230‑1 (MQ=255)
                                                                                                                                                                                                                           gcgAAT‑GAAAAAAAAAACACCGGGCAACTCGCCCACCCGGACATAACTACGATGAGATGGATGATGCAGAGCAACATACTGAAGCAGCAGGTCAGGTGAAGGTGGCAGCGAAAGATAATCAAACCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa       <  2:407925/238‑1 (MQ=255)
                                                                                                                                                                                                                           ccaaaT‑GAAAGAAAACA‑ACCCCTCACCCCCCCCCTCCTGCCAAAACTACGATGCAATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTAc      <  2:255173/235‑1 (MQ=255)
                                                                                                                                                                                                                               aT‑GAAAACAAACA‑CCCCGGCAACCTGCAGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTa  <  1:473529/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                              
GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATA‑TTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAAT‑GAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTA  >  NC_000913/1058252‑1058708

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