Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058,486 T→C 19.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,4860TC19.5% 99.6 / 6.8 58S135P (TCG→CCG) torCtrimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base T (31/11);  new base C (3/9);  total (34/24)
Fisher's exact test for biased strand distribution p-value = 4.85e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

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