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

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)
                                            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)
                                                            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)
                                                                                                                            cAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                                  >  1:624729/1‑143 (MQ=255)
                                                                                                                                  aatCCACCAGGCCTTTATTGCTCACTCCATTGAAACCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCTGGAATGGGCGCGAATGAAAGAAAACAACGCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                              <  1:726531/236‑1 (MQ=255)
                                                                                                                                       aCCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCAGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGCTGGCAGCGAAAGATAATCAAGCCTGCATCAACTGTCATAAAGGGATTGCcaacccg                                                                                       >  1:854315/1‑237 (MQ=255)
                                                                                                                                            aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGCCTGTCATAAAGGAATTGCCCACCAGTTAcc                                                                                  >  1:534168/1‑239 (MQ=255)
                                                                                                                                                       tCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCCGATGAAGGTGCCAGCGAAAGATAATCAATCCTCCATCGCCTGTCAAAAAGGAATTGCCCCCCAGTTACCGGAAATGAgaa                                                                       >  2:900990/1‑237 (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)
                                                                                                                                                                       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)
                                                                                                                                                                                                     aCTTGCCGAGCTGGAAGGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                           <  1:14760/238‑1 (MQ=255)
                                                                                                                                                                                                      cTTGCCGGGCGTGAAGGGGCGCGAATGAAAGAAAACAACTTGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGa                         <  1:75687/239‑1 (MQ=255)
                                                                                                                                                                                                              gCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCTTGAAGCAGCACGTCAGATGAAGGTgg                                                                                                                                            >  1:858452/1‑116 (MQ=255)
                                                                                                                                                                                                                 gaaaTGGTCGCAAATCAGAAAAATCACCTGGGCAACCGGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTTGTGaca              <  1:344165/237‑1 (MQ=255)
                                                                                                                                                                                                                    aTGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAACAGGTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGc           >  2:640545/1‑239 (MQ=255)
                                                                                                                                                                                                                    aTGGGCCAGAATGAAAGAAAAAAAGCCGGCACCCTGCCTGCCCTGGCATAACTACGATGCGATGGATCAGGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg            <  1:555004/238‑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 < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: