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

ATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCT  >  NC_000913/1058263‑1058700
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                       
atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCTGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAAAAAAAAACTCgg                                                                                                                                                                                                                      >  1:657249/1‑227 (MQ=255)
atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCCCGGAACTTGCCGACCGTGAATGGG‑CGCGAATGAAAGAAAAAAACTCgg                                                                                                                                                                                                                      >  2:449064/1‑227 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGCTATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAATCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAAAAACTCGGAAAACTCCCg                                                                                                                                                                                                            >  1:515503/1‑228 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAATCGAGCGATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCCGGGGGGCACGGG‑CGGGAATGAAAGAAAAAAAATCGCCAACccgccg                                                                                                                                                                                                            >  1:641149/1‑228 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGGGAATGGC‑GGCGAATGAAAGAAAACAACTCGCCAAccgcccg                                                                                                                                                                                                            >  2:133233/1‑223 (MQ=255)
          tCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGCGCGTGAATGGG‑CGCGAATGAAAAAAAAAAACCCGGAAACCTGCCgg                                                                                                                                                                                                           >  1:282129/1‑227 (MQ=255)
           cGGTGCATTGCCAGACCGCCGCCGGCGTGCGAGCTGAAGGCCATGACTGTCATCTCCCGCAGGATATCCCAGGCATGGTGAAGCGAAAACTGGAAGCGAGCAATGATATCTACCAGTCCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCAAATGAAAGAAAACAACTCGGCAACCTGCCGCt                                                                                                                                                                                                          <  2:190932/228‑1 (MQ=255)
                    tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                 >  2:497711/1‑228 (MQ=255)
                      ccgaccGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGCATACCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACC‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCAt                                                                                                                                                                                                  <  2:488565/220‑1 (MQ=255)
                        gAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCGTAACTa                                                                                                                                                                                             >  2:208547/1‑228 (MQ=255)
                           cGCCTCCGGCGGGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                           <  2:87488/227‑1 (MQ=255)
                                        cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGC‑CTCGAATGAAAGAAAACAACCCGGCAACCCGCCGGCCCTGCCCTAACTACGAAGCAAAGGAACa                                                                                                                                                                              >  1:119057/1‑227 (MQ=255)
                                                    caggcgtGTCATAACCCGCCGGTTATTCCAGGCATGGGGAAGCAAAACGTGGAAGCGAGCAGTGATATCTGCCAGACCTTTATTGGCCACTCCAGTGATGCA‑CCTGAAAAATTCGAAGCCAAGC‑GCCCGGACCTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACCCGGCAACCTGCCTCTCCTGCCATAACTACGATGCGATGGATCATGCGAAgcagca                                                                                                                                                                  <  2:515503/221‑1 (MQ=255)
                                                            tCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGTTGCGATGGATCATGCGAAGCAGCATCCTGAagc                                                                                                                                                         >  2:561719/1‑228 (MQ=255)
                                                                    cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGGGAAGCAGCATCCTGAAGCAGCA‑CGTc                                                                                                                                                >  2:558378/1‑228 (MQ=255)
                                                                    cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAACCAGCATCCTGAAGCAGCC‑CGTc                                                                                                                                                >  2:205604/1‑228 (MQ=255)
                                                                                   gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGAAGGACCATTCGAAGCACCATCCTAAAGCAGCA‑CGTCAGATaaaagaggaag                                                                                                                                 >  1:635779/1‑220 (MQ=255)
                                                                                       ggTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCGGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGGTCATGCGAAGCAGCATCCTGAAGCAGAC‑CGTCAGATGAAGGTGGCAGCGa                                                                                                                              >  2:41716/1‑227 (MQ=255)
                                                                                        gTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGTG‑CGCGAATGCAAGAAAACAACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCCATGGATCATGCAAAGCAGCACCCT‑TCCCAACACCGCCAGATGAAGGGGGCGGCGaaa                                                                                                                            >  1:59041/1‑228 (MQ=255)
                                                                                        gTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGAGTGAATGGG‑CGCGAATGAAAGAACACAACTCGGCAACCTGCCGCTCCTGCCTTAACTACGATGCGGTGGATCATGGGAACCAGCATCCTGAAGAAGAA‑Cgccggcgcaggggggacgggaa                                                                                                                             >  1:134358/1‑206 (MQ=255)
                                                                                             gcgcAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGCTa                                                                                                                        >  2:413620/1‑227 (MQ=255)
                                                                                              cgcAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGAACATGCGAAGCAGCATCCTGAACAACCA‑CCCCAGAGAAAGGTGGCAGCGAAAGATaat                                                                                                                      >  1:337363/1‑228 (MQ=255)
                                                                                                   aCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGCCAGATGAAGGTGGCAGCGAAAGATaatcaatc                                                                                                                 >  2:24285/1‑228 (MQ=255)
                                                                                                   aCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGACCATGCGAAGCAGCAACCTGAAGCACCA‑CGCCAACGGAAGGGGGACGCGAAAAAAaataaaac                                                                                                                 >  1:207114/1‑226 (MQ=255)
                                                                                                             gacattaaatatcacaaaacttttATTGC‑CACTCCCGTTAAACACCCTAAAAAAACCGAAACAAACA‑GCGGGGAACTACCCG‑GCGGGAAGGGGCCGACAATGAAAGAAACAAACCCGGCAACCTCCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                       <  1:131781/207‑1 (MQ=255)
                                                                                                             atcAATGATATCTACCAGACCTTTATTGCACACTCCATTGACACA‑CCTGAAAAATTCGAAGCCAACC‑GCGCGGAACTTGCCGAGCGTGAAGGGG‑CGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                       <  2:119057/226‑1 (MQ=255)
                                                                                                             aGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATCCGAAGCAGCATCCTGAAGCAGCA‑CGGCAGATAAAGGTGGCACCGAAGAAAAACCAATCCTgaacagcc                                                                                                       >  2:274268/1‑221 (MQ=255)
                                                                                                                    aaacgcacaaGACCTTTAATGCTCACTCAATTGATACC‑ACGGAAAAATGCAAAGCAAAAC‑GCGGGGAACTGGCCGAGCGTGAAGGGG‑CGCGAATGAAAGAAAACAACACGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATa                                                                                                <  2:341194/219‑1 (MQ=255)
                                                                                                                            cagccttttATTGCTCACTCCATTGATACA‑CCTAAAAAATCCGAACCCAACC‑CGGCCGAACTTGCCGGGCGGGAAGGGG‑CGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTa                                                                                          <  2:237149/220‑1 (MQ=255)
                                                                                                                              gACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCGGCGAAAGATATTCAAACTTGCACCGACTGTCATACAGGAAtgg                                                                                       >  2:257071/1‑225 (MQ=255)
                                                                                                                              gACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGCCAGCTGAAGGGGGCAGCGAAAGATAATCAATCATGCATGCACTTTCAAAAAGGAAATGc                                                                                      >  1:672853/1‑228 (MQ=255)
                                                                                                                               aCCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAACGAAAACAACTCGGCAACCTGCCGCTCCTGCCCTAACTACGATGCGATGGATCATGCGACGCAGCATCCTGAAACAGCC‑CGTCCGATGAAGGTGTCCGCGCAAGATAATCAATCCCGCCTCGCCCGTCATAAAGGAATTGcc                                                                                     >  1:341194/1‑228 (MQ=255)
                                                                                                                               aCCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGc                                                                                                                                                                                                                     >  1:50030/1‑101 (MQ=255)
                                                                                                                                ccTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCa                                                                                                                                                                                                                    <  2:50030/101‑1 (MQ=255)
                                                                                                                                  tttATTGCTCAATCCATGAATACC‑CCAAAAAAATTAGAAGCCAAAC‑CCGCGGACTTTGCCGGGCGTGAAGGGG‑AGCAAAGGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCccac                                                                                  <  2:635779/228‑1 (MQ=255)
                                                                                                                                     aTTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAg                                                                               >  2:180386/1‑228 (MQ=255)
                                                                                                                                         cccccTCCATTGATACA‑CCAGAAAAATTCAAAGCCAAAC‑GCGCGGAACTGCCGGGGGGTGATGGGG‑GGAAAAAGAAAGAAAACACCCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  1:497711/224‑1 (MQ=255)
                                                                                                                                          cgcacaccatgaatccccttGAAAAAGTCAAAGCAAACC‑GCGCGGACCTTCCCGAGCGTGAGTGGG‑CGCGAAGAAAAGAAAACAACCCGGCACCCTGCCGCTCCTGCCATAACTGCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  1:180386/209‑1 (MQ=255)
                                                                                                                                           caacacATTGAAACA‑CCGAAAAAATTCGAAGCCAACC‑GCGCGGAACTGGCCGAGCGGGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  2:657249/221‑1 (MQ=255)
                                                                                                                                                      gacacacgaaaaaaaTTGAAGCCCAACAGCGCGGAACTGGCCGAGCTTGAAGGGG‑CTAGAAGGAAAGAAAACAACGCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGt                                                                <  1:468716/219‑1 (MQ=255)
                                                                                                                                                        ca‑cCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGc                                                             >  2:457288/1‑227 (MQ=255)
                                                                                                                                                        ca‑cCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGCCATCAAGGAATTGCCCACCGGTAACCGGatttaagaacc                                                             >  2:436147/1‑218 (MQ=255)
                                                                                                                                                           ccTGAACAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGTTAATCACTCCTGCATCGCCTGTCGTAAAGGTATGGCCCACCAGTTCCCGGCTAGgaggagacga                                                          >  2:306624/1‑224 (MQ=255)
                                                                                                                                                           ccTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAACGGGATTGCCCACCCGTTACCGGCTATGAg                                                                 >  2:468716/1‑221 (MQ=255)
                                                                                                                                                           ccTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGCTATGAGTAGCGGc                                                          >  2:481848/1‑228 (MQ=255)
                                                                                                                                                                                        aaCTTGCCGGCGGTGAATGGG‑CGAGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcg                                <  2:534916/226‑1 (MQ=255)
                                                                                                                                                                                           ctGCCGGGCGGAAAGGGG‑CGGAAATGAAAGAAACAAGCTCGCACACCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCa                           <  2:672853/227‑1 (MQ=255)
                                                                                                                                                                                                    gggAATGGG‑CACGAATGAAAAAAAACACATCGGCAACCTGCCGCCACTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                   <  2:282129/225‑1 (MQ=255)
                                                                                                                                                                                                     gggaagggggggaataaaaaaaaaaCACCCGGGAAACCCGCCCTCCCTGCCAAAACTACGATGCGATGGATCAAGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                   <  1:24285/210‑1 (MQ=255)
                                                                                                                                                                                                       aaTGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                      >  2:176083/1‑157 (MQ=255)
                                                                                                                                                                                                              cgcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtgg           >  2:471753/1‑226 (MQ=255)
                                                                                                                                                                                                              cgcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGGCATAAAGGTATTGCCAACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGTTGCGCGCCAGTGCTAATGACCgtgg           >  2:658016/1‑226 (MQ=255)
                                                                                                                                                                                                               gtaaaggaaagaaaACAACCCGGCACCCTGCCGCCGCTGCCATAACTCCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCTAAAGATAATCAATCCTTCACCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtggtg         <  2:641149/221‑1 (MQ=255)
                                                                                                                                                                                                               gggAAGGAAAAAAAACGACGGGTCACCCGCCCGCCCCTGCCATAACTACGAGGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtggtg         <  1:273534/225‑1 (MQ=255)
                                                                                                                                                                                                                  aaTAAAAGAAAACAACGCGGCACTCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAAGCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGaca      <  2:259712/227‑1 (MQ=255)
                                                                                                                                                                                                                    tGAAAAAAAACAACTCGGCAACCTGCAGCTCCTGACATAACTACGATGCGATGGATCAGGCGAAGCAGCAGCGTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg    <  2:337363/227‑1 (MQ=255)
                                                                                                                                                                                                                    ataaaaaaaaCAACTCGGCACCCTGCCCCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg    <  2:207114/225‑1 (MQ=255)
                                                                                                                                                                                                                     ccaagaaaACACCCCGGCCACCGCCCCGCCCTGCCATAACTACGATGCGATGGAACATGCGAAGCAGCATCCTGAAGCAGCC‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCTGGTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCt  <  1:608662/226‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                       
ATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACA‑CCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGG‑CGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCT  >  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: ‑