Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,442,061 A→G V39A (GTA→GCA)  msrA ← methionine sulfoxide reductase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,442,0610AG96.2% 85.7 / ‑6.5 26V39A (GTA→GCA) msrAmethionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (12/13);  minor base T (0/1);  total (12/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCA  >  NC_000913/4441926‑4442175
                                                                                                                                       |                                                                                                                  
cGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCAt                                                                                                                 <  2:81856/139‑1 (MQ=255)
    aaGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGt                                                                                                             <  1:392496/139‑1 (MQ=255)
           aTTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGATCAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAg                                                                                                      <  1:358550/139‑1 (MQ=255)
             ttGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTg                                                                                                    <  2:72465/139‑1 (MQ=255)
                                              tGCTGTAAACGCCGGGTAACTGCCAGAACAGACCCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCAt                                                                   >  2:66804/1‑139 (MQ=255)
                                                cTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCg                                                                 <  1:275666/139‑1 (MQ=255)
                                                   tAAACGCCGGGTAACTGCCAGAACAGACTCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCgggg                                                              >  1:363305/1‑139 (MQ=255)
                                                        gCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGt                                                                                               >  2:42863/1‑101 (MQ=255)
                                                        gCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGt                                                                                               <  1:42863/101‑1 (MQ=255)
                                                          cGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTg                                                                                                        <  2:457241/90‑1 (MQ=255)
                                                          cGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTg                                                                                                        >  1:457241/1‑90 (MQ=255)
                                                                aCTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTTCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg                                                 <  2:363305/139‑1 (MQ=255)
                                                                  tGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGAcc                                                                                          <  2:120472/96‑1 (MQ=255)
                                                                  tGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGAcc                                                                                          >  1:120472/1‑96 (MQ=255)
                                                                                        cccAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAg                         >  2:130614/1‑139 (MQ=255)
                                                                                             aaaCAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATg                                                                                     <  2:399288/74‑1 (MQ=255)
                                                                                             aaaCAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATg                                                                                     >  1:399288/1‑74 (MQ=255)
                                                                                                  acCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTa               >  2:332904/1‑139 (MQ=255)
                                                                                                  acCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTa               >  2:168711/1‑139 (MQ=255)
                                                                                                   cccATCGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTAt              <  1:168711/139‑1 (MQ=255)
                                                                                                        cGCAAAAATGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGCGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAAt         >  2:195026/1‑139 (MQ=255)
                                                                                                               aTGGCAATCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCa  <  1:426669/139‑1 (MQ=255)
                                                                                                                    aaTCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg                                                 <  2:271425/87‑1 (MQ=255)
                                                                                                                    aaTCTCCATTCCGTCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCcagg                                                 >  1:271425/1‑87 (MQ=255)
                                                                                                                                 tCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCa  <  1:303365/121‑1 (MQ=255)
                                                                                                                                 tCAGGTGCATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTATTTATCAAATAAACTCa  >  2:303365/1‑121 (MQ=255)
                                                                                                                                       |                                                                                                                  
CGATAAGTCGGATTTGGCGTATAGCCGCCGGTATAGCCTGCGGCGGTGCTGTAAACGCCGGGTAACTGCCAGAACAGACGCTCCACACCCCAGAAACAACCCATCGCAAAAATGGCAATCTCCATTCCGTCAGGTACATTGGTCATTGAGTGACCGTTGACCGCATGCAGCGTGGCTACGGGCATCGGGGTGTTACGTCCAGGCAGGGCATCGGCGGGGGAAACCAGATGCTTTTTATCAAATAAACTCA  >  NC_000913/4441926‑4442175

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