Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 432,015 A→G 15.4% F101S (TTC→TCC)  PP_0354 ← CBS domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947432,0150AG15.4% 51.2 / 3.5 26F101S (TTC→TCC) PP_0354CBS domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (9/13);  new base G (2/2);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

ACCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGA  >  NC_002947/431870‑432144
                                                                                                                                                 |                                                                                                                                  
acCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGGAgtg                                                                                                                                >  2:326837/1‑148 (MQ=255)
 ccGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCAc                                                                                                                               >  1:218230/1‑150 (MQ=255)
            aaaaTCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGt                                                                                                                    <  2:72451/150‑1 (MQ=255)
                  ccGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGGAGTGCTTCCCGCTggcagccct                                                                                                              <  1:326837/150‑8 (MQ=255)
                    gAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGa                                                                                                            <  1:365810/150‑1 (MQ=255)
                      aaCACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAAc                                                                                                          <  2:264086/150‑1 (MQ=255)
                                         acacGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGtt                                                                                       <  1:135336/150‑1 (MQ=255)
                                           acGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGAGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGAACTCGAACGTGGTTTCGACACCCGGTTTg                                                                                     <  2:286303/150‑1 (MQ=255)
                                                   gAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTg                                                                             >  2:469185/1‑150 (MQ=255)
                                                   gAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTACCAGCGGGGAGTGCTTCCCGCTgg                                                                                                                     >  2:459021/1‑110 (MQ=255)
                                                   gAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTACCAGCGGGGAGTGCTTCCCGCTgg                                                                                                                     <  1:459021/110‑1 (MQ=255)
                                                       gCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAGCGTGGTTTCGACACCCGGTTTGGCCAGGTGCTgg                                                                         <  1:297000/150‑1 (MQ=255)
                                                             ttGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCg                                                                   >  2:508325/1‑150 (MQ=255)
                                                                           gcagAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGcc                                                     <  1:469185/150‑1 (MQ=255)
                                                                               aaGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCAGc                                                 <  1:508325/150‑1 (MQ=255)
                                                                                  gtgtCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGt                                                                                                                    <  2:60895/80‑1 (MQ=255)
                                                                                  gtgtCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGt                                                                                                                    >  1:60895/1‑80 (MQ=255)
                                                                                   tgtCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGa                                             <  1:45485/150‑1 (MQ=255)
                                                                                       ttCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGAt                                         >  2:273346/1‑150 (MQ=255)
                                                                                           ccGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGttt                                                                                      <  1:407205/101‑1 (MQ=255)
                                                                                           ccGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGttt                                                                                      >  2:407205/1‑101 (MQ=255)
                                                                                                   gTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTc                            <  2:218230/150‑1 (MQ=255)
                                                                                                       tCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTTCCAGCGTGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTACTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGAcc                        <  1:192537/150‑1 (MQ=255)
                                                                                                       tCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGACGTAGAAGTGCTCGAcc                        >  2:264539/1‑150 (MQ=255)
                                                                                                                     ccGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAgg          >  1:59292/1‑150 (MQ=255)
                                                                                                                             ccAGGAGGGCTGACAGCAGGAAGATCACCCCGCTGAAGATCGCAAGCGTGGTTTCGATACCCGGTTTCGCCAGCTGCTGGCGCTCGTCAACCCCGCCACTCTGCCTTAGAAGTCGGACGTGAGCGGCGGGCCCGTCGGCAGGTTTGACGa  >  2:364902/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                  
ACCGCGCAGGGCAAAATCCCGAAACACATCAGACATCGAGAACACGCGGATGAATGCGGCCTTGTTCAGTTGCAGGCAGAAGGTGTCTTCCCCGGCCAGGTGCTCGGTGCGGGTTGCCCGCTCACCCAGCAGGGCTGCCAGCGGGAAGCACTCCCCGCTGGTGATCTCGAACGTGGTTTCGACACCCGGTTTGGCCAGGTGCTGGCGCTCGCCAACCACGCGGCCCTGCTTGACGATGTAGAAGTGCTCGACCGGCCCGTCGGCAGGTTTGACGA  >  NC_002947/431870‑432144

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