Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,306,533 A→G 13.8% P368P (CCT→CCC lmrS ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,306,5330AG13.8% 66.1 / 6.0 29P368P (CCT→CCClmrShypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base A (15/10);  new base G (2/2);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01

GTATATTAACTGCGATATACTAAGTGCCCCAATAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACCACATATCAA  >  CP000730/2306398‑2306666
                                                                                                                                       |                                                                                                                                     
gTATATTAACTGCGATATACTAAGTGCCCCAATAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCGGGTAt                                                                                                                                  >  2:254149/1‑141 (MQ=255)
gTATATTAACTGCGATATACTAAGTGCCCCAATAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCGGGTAt                                                                                                                                  >  2:254102/1‑141 (MQ=255)
              atatACTAAGTGCCCCAATAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTa                                                                                                                    <  2:8444/141‑1 (MQ=255)
                  aCTAAGTGCCCCAATAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATAGAGAAGGTAAAAGCAGATCGGCAAAACCGATAATAATCGGTGGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTaaaaa                                                                                                                >  1:158707/1‑141 (MQ=255)
                              aaTAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAgg                                                                                                    <  2:87575/141‑1 (MQ=255)
                                                   aTTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACaaaaaa                                                                               >  2:244555/1‑141 (MQ=255)
                                                            acCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTa                                                                      <  1:242436/141‑1 (MQ=255)
                                                            acCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTa                                                                      <  1:179550/141‑1 (MQ=255)
                                                                 cAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCg                                                                 <  2:158707/141‑1 (MQ=255)
                                                                       atatAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCaataa                                                           >  2:241307/1‑141 (MQ=255)
                                                                          taGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATgg                                                        >  1:141994/1‑141 (MQ=255)
                                                                          aaGAAGGTAAAAACATATCGGAAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATAAACTCTAAAAATGGAACAAAAGGTTTAAAAAAAAAAAAAAAAAAGGGCGTGTAGTTCGCAATAATgg                                                        <  2:126943/140‑1 (MQ=255)
                                                                           aGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGATTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGt                                                       >  2:167341/1‑141 (MQ=255)
                                                                             aaGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAg                                                     >  2:126349/1‑141 (MQ=255)
                                                                                      aCATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTAc                                            >  2:84643/1‑141 (MQ=255)
                                                                                         aatCGGCAAAACCCATAATAATCGGTTGAAAAGCCTGTGCCGCTTCAAGTATTTTCAATAACCCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCaa                                         <  1:102654/140‑1 (MQ=255)
                                                                                                aaaaCCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCAt                                  >  2:74242/1‑141 (MQ=255)
                                                                                                  aaCCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATta                                >  2:209837/1‑141 (MQ=255)
                                                                                                   aCCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATtac                               >  2:47862/1‑141 (MQ=255)
                                                                                                          ataatCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGc                        >  1:155134/1‑141 (MQ=255)
                                                                                                          ataatCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGc                        >  2:55925/1‑141 (MQ=255)
                                                                                                           taatCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCa                       >  2:82778/1‑141 (MQ=255)
                                                                                                                  ttGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAAc                >  2:212676/1‑141 (MQ=255)
                                                                                                                     agaAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACaaa             <  2:262686/141‑1 (MQ=255)
                                                                                                                        aGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACca          <  2:98259/141‑1 (MQ=255)
                                                                                                                           cTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACCACat       <  1:2059/141‑1 (MQ=255)
                                                                                                                              tgCCGCTTCGGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACCACATATc    <  1:254149/141‑1 (MQ=255)
                                                                                                                              tgCCGCTTCGGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACCACATATc    <  1:254102/141‑1 (MQ=255)
                                                                                                                                ccGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACCACATATCaa  >  2:235368/1‑141 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GTATATTAACTGCGATATACTAAGTGCCCCAATAACAAAACGAGTTATATCATTTTGTACACCTTCAATTAATATAGAAGGTAAAAACATATCGGCAAAACCGATAATAATCGTTTGAGAAGCCTGTGCCGCTTCAGGTATTTGCATTAACTCTAAAAATGGAACAAAAGGTTTACCCAATATGACAAAAAAGGGCGTGTAGTTCGCAATAATGGTAGCAATAGTACCAATACTCATTACTACAGGCAAAATAACAAACCACATATCAA  >  CP000730/2306398‑2306666

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