Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,486 T→C 15.5% 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,4860TC15.5% 104.1 / 3.0 54S135P (TCG→CCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base T (30/11);  new base C (3/6);  total (33/21)
Fisher's exact test for biased strand distribution p-value = 4.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCT  >  NC_000913/1058263‑1058700
                                                                                                                                                                                                                                |                                                                                                                                                                                                                       
atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCTGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAAAAAAAAACTCgg                                                                                                                                                                                                                      >  2:657249/1‑227 (MQ=255)
atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCCCGGAACTTGCCGACCGTGAATGGGCGCGAATGAAAGAAAAAAACTCgg                                                                                                                                                                                                                      >  1:449064/1‑227 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGCTATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAATCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAAAAACTCGGAAAACTCCCg                                                                                                                                                                                                            >  2:515503/1‑228 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAATCGAGCGATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCCGGGGGGCACGGGCGGGAATGAAAGAAAAAAAATCGCCAACccgccg                                                                                                                                                                                                            >  2:641149/1‑228 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGGGAATGGCGGCGAATGAAAGAAAACAACTCGCCAAccgcccg                                                                                                                                                                                                            >  1:133233/1‑223 (MQ=255)
          tCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGCGCGTGAATGGGCGCGAATGAAAAAAAAAAACCCGGAAACCTGCCgg                                                                                                                                                                                                           >  2:282129/1‑227 (MQ=255)
           cGGTGCATTGCCAGACCGCCGCCGGCGTGCGAGCTGAAGGCCATGACTGTCATCTCCCGCAGGATATCCCAGGCATGGTGAAGCGAAAACTGGAAGCGAGCAATGATATCTACCAGTCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCAAATGAAAGAAAACAACTCGGCAACCTGCCGCt                                                                                                                                                                                                          <  1:190932/228‑1 (MQ=255)
                    tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                 >  1:497711/1‑228 (MQ=255)
                      ccgaccGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGCATACCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACCCCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCAt                                                                                                                                                                                                  <  1:488565/220‑1 (MQ=255)
                        gAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCGTAACTa                                                                                                                                                                                             >  1:208547/1‑228 (MQ=255)
                           cGCCTCCGGCGGGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                           <  1:87488/227‑1 (MQ=255)
                                        cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGCCTCGAATGAAAGAAAACAACCCGGCAACCCGCCGGCCCTGCCCTAACTACGAAGCAAAGGAACa                                                                                                                                                                              >  2:119057/1‑227 (MQ=255)
                                                    caggcgtGTCATAACCCGCCGGTTATTCCAGGCATGGGGAAGCAAAACGTGGAAGCGAGCAGTGATATCTGCCAGACCTTTATTGGCCACTCCAGTGATGCACCTGAAAAATTCGAAGCCAAGCGCCCGGACCTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCTCTCCTGCCATAACTACGATGCGATGGATCATGCGAAgcagca                                                                                                                                                                  <  1:515503/221‑1 (MQ=255)
                                                            tCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGTTGCGATGGATCATGCGAAGCAGCATCCTGAagc                                                                                                                                                         >  1:561719/1‑228 (MQ=255)
                                                                    cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGGGAAGCAGCATCCTGAAGCAGCACGTc                                                                                                                                                >  1:558378/1‑228 (MQ=255)
                                                                    cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAACCAGCATCCTGAAGCAGCCCGTc                                                                                                                                                >  1:205604/1‑228 (MQ=255)
                                                                                   gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGAAGGACCATTCGAAGCACCATCCTAAAGCAGCACGTCAGATaaaagaggaag                                                                                                                                 >  2:635779/1‑220 (MQ=255)
                                                                                       ggTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCGGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGGTCATGCGAAGCAGCATCCTGAAGCAGACCGTCAGATGAAGGTGGCAGCGa                                                                                                                              >  1:41716/1‑227 (MQ=255)
                                                                                        gTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGTGCGCGAATGCAAGAAAACAACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCCATGGATCATGCAAAGCAGCACCCTTCCCAACACCGCCAGATGAAGGGGGCGGCGaaa                                                                                                                            >  2:59041/1‑228 (MQ=255)
                                                                                             gcgcAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGCTa                                                                                                                        >  1:413620/1‑227 (MQ=255)
                                                                                              cgcAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGAACATGCGAAGCAGCATCCTGAACAACCACCCCAGAGAAAGGTGGCAGCGAAAGATaat                                                                                                                      >  2:337363/1‑228 (MQ=255)
                                                                                                   aCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGCCAGATGAAGGTGGCAGCGAAAGATaatcaatc                                                                                                                 >  1:24285/1‑228 (MQ=255)
                                                                                                   aCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGACCATGCGAAGCAGCAACCTGAAGCACCACGCCAACGGAAGGGGGACGCGAAAAAAaataaaac                                                                                                                 >  2:207114/1‑226 (MQ=255)
                                                                                                             atcAATGATATCTACCAGACCTTTATTGCACACTCCATTGACACACCTGAAAAATTCGAAGCCAACCGCGCGGAACTTGCCGAGCGTGAAGGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                       <  1:119057/226‑1 (MQ=255)
                                                                                                             aGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATCCGAAGCAGCATCCTGAAGCAGCACGGCAGATAAAGGTGGCACCGAAGAAAAACCAATCCTgaacagcc                                                                                                       >  1:274268/1‑221 (MQ=255)
                                                                                                                    aaacgcacaaGACCTTTAATGCTCACTCAATTGATACCACGGAAAAATGCAAAGCAAAACGCGGGGAACTGGCCGAGCGTGAAGGGGCGCGAATGAAAGAAAACAACACGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATa                                                                                                <  1:341194/219‑1 (MQ=255)
                                                                                                                            cagccttttATTGCTCACTCCATTGATACACCTAAAAAATCCGAACCCAACCCGGCCGAACTTGCCGGGCGGGAAGGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTa                                                                                          <  1:237149/220‑1 (MQ=255)
                                                                                                                              gACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCGGCGAAAGATATTCAAACTTGCACCGACTGTCATACAGGAAtgg                                                                                       >  1:257071/1‑225 (MQ=255)
                                                                                                                              gACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGCCAGCTGAAGGGGGCAGCGAAAGATAATCAATCATGCATGCACTTTCAAAAAGGAAATGc                                                                                      >  2:672853/1‑228 (MQ=255)
                                                                                                                               aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAACGAAAACAACTCGGCAACCTGCCGCTCCTGCCCTAACTACGATGCGATGGATCATGCGACGCAGCATCCTGAAACAGCCCGTCCGATGAAGGTGTCCGCGCAAGATAATCAATCCCGCCTCGCCCGTCATAAAGGAATTGcc                                                                                     >  2:341194/1‑228 (MQ=255)
                                                                                                                               aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGc                                                                                                                                                                                                                     >  2:50030/1‑101 (MQ=255)
                                                                                                                                ccTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCa                                                                                                                                                                                                                    <  1:50030/101‑1 (MQ=255)
                                                                                                                                  tttATTGCTCAATCCATGAATACCCCAAAAAAATTAGAAGCCAAACCCGCGGACTTTGCCGGGCGTGAAGGGGAGCAAAGGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCccac                                                                                  <  1:635779/228‑1 (MQ=255)
                                                                                                                                     aTTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAg                                                                               >  1:180386/1‑228 (MQ=255)
                                                                                                                                         cccccTCCATTGATACACCAGAAAAATTCAAAGCCAAACGCGCGGAACTGCCGGGGGGTGATGGGGGGAAAAAGAAAGAAAACACCCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  2:497711/224‑1 (MQ=255)
                                                                                                                                           caacacATTGAAACACCGAAAAAATTCGAAGCCAACCGCGCGGAACTGGCCGAGCGGGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  1:657249/221‑1 (MQ=255)
                                                                                                                                                     gacacacgaaaaaaaTTGAAGCCCAACAGCGCGGAACTGGCCGAGCTTGAAGGGGCTAGAAGGAAAGAAAACAACGCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGt                                                                <  2:468716/219‑1 (MQ=255)
                                                                                                                                                        cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGc                                                             >  1:457288/1‑227 (MQ=255)
                                                                                                                                                        cacCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGCCATCAAGGAATTGCCCACCGGTAACCGGatttaagaacc                                                             >  1:436147/1‑218 (MQ=255)
                                                                                                                                                          ccTGAACAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGTTAATCACTCCTGCATCGCCTGTCGTAAAGGTATGGCCCACCAGTTCCCGGCTAGgaggagacga                                                          >  1:306624/1‑224 (MQ=255)
                                                                                                                                                          ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAACGGGATTGCCCACCCGTTACCGGCTATGAg                                                                 >  1:468716/1‑221 (MQ=255)
                                                                                                                                                          ccTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGCTATGAGTAGCGGc                                                          >  1:481848/1‑228 (MQ=255)
                                                                                                                                                                                       aaCTTGCCGGCGGTGAATGGGCGAGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcg                                <  1:534916/226‑1 (MQ=255)
                                                                                                                                                                                          ctGCCGGGCGGAAAGGGGCGGAAATGAAAGAAACAAGCTCGCACACCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCa                           <  1:672853/227‑1 (MQ=255)
                                                                                                                                                                                                   gggAATGGGCACGAATGAAAAAAAACACATCGGCAACCTGCCGCCACTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                   <  1:282129/225‑1 (MQ=255)
                                                                                                                                                                                                      aaTGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                      >  1:176083/1‑157 (MQ=255)
                                                                                                                                                                                                            cgcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtgg           >  1:471753/1‑226 (MQ=255)
                                                                                                                                                                                                            cgcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGGCATAAAGGTATTGCCAACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGTTGCGCGCCAGTGCTAATGACCgtgg           >  1:658016/1‑226 (MQ=255)
                                                                                                                                                                                                             gtaaaggaaagaaaACAACCCGGCACCCTGCCGCCGCTGCCATAACTCCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCTAAAGATAATCAATCCTTCACCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtggtg         <  1:641149/221‑1 (MQ=255)
                                                                                                                                                                                                             gggAAGGAAAAAAAACGACGGGTCACCCGCCCGCCCCTGCCATAACTACGAGGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtggtg         <  2:273534/225‑1 (MQ=255)
                                                                                                                                                                                                                aaTAAAAGAAAACAACGCGGCACTCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAAGCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGaca      <  1:259712/227‑1 (MQ=255)
                                                                                                                                                                                                                  tGAAAAAAAACAACTCGGCAACCTGCAGCTCCTGACATAACTACGATGCGATGGATCAGGCGAAGCAGCAGCGTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg    <  1:337363/227‑1 (MQ=255)
                                                                                                                                                                                                                  ataaaaaaaaCAACTCGGCACCCTGCCCCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg    <  1:207114/225‑1 (MQ=255)
                                                                                                                                                                                                                   ccaagaaaACACCCCGGCCACCGCCCCGCCCTGCCATAACTACGATGCGATGGAACATGCGAAGCAGCATCCTGAAGCAGCCCGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCTGGTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCt  <  2:608662/226‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                       
ATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAACGCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCACGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCT  >  NC_000913/1058263‑1058700

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: