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

ATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCT  >  NC_000913/1058263‑1058700
                                                                                                                                                                                                                                |                                                                                                                                                                                                                       
atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCTGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAAAAAAAAACTCgg                                                                                                                                                                                                                      >  2:657249/1‑227 (MQ=255)
atatAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCCCGGAACTTGCCGACCGTGAATGGGCGCGAATGAAAGAAAAAAACTCgg                                                                                                                                                                                                                      >  1:449064/1‑227 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGCTATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAATCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAAAAACTCGGAAAACTCCCg                                                                                                                                                                                                            >  2:515503/1‑228 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAATCGAGCGATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCCGGGGGGCACGGGCGGGAATGAAAGAAAAAAAATCGCCAACccgccg                                                                                                                                                                                                            >  2:641149/1‑228 (MQ=255)
         gTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGGGAATGGCGGCGAATGAAAGAAAACAACTCGCCAAccgcccg                                                                                                                                                                                                            >  1:133233/1‑223 (MQ=255)
          tCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGCGCGTGAATGGGCGCGAATGAAAAAAAAAAACCCGGAAACCTGCCgg                                                                                                                                                                                                           >  2:282129/1‑227 (MQ=255)
           cGGTGCATTGCCAGACCGCCGCCGGCGTGCGAGCTGAAGGCCATGACTGTCATCTCCCGCAGGATATCCCAGGCATGGTGAAGCGAAAACTGGAAGCGAGCAATGATATCTACCAGTCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCAAATGAAAGAAAACAACTCGGCAACCTGCCGCt                                                                                                                                                                                                          <  1:190932/228‑1 (MQ=255)
                    tCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATa                                                                                                                                                                                                 >  1:497711/1‑228 (MQ=255)
                      ccgaccGCCTCCGGCGTGCGAGCTGAATGCCATGCCTGGCATACCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACCCCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCAt                                                                                                                                                                                                  <  1:488565/220‑1 (MQ=255)
                        gAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCGTAACTa                                                                                                                                                                                             >  1:208547/1‑228 (MQ=255)
                           cGCCTCCGGCGGGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTAcg                                                                                                                                                                                           <  1:87488/227‑1 (MQ=255)
                                        cGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGCCTCGAATGAAAGAAAACAACCCGGCAACCCGCCGGCCCTGCCCTAACTACGAAGCAAAGGAACa                                                                                                                                                                              >  2:119057/1‑227 (MQ=255)
                                                    caggcgtGTCATAACCCGCCGGTTATTCCAGGCATGGGGAAGCAAAACGTGGAAGCGAGCAGTGATATCTGCCAGACCTTTATTGGCCACTCCAGTGATGCACCTGAAAAATTCGAAGCCAAGC‑GCCCGGACCTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCTCTCCTGCCATAACTACGATGCGATGGATCATGCGAAgcagca                                                                                                                                                                  <  1:515503/221‑1 (MQ=255)
                                                            tCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGTTGCGATGGATCATGCGAAGCAGCATCCTGAagc                                                                                                                                                         >  1:561719/1‑228 (MQ=255)
                                                                    cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGGGAAGCAGCATCCTGAAGCAGCA‑CGTc                                                                                                                                                >  1:558378/1‑228 (MQ=255)
                                                                    cgccgGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAACCAGCATCCTGAAGCAGCC‑CGTc                                                                                                                                                >  1:205604/1‑228 (MQ=255)
                                                                                   gCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGAAGGACCATTCGAAGCACCATCCTAAAGCAGCA‑CGTCAGATaaaagaggaag                                                                                                                                 >  2:635779/1‑220 (MQ=255)
                                                                                       ggTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCGGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGGTCATGCGAAGCAGCATCCTGAAGCAGAC‑CGTCAGATGAAGGTGGCAGCGa                                                                                                                              >  1:41716/1‑227 (MQ=255)
                                                                                        gTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGTGCGCGAATGCAAGAAAACAACCCGGCAACCGGCCGCTCCTGCCATAACTACGATGCCATGGATCATGCAAAGCAGCACCCT‑TCCCAACACCGCCAGATGAAGGGGGCGGCGaaa                                                                                                                            >  2:59041/1‑228 (MQ=255)
                                                                                             gcgcAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGCTa                                                                                                                        >  1:413620/1‑227 (MQ=255)
                                                                                              cgcAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGAACATGCGAAGCAGCATCCTGAACAACCA‑CCCCAGAGAAAGGTGGCAGCGAAAGATaat                                                                                                                      >  2:337363/1‑228 (MQ=255)
                                                                                                   aCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGCCAGATGAAGGTGGCAGCGAAAGATaatcaatc                                                                                                                 >  1:24285/1‑228 (MQ=255)
                                                                                                   aCTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGACCATGCGAAGCAGCAACCTGAAGCACCA‑CGCCAACGGAAGGGGGACGCGAAAAAAaataaaac                                                                                                                 >  2:207114/1‑226 (MQ=255)
                                                                                                             atcAATGATATCTACCAGACCTTTATTGCACACTCCATTGACACACCTGAAAAATTCGAAGCCAACC‑GCGCGGAACTTGCCGAGCGTGAAGGGGCGCGAATGAAAGAAACCAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGAc                                                                                                       <  1:119057/226‑1 (MQ=255)
                                                                                                             aGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATCCGAAGCAGCATCCTGAAGCAGCA‑CGGCAGATAAAGGTGGCACCGAAGAAAAACCAATCCTgaacagcc                                                                                                       >  1:274268/1‑221 (MQ=255)
                                                                                                                    aaacgcacaaGACCTTTAATGCTCACTCAATTGATACCACGGAAAAATGCAAAGCAAAAC‑GCGGGGAACTGGCCGAGCGTGAAGGGGCGCGAATGAAAGAAAACAACACGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATa                                                                                                <  1:341194/219‑1 (MQ=255)
                                                                                                                            cagccttttATTGCTCACTCCATTGATACACCTAAAAAATCCGAACCCAACC‑CGGCCGAACTTGCCGGGCGGGAAGGGGCGCGAATGAAAGAAAACAACCCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTa                                                                                          <  1:237149/220‑1 (MQ=255)
                                                                                                                              gACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCGGCGAAAGATATTCAAACTTGCACCGACTGTCATACAGGAAtgg                                                                                       >  1:257071/1‑225 (MQ=255)
                                                                                                                              gACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGCCAGCTGAAGGGGGCAGCGAAAGATAATCAATCATGCATGCACTTTCAAAAAGGAAATGc                                                                                      >  2:672853/1‑228 (MQ=255)
                                                                                                                               aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAACGAAAACAACTCGGCAACCTGCCGCTCCTGCCCTAACTACGATGCGATGGATCATGCGACGCAGCATCCTGAAACAGCC‑CGTCCGATGAAGGTGTCCGCGCAAGATAATCAATCCCGCCTCGCCCGTCATAAAGGAATTGcc                                                                                     >  2:341194/1‑228 (MQ=255)
                                                                                                                               aCCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGc                                                                                                                                                                                                                     >  2:50030/1‑101 (MQ=255)
                                                                                                                                ccTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCa                                                                                                                                                                                                                    <  1:50030/101‑1 (MQ=255)
                                                                                                                                  tttATTGCTCAATCCATGAATACCCCAAAAAAATTAGAAGCCAAAC‑CCGCGGACTTTGCCGGGCGTGAAGGGGAGCAAAGGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCccac                                                                                  <  1:635779/228‑1 (MQ=255)
                                                                                                                                     aTTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAg                                                                               >  1:180386/1‑228 (MQ=255)
                                                                                                                                         cccccTCCATTGATACACCAGAAAAATTCAAAGCCAAAC‑GCGCGGAACTGCCGGGGGGTGATGGGGGGAAAAAGAAAGAAAACACCCCGGCACCCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  2:497711/224‑1 (MQ=255)
                                                                                                                                           caacacATTGAAACACCGAAAAAATTCGAAGCCAACC‑GCGCGGAACTGGCCGAGCGGGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTAc                                                                           <  1:657249/221‑1 (MQ=255)
                                                                                                                                                     gacacacgaaaaaaaTTGAAGCCCAACAGCGCGGAACTGGCCGAGCTTGAAGGGGCTAGAAGGAAAGAAAACAACGCGGCAACCGGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGt                                                                <  2:468716/219‑1 (MQ=255)
                                                                                                                                                        cacCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGc                                                             >  1:457288/1‑227 (MQ=255)
                                                                                                                                                        cacCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGCCATCAAGGAATTGCCCACCGGTAACCGGatttaagaacc                                                             >  1:436147/1‑218 (MQ=255)
                                                                                                                                                          ccTGAACAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGTTAATCACTCCTGCATCGCCTGTCGTAAAGGTATGGCCCACCAGTTCCCGGCTAGgaggagacga                                                          >  1:306624/1‑224 (MQ=255)
                                                                                                                                                          ccTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAACGGGATTGCCCACCCGTTACCGGCTATGAg                                                                 >  1:468716/1‑221 (MQ=255)
                                                                                                                                                          ccTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGAATTGCCCACCAGTTACCGGCTATGAGTAGCGGc                                                          >  1:481848/1‑228 (MQ=255)
                                                                                                                                                                                       aaCTTGCCGGCGGTGAATGGGCGAGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTgcg                                <  1:534916/226‑1 (MQ=255)
                                                                                                                                                                                          ctGCCGGGCGGAAAGGGGCGGAAATGAAAGAAACAAGCTCGCACACCTGCCGCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCa                           <  1:672853/227‑1 (MQ=255)
                                                                                                                                                                                                   gggAATGGGCACGAATGAAAAAAAACACATCGGCAACCTGCCGCCACTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAAt                   <  1:282129/225‑1 (MQ=255)
                                                                                                                                                                                                      aaTGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGc                                                                                      >  1:176083/1‑157 (MQ=255)
                                                                                                                                                                                                            cgcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtgg           >  1:471753/1‑226 (MQ=255)
                                                                                                                                                                                                            cgcgAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGGCATAAAGGTATTGCCAACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGTTGCGCGCCAGTGCTAATGACCgtgg           >  1:658016/1‑226 (MQ=255)
                                                                                                                                                                                                             gtaaaggaaagaaaACAACCCGGCACCCTGCCGCCGCTGCCATAACTCCGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCTAAAGATAATCAATCCTTCACCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtggtg         <  1:641149/221‑1 (MQ=255)
                                                                                                                                                                                                             gggAAGGAAAAAAAACGACGGGTCACCCGCCCGCCCCTGCCATAACTACGAGGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAgtggtg         <  2:273534/225‑1 (MQ=255)
                                                                                                                                                                                                                aaTAAAAGAAAACAACGCGGCACTCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAAGCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGaca      <  1:259712/227‑1 (MQ=255)
                                                                                                                                                                                                                  tGAAAAAAAACAACTCGGCAACCTGCAGCTCCTGACATAACTACGATGCGATGGATCAGGCGAAGCAGCAGCGTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg    <  1:337363/227‑1 (MQ=255)
                                                                                                                                                                                                                  ataaaaaaaaCAACTCGGCACCCTGCCCCCCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCGGCTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACg    <  1:207114/225‑1 (MQ=255)
                                                                                                                                                                                                                   ccaagaaaACACCCCGGCCACCGCCCCGCCCTGCCATAACTACGATGCGATGGAACATGCGAAGCAGCATCCTGAAGCAGCC‑CGTCAGATGAAGGTGGCAGCGAAAGATAATCAATCCTGCATCGACTGTCATAAAGGTATTGCCCACCAGTTACCGGATATGAGTAGCTGGTTCCGTAAGCAGTTCGATGAGCTGCGCGCCAGTGCTAATGACAGTGGTGACACGCt  <  2:608662/226‑1 (MQ=255)
                                                                                                                                                                                                                                |                                                                                                                                                                                                                       
ATATAAACAGTCGGTGCATTTCCAGAACGCCTCCGGCGTGCGAGCTGAATGCCATGACTGTCATATCCCGCCGGATATTCCAGGCATGGTGAAGCGCAAACTGGAAGCGAGCAATGATATCTACCAGACCTTTATTGCTCACTCCATTGATACACCTGAAAAATTCGAAGCCAAAC‑GCGCGGAACTTGCCGAGCGTGAATGGGCGCGAATGAAAGAAAACAACTCGGCAACCTGCCGCTCCTGCCATAACTACGATGCGATGGATCATGCGAAGCAGCATCCTGAAGCAGCA‑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: ‑