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% 103.3 / 5.3 57S135P (TCG→CCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base T (34/8);  new base C (2/9);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 2.10e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTA  >  NC_000913/1058252‑1058708
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGCAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGTTAAACCTGAAAAATTCCAACCCAAACGC‑GCGGAACTTGCCGAGCCTGATTGGGCGCGAATACAAGAAAACA‑ACTCGGc                                                                                                                                                                                                                            >  1:538806/1‑239 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑CAGGACATTGCCGAACGTGACTGGGCGCGAATGAAAGAAAACA‑ACTCgcg                                                                                                                                                                                                                            >  1:871094/1‑237 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACACGGCACCCGGCCGc                                                                                                                                                                                                                  >  1:203115/1‑238 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAAAA‑AACCGGCAACCGGCCGc                                                                                                                                                                                                                  >  2:726531/1‑238 (MQ=255)
            tataAACGTTCGGTGC‑TTTTCAGAACCGCCTCGGGCGGGCGGGCGAAATGCCATGACTGTCATACCCCGCCGGATATTCCAGGCATGGGGAGGCCAAACCTGA‑AAGGAAGAAATAATATCTACCAGCCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCAAACCGC‑GGGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCACCCTGCCGCTc                                                                                                                                                                                                                <  2:4340/239‑1 (MQ=255)
             ataAACGGGGGGTTCATTTCCAGAC‑CCCCTCCGGCGTGCGAGCTAAATGCAATGCCTGTAATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACGGA‑AAGGGGACAATGATAACTCCCAGCCCTTTATTGCTCACCCCATTGATACACCGGAAAAATTCGAAGCCAAACGC‑GGGGAACTTGCCGGGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACCCGGCACCCTGCCGcccc                                                                                                                                                                                                               <  2:256085/239‑1 (MQ=255)
                 aCAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGcc                                                                                                                                                                                                           >  1:800379/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACCCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                        >  1:154785/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCCCGAATGAAAGAAAACA‑ACTCGGCAACCCGCCGCCCCTGCCATa                                                                                                                                                                                                        >  2:122203/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACCC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGACTGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCAtc                                                                                                                                                                                                        >  1:174589/1‑238 (MQ=255)
                       ggTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGATTGGGCCGGAATGAAAAAAAACA‑ACTCGGCAACCTGACGCACCCGCCATaaaa                                                                                                                                                                                                     >  2:294726/1‑237 (MQ=255)
                           cATTTCCAAAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATACCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  <  1:122203/238‑1 (MQ=255)
                                            ccGGCGGGCGGTTTGAATGCCCTGCCTGACATATCCGCCCGGAAATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GAGGAACTTGCCGAGCGTGAAGGGGCGCGAATGAAAGAAAACA‑ACCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCg                                                                                                                                                                                 <  1:294726/239‑1 (MQ=255)
                                                  tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                           >  1:392546/1‑239 (MQ=255)
                                                  tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAACCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAACCA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                           >  1:662563/1‑239 (MQ=255)
                                                    cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGa                                                                                                                                                                                        >  1:765410/1‑224 (MQ=255)
                                                    cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcagc                                                                                                                                                                          >  1:824582/1‑238 (MQ=255)
                                                            cttccATTACTTTCATCTCCCGCCGATTTTCCCAGGAAGGGTGAAGCAAAAACTGT‑AAGCGAGCAATAATATCCACCAGACTTATATGGCTCACTCCATCGATCACCCTGAAAAATTAGAACCCAAACGG‑GCGCACCTTGCCGAGCGTGAAGGGGAGCAAAAGAACGAAAACA‑ACTCGCCACCCCCCTGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAag                                                                                                                                                                 <  2:487258/236‑1 (MQ=255)
                                                                   gtttGTCCTACCCCGCGGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGAAATGATATCCACCAGACCTTTATGGCTCACTCAATTGATACACCTGAAAAATTCGAAGCAAAACGC‑GCGAAACTTGCCGGGCGTGAATGGGCGAGAATGAAAGAAAACA‑ACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACg                                                                                                                                                          <  2:897240/236‑1 (MQ=255)
                                                                     cTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGATAGAAAACA‑ATTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGTACCTCGGAAGCCGCACCCTGAAGAGGCACGTc                                                                                                                                                        >  1:31988/1‑239 (MQ=255)
                                                                         cATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGCTGGATCATGCGAACCAGCATCCTGAAGCAGCACGTCAGAt                                                                                                                                                    >  1:450797/1‑239 (MQ=255)
                                                                                  ccgGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAAATAAAGggggg                                                                                                                                           >  2:871632/1‑238 (MQ=255)
                                                                                               gCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGCAAGAAAACA‑ACTCGCCAACCTGCCGCTCCTGCCATAACTACGATGCGAGGGAGCATGCGAACCAGAACACTCAAGCACCACGCAGGAGGAAGGGGGGAGGGAAAGaaaaa                                                                                                                              >  2:465372/1‑238 (MQ=255)
                                                                                               gCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaa                                                                                                                               >  1:842668/1‑238 (MQ=255)
                                                                                               gCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTCCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATACTTAAGCAGCACGTCAGGTGTAGGGAGCAGCGAAAGATaat                                                                                                                              >  1:842674/1‑239 (MQ=255)
                                                                                                            cAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCg                                                                                                                 >  1:456536/1‑239 (MQ=255)
                                                                                                               aCTGGAAAGCGAGAAATAATATCTCCCAGCCCTTTTTTTCTCATCCCATTGATACACCTGAAAAATTCAAACCAAAACGC‑GCGGACGTGGCGGGGCGTGAAGGGGCGCGAAAAAAAGAAAACA‑ACTCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                               <  2:203115/239‑1 (MQ=255)
                                                                                                                            cAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  >  1:624729/1‑143 (MQ=255)
                                                                                                                                  aatCCACCAGGCCTTTATTGCTCACTCCATTGAAACCCCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCTGGAATGGGCGCGAATGAAAGAAAACA‑ACGCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                              <  1:726531/236‑1 (MQ=255)
                                                                                                                                       aCCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCAGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGCTGGCAGCGAAAGATAATCAAGCCTGCATCAACTGTCATAAAGGGATTGCcaacccg                                                                                       >  1:854315/1‑237 (MQ=255)
                                                                                                                                            aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTAcc                                                                                  >  1:534168/1‑239 (MQ=255)
                                                                                                                                                       tCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCCGATGAAGGTGCCAGCGAAAGATAATCAATCCTCCATCGCCTGTCAAAAAGGAATTGCCCCCCAGTTACCGGAAATGAgaa                                                                       >  2:900990/1‑237 (MQ=255)
                                                                                                                                                           tctagtagtACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAAAAAACA‑ACGCGTCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                   <  1:115657/231‑1 (MQ=255)
                                                                                                                                                           tCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATCACTACGATGCGATGGATCATGCGAACCAGCATCCTGACGCAGAACGTCCGTGGAAGGTGGCAGGGAAAGATAATCAATCCTGCATCACCTGTCATAAAGGAATTTCCAACCAGTTACCGGATATCAGAAACgg                                                                   >  2:723751/1‑239 (MQ=255)
                                                                                                                                                                  aTACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGACAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATCTGAGTACGGGCTTCCGt                                                            >  1:345353/1‑239 (MQ=255)
                                                                                                                                                                   tACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAgg                                                                                                    >  1:138324/1‑198 (MQ=255)
                                                                                                                                                                   tACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCCCCAGTTACCGGATATGAGTAGCGGCTTCCGTa                                                           >  1:413012/1‑239 (MQ=255)
                                                                                                                                                                     ccccTAAAAAATTCAAAG‑CAAACGCGGCGAAACTTGCCGACGGTGAATGGGCCGGAATGAAAGAAACCA‑ACCCGGAAACCTGCCGTCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAg                                                         <  1:358177/237‑1 (MQ=255)
                                                                                                                                                                     cacCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGCGGCAGCGAAAGATAATCAATCCTTCATCGACTGTCATAAAGGTATAGCCCACACGTTACAGGATATGAGTAGCGGCGTCCGTaa                                                          >  1:703837/1‑238 (MQ=255)
                                                                                                                                                                     cacCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCgg                                                                                >  1:357050/1‑216 (MQ=255)
                                                                                                                                                                     cacCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCCGTTACCGGATATGAGTAGCGGCTTCCGTaac                                                         >  1:255173/1‑238 (MQ=255)
                                                                                                                                                                       ccTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCCTCCGTaacaa                                                       >  1:916655/1‑236 (MQ=255)
                                                                                                                                                                       ccTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCCCCAGTTACCGGATATGATAAGCGGCTCCCGTAAGCa                                                       >  2:14760/1‑239 (MQ=255)
                                                                                                                                                                                ttCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCCCCAGTTACCGGAAATGAGCAGCGGCTTCCGTAAAAAGTTCGCTGa                                              >  2:58496/1‑239 (MQ=255)
                                                                                                                                                                                  cGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGCTGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGTGGCAGCGAACGATAATCAATCCTGCATCGACTGTCATAAAGGACTTGCCCACCAGTTACCGGATATGAGCAACGGCATCCGAAAGAAGTGCGAAGAGc                                            >  2:52195/1‑239 (MQ=255)
                                                                                                                                                                                          aaaccggcgGAAATTGCCGGGCGTAAATGGGCGCGAATGAAAGAAAACA‑ACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcgcgc                                     <  1:58496/233‑1 (MQ=255)
                                                                                                                                                                                                     aCTTGCCGAGCTGGAAGGGGCGCGAATGAAAGAAACCA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                           <  1:14760/238‑1 (MQ=255)
                                                                                                                                                                                                      cTTGCCGGGCGTGAAGGGGCGCGAATGAAAGAAAACA‑ACTTGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                         <  1:75687/239‑1 (MQ=255)
                                                                                                                                                                                                              gCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCTTGAAGCAGCACGTCAGATGAAGGTgg                                                                                                                                            >  1:858452/1‑116 (MQ=255)
                                                                                                                                                                                                                 gaaaTGGTCGCAAATCAGAAAAATCA‑CCTGGGCAACCGGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTTGTGaca              <  1:344165/237‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAACAGGTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGc           >  2:640545/1‑239 (MQ=255)
                                                                                                                                                                                                                    aTGGGCCAGAATGAAAGAAAAAA‑AGCCGGCACCCTGCCTGCCCTGGCATAACTACGATGCGATGGATCAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg            <  1:555004/238‑1 (MQ=255)
                                                                                                                                                                                                                      gggaggaaaTGAAAGAAAACC‑ACCTGGACACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTg         <  1:871632/232‑1 (MQ=255)
                                                                                                                                                                                                                         ggaaaaaaaaaaaaaaAC‑ACGCGGAAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa       <  2:413012/230‑1 (MQ=255)
                                                                                                                                                                                                                          gcgAATGAAAAAAAAAACACCGGGCAACTCGCCCACCCGGACATAACTACGATGAGATGGATGATGCAGAGCAACATACTGAAGCAGCAGGTCAGGTGAAGGTGGCAGCGAAAGATAATCAAACCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa       <  2:407925/238‑1 (MQ=255)
                                                                                                                                                                                                                          ccaaaTGAAAGAAAACA‑ACCCCTCACCCCCCCCCTCCTGCCAAAACTACGATGCAATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTAc      <  2:255173/235‑1 (MQ=255)
                                                                                                                                                                                                                              aTGAAAACAAACA‑CCCCGGCAACCTGCAGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTa  <  1:473529/238‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                              
GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAA‑CGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGG‑AAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGC‑GCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACA‑ACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTA  >  NC_000913/1058252‑1058708

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