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

GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTA  >  NC_000913/1058252‑1058708
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                              
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGCAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGTTAAACCTGAAAAATTCCAACCCAAACGCGCGGAACTTGCCGAGCCTGATTGGGCGCGAATACAAGAAAACAACTCGGc                                                                                                                                                                                                                            >  1:538806/1‑239 (MQ=255)
gtgtATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCCAGGACATTGCCGAACGTGACTGGGCGCGAATGAAAGAAAACAACTCgcg                                                                                                                                                                                                                            >  1:871094/1‑237 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACACGGCACCCGGCCGc                                                                                                                                                                                                                  >  1:203115/1‑238 (MQ=255)
           atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAAAAAACCGGCAACCGGCCGc                                                                                                                                                                                                                  >  2:726531/1‑238 (MQ=255)
            tataAACGTTCGGTGCTTTTCAGAACCGCCTCGGGCGGGCGGGCGAAATGCCATGACTGTCATACCCCGCCGGATATTCCAGGCATGGGGAGGCCAAACCTGAAAGGAAGAAATAATATCTACCAGCCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCAAACCGCGGGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCACCCTGCCGCTc                                                                                                                                                                                                                <  2:4340/239‑1 (MQ=255)
             ataAACGGGGGGTTCATTTCCAGACCCCCTCCGGCGTGCGAGCTAAATGCAATGCCTGTAATATCCCGCCGGATATTCCAGGCATGGGGAAGCGCAAACGGAAAGGGGACAATGATAACTCCCAGCCCTTTATTGCTCACCCCATTGATACACCGGAAAAATTCGAAGCCAAACGCGGGGAACTTGCCGGGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCACCCTGCCGcccc                                                                                                                                                                                                               <  2:256085/239‑1 (MQ=255)
                 aCAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGcc                                                                                                                                                                                                           >  1:800379/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                        >  1:154785/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCCCGAATGAAAGAAAACAACTCGGCAACCCGCCGCCCCTGCCATa                                                                                                                                                                                                        >  2:122203/1‑239 (MQ=255)
                    gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACCCGCGGAACTTGCCGAGCGTGAATGGGCGCGACTGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCAtc                                                                                                                                                                                                        >  1:174589/1‑238 (MQ=255)
                       ggTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGATTGGGCCGGAATGAAAAAAAACAACTCGGCAACCTGACGCACCCGCCATaaaa                                                                                                                                                                                                     >  2:294726/1‑237 (MQ=255)
                           cATTTCCAAAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATACCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  <  1:122203/238‑1 (MQ=255)
                                          ccccggggtgggtgtgaagtcccATGCCTGTAATATCCCGCCGGATTTCCCGGGCATGGGAAACCCAAAACTGGAAGCGAGAAATGATATCTACCAGACCTTTATTGCGCACTCCATTGATACCCCTGAAAAATTCGAAGCAAAACGCGCGGAATTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATg                                                                                                                                                                                   <  2:392546/218‑1 (MQ=255)
                                            ccGGCGGGCGGTTTGAATGCCCTGCCTGACATATCCGCCCGGAAATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGAGGAACTTGCCGAGCGTGAAGGGGCGCGAATGAAAGAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCg                                                                                                                                                                                 <  1:294726/239‑1 (MQ=255)
                                                  tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                           >  1:392546/1‑239 (MQ=255)
                                                  tGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAACCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcag                                                                                                                                                                           >  1:662563/1‑239 (MQ=255)
                                                    cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGa                                                                                                                                                                                        >  1:765410/1‑224 (MQ=255)
                                                    cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAagcagc                                                                                                                                                                          >  1:824582/1‑238 (MQ=255)
                                                       ggttaattgccctgacgtgcaatCCCGCCGGGTATTTCCAGCCGTGGGAAAGCAAAAACGGAAAGGGTACAAGAATATTTACCAACCCTTTATTGCCCACCCCATTGATACCCCTGAAAAAGTCGAAGCAAAACGCGCGGAACTTGCCGGGCGTGAAGGGGCGCGAAGGAAAGAAAACAACCCGGCAGCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATc                                                                                                                                                                       <  2:869292/218‑1 (MQ=255)
                                                            cttccATTACTTTCATCTCCCGCCGATTTTCCCAGGAAGGGTGAAGCAAAAACTGTAAGCGAGCAATAATATCCACCAGACTTATATGGCTCACTCCATCGATCACCCTGAAAAATTAGAACCCAAACGGGCGCACCTTGCCGAGCGTGAAGGGGAGCAAAAGAACGAAAACAACTCGCCACCCCCCTGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAag                                                                                                                                                                 <  2:487258/236‑1 (MQ=255)
                                                                   gtttGTCCTACCCCGCGGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGAAATGATATCCACCAGACCTTTATGGCTCACTCAATTGATACACCTGAAAAATTCGAAGCAAAACGCGCGAAACTTGCCGGGCGTGAATGGGCGAGAATGAAAGAAAACAACTCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACg                                                                                                                                                          <  2:897240/236‑1 (MQ=255)
                                                                     cTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGATAGAAAACAATTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGTACCTCGGAAGCCGCACCCTGAAGAGGCACGTc                                                                                                                                                        >  1:31988/1‑239 (MQ=255)
                                                                         cATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGCTGGATCATGCGAACCAGCATCCTGAAGCAGCACGTCAGAt                                                                                                                                                    >  1:450797/1‑239 (MQ=255)
                                                                                  ccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAAATAAAGggggg                                                                                                                                           >  2:871632/1‑238 (MQ=255)
                                                                                                gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGCAAGAAAACAACTCGCCAACCTGCCGCTCCTGCCATAACTACGATGCGAGGGAGCATGCGAACCAGAACACTCAAGCACCACGCAGGAGGAAGGGGGGAGGGAAAGaaaaa                                                                                                                              >  2:465372/1‑238 (MQ=255)
                                                                                                gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATaa                                                                                                                               >  1:842668/1‑238 (MQ=255)
                                                                                                gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTCCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATACTTAAGCAGCACGTCAGGTGTAGGGAGCAGCGAAAGATaat                                                                                                                              >  1:842674/1‑239 (MQ=255)
                                                                                                             cAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCg                                                                                                                 >  1:456536/1‑239 (MQ=255)
                                                                                                                aCTGGAAAGCGAGAAATAATATCTCCCAGCCCTTTTTTTCTCATCCCATTGATACACCTGAAAAATTCAAACCAAAACGCGCGGACGTGGCGGGGCGTGAAGGGGCGCGAAAAAAAGAAAACAACTCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                               <  2:203115/239‑1 (MQ=255)
                                                                                                                           aGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGCTGCGATGGACCATCCGAAGCAGCATCCTGAACCAGCCCGTCAGATGCAAGGGGCAGCGAAAGGAAATCAAACCGACATCggcgttaataaggg                                                                                                    >  2:518299/1‑226 (MQ=255)
                                                                                                                             cAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  >  1:624729/1‑143 (MQ=255)
                                                                                                                             cAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGTGGCGATGGATCATGCGAAGCAACATCCTGCAGCAGCCCGCCGGCTGAAGGTGGCAGCGAAAGATAATCAAAccctcaccgcacgaccaaaaggaa                                                                                                  >  2:797006/1‑217 (MQ=255)
                                                                                                                                  tataagaaccagcttttATTGCTCACTCCATTGAACCCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGcc                                                                                             <  1:784899/225‑1 (MQ=255)
                                                                                                                                   aatCCACCAGGCCTTTATTGCTCACTCCATTGAAACCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCTGGAATGGGCGCGAATGAAAGAAAACAACGCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                              <  1:726531/236‑1 (MQ=255)
                                                                                                                                        aCCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCAGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGCTGGCAGCGAAAGATAATCAAGCCTGCATCAACTGTCATAAAGGGATTGCcaacccg                                                                                       >  1:854315/1‑237 (MQ=255)
                                                                                                                                            gccctttattccctcttCCTTTGATACACTTAAAAAATTCGAACACAAAAGCGGGGAACTCGCGGAGCTGGAATGGGCGAAAATGAAAAAAAACACCCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                                   <  2:456536/223‑1 (MQ=255)
                                                                                                                                             aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTAcc                                                                                  >  1:534168/1‑239 (MQ=255)
                                                                                                                                                       caaccgcaatgaaacaacCTGAAAAATGAGAGGCAAAACGCGTGAAACTTGCCGAGCGATAATGGGCAGAATGGAAAGAAACAATCTCGGAAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCGCCAGTTACCGGATATGAGt                                                                        <  1:137999/223‑1 (MQ=255)
                                                                                                                                                        tCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCCGATGAAGGTGCCAGCGAAAGATAATCAATCCTCCATCGCCTGTCAAAAAGGAATTGCCCCCCAGTTACCGGAAATGAgaa                                                                       >  2:900990/1‑237 (MQ=255)
                                                                                                                                                         cACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAACAGCATCCTGAAGCAGCACGTAAGATGAAGGGGGCCGCGAAAGATAATCAATCCTGCATCGACTGCCATAAAGGAATTGCCAACCAGTTAccccgtcccttata                                                                      >  2:784899/1‑227 (MQ=255)
                                                                                                                                                            tctagtagtACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAAAAAACAACGCGTCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCgg                                                                   <  1:115657/231‑1 (MQ=255)
                                                                                                                                                            tCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATCACTACGATGCGATGGATCATGCGAACCAGCATCCTGACGCAGAACGTCCGTGGAAGGTGGCAGGGAAAGATAATCAATCCTGCATCACCTGTCATAAAGGAATTTCCAACCAGTTACCGGATATCAGAAACgg                                                                   >  2:723751/1‑239 (MQ=255)
                                                                                                                                                                   aTACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGACAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATCTGAGTACGGGCTTCCGt                                                            >  1:345353/1‑239 (MQ=255)
                                                                                                                                                                    tACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAgg                                                                                                    >  1:138324/1‑198 (MQ=255)
                                                                                                                                                                    tACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCCCCAGTTACCGGATATGAGTAGCGGCTTCCGTa                                                           >  1:413012/1‑239 (MQ=255)
                                                                                                                                                                      ccccTAAAAAATTCAAAGCAAACGCGGCGAAACTTGCCGACGGTGAATGGGCCGGAATGAAAGAAACCAACCCGGAAACCTGCCGTCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAg                                                         <  1:358177/237‑1 (MQ=255)
                                                                                                                                                                      cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGCGGCAGCGAAAGATAATCAATCCTTCATCGACTGTCATAAAGGTATAGCCCACACGTTACAGGATATGAGTAGCGGCGTCCGTaa                                                          >  1:703837/1‑238 (MQ=255)
                                                                                                                                                                      cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCgg                                                                                >  1:357050/1‑216 (MQ=255)
                                                                                                                                                                      cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCCGTTACCGGATATGAGTAGCGGCTTCCGTaac                                                         >  1:255173/1‑238 (MQ=255)
                                                                                                                                                                        ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGACTGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAGGCAGCACGTCAGATGAAGGTGGCAGCGAACGAAAATCAATCCTGCATCGACTGTCATAAGGGTATTGCCCACCAGTTACCGAAAATGAGCAGAggacaccgaaagga                                                       >  2:323687/1‑227 (MQ=255)
                                                                                                                                                                        ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTACCGGATATGAGTAGCGGCCTCCGTaacaa                                                       >  1:916655/1‑236 (MQ=255)
                                                                                                                                                                        ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCCCCAGTTACCGGATATGATAAGCGGCTCCCGTAAGCa                                                       >  2:14760/1‑239 (MQ=255)
                                                                                                                                                                                 ttCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCCCCAGTTACCGGAAATGAGCAGCGGCTTCCGTAAAAAGTTCGCTGa                                              >  2:58496/1‑239 (MQ=255)
                                                                                                                                                                                   cGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGCTGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGCTGAAGGTGGCAGCGAACGATAATCAATCCTGCATCGACTGTCATAAAGGACTTGCCCACCAGTTACCGGATATGAGCAACGGCATCCGAAAGAAGTGCGAAGAGc                                            >  2:52195/1‑239 (MQ=255)
                                                                                                                                                                                           aaaccggcgGAAATTGCCGGGCGTAAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcgcgc                                     <  1:58496/233‑1 (MQ=255)
                                                                                                                                                                                                   gaacctgcccggggggaaaTGGGCACGAAGGAAAAAAAACACCTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCt                              <  1:465372/221‑1 (MQ=255)
                                                                                                                                                                                                      aCTTGCCGAGCTGGAAGGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                           <  1:14760/238‑1 (MQ=255)
                                                                                                                                                                                                       cTTGCCGGGCGTGAAGGGGCGCGAATGAAAGAAAACAACTTGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                         <  1:75687/239‑1 (MQ=255)
                                                                                                                                                                                                               gCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCTTGAAGCAGCACGTCAGATGAAGGTgg                                                                                                                                            >  1:858452/1‑116 (MQ=255)
                                                                                                                                                                                                                 tttatgggggggaaaTAAAAGAAACCCCCTCGGAACCCCGCCGCTCCTGACATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac               <  2:154785/226‑1 (MQ=255)
                                                                                                                                                                                                                  gaaaTGGTCGCAAATCAGAAAAATCACCTGGGCAACCGGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTTGTGaca              <  1:344165/237‑1 (MQ=255)
                                                                                                                                                                                                                   taaaggggggccatgaaaaaaaaaCAACTCGGCAACCGCCCGCTCCTGCCATAACCACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGac               <  2:916655/222‑1 (MQ=255)
                                                                                                                                                                                                                     aTGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAACAGGTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGc           >  2:640545/1‑239 (MQ=255)
                                                                                                                                                                                                                     aTGGGCCAGAATGAAAGAAAAAAAGCCGGCACCCTGCCTGCCCTGGCATAACTACGATGCGATGGATCAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg            <  1:555004/238‑1 (MQ=255)
                                                                                                                                                                                                                      agggggcagaaggaagaaaaaCAACCCGCCACCCTGCCGGTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTGCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg            <  1:723751/221‑1 (MQ=255)
                                                                                                                                                                                                                       gggaggaaaTGAAAGAAAACCACCTGGACACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTg         <  1:871632/232‑1 (MQ=255)
                                                                                                                                                                                                                           ggaaaaaaaaaaaaaaACACGCGGAAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa       <  2:413012/230‑1 (MQ=255)
                                                                                                                                                                                                                           gcgAATGAAAAAAAAAACACCGGGCAACTCGCCCACCCGGACATAACTACGATGAGATGGATGATGCAGAGCAACATACTGAAGCAGCAGGTCAGGTGAAGGTGGCAGCGAAAGATAATCAAACCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTa       <  2:407925/238‑1 (MQ=255)
                                                                                                                                                                                                                           ccaaaTGAAAGAAAACAACCCCTCACCCCCCCCCTCCTGCCAAAACTACGATGCAATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTAc      <  2:255173/235‑1 (MQ=255)
                                                                                                                                                                                                                               aTGAAAACAAACACCCCGGCAACCTGCAGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTa  <  1:473529/238‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                              
GTGTATGAAGAATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCTGTACTCTA  >  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: