Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 757,227 C→A 14.8% L441M (CTG→ATG)  sdhA → succinate dehydrogenase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913757,2270CA14.8% 74.4 / 6.7 27L441M (CTG→ATG) sdhAsuccinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (3/20);  new base A (2/2);  total (5/22)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGA  >  NC_000913/757079‑757360
                                                                                                                                                    |                                                                                                                                     
tgtATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTctct                                                                                                                                      <  2:114974/150‑1 (MQ=255)
        tACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTATGGATCGcc                                                                                                                              >  2:39056/1‑150 (MQ=255)
        tACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTATGGATCGcc                                                                                                                              >  2:39058/1‑150 (MQ=255)
                gcTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGc                                                                                                                      <  2:194816/150‑1 (MQ=255)
                gcTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGc                                                                                                                      <  2:186782/150‑1 (MQ=255)
                gcTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGc                                                                                                                      <  2:194817/150‑1 (MQ=255)
                gcTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGc                                                                                                                      <  2:194819/150‑1 (MQ=255)
                gcTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGc                                                                                                                      <  2:205609/150‑1 (MQ=255)
                    aCCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGa                                                                                                                  <  2:134003/150‑1 (MQ=255)
                    aCCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGa                                                                                                                  <  2:140683/150‑1 (MQ=255)
                    aCCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGa                                                                                                                  <  2:140684/150‑1 (MQ=255)
                                                  tggtggTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTgg                                                                                    <  2:2382/150‑1 (MQ=255)
                                                           ttGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTa                                                                           <  2:11472/150‑1 (MQ=255)
                                                                                 tctgcaAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATgca                                                     <  2:52619/150‑1 (MQ=255)
                                                                                        gagTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAAc                                              >  1:234846/1‑150 (MQ=255)
                                                                                        gagTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAAc                                              >  1:234838/1‑150 (MQ=255)
                                                                                                              cACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTg                        <  1:26849/150‑1 (MQ=255)
                                                                                                              cACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTg                        <  1:26832/150‑1 (MQ=255)
                                                                                                              cACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTg                        <  1:26829/150‑1 (MQ=255)
                                                                                                                cTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGAt                      <  2:15181/150‑1 (MQ=255)
                                                                                                                cTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTATGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGAt                      <  1:39058/150‑1 (MQ=255)
                                                                                                                cTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTATGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGAt                      <  1:39056/150‑1 (MQ=255)
                                                                                                                  gcgcGAGGCCAGCGCGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACACTAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGc                    <  2:17971/150‑1 (MQ=255)
                                                                                                                             gCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATg                <  1:77178/143‑1 (MQ=255)
                                                                                                                             gCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATg                >  2:77178/1‑143 (MQ=255)
                                                                                                                                  tCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCtt    <  1:124738/150‑1 (MQ=255)
                                                                                                                                    tGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGa  <  2:167891/150‑1 (MQ=255)
                                                                                                                                    tGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGa  <  2:167872/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                     
TGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAACTTCTCGGTCTTCCGTGAAGGTGATGCGATGGCGAAAGGGCTTGA  >  NC_000913/757079‑757360

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: