Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,251,102 C→T 12.9% A228T (GCA→ACA)  hmrA ← M20D subfamily peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,251,1020CT12.9% 71.5 / 4.2 31A228T (GCA→ACA) hmrAM20D subfamily peptidase
Reads supporting (aligned to +/- strand):  ref base C (13/14);  new base T (2/2);  total (15/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.22e-01

TTCATAATCACAACCAGTCTGTATCGCAGCTCCACGTGCGATTTGATTTACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCAGAGGCATGAGC  >  CP000730/2250965‑2251240
                                                                                                                                         |                                                                                                                                          
tccATAATCACAACCAGTCTGTATCGCAGCTCCACGTGCGATTTGATTTACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCttt                                                                                                                                         <  1:321717/139‑1 (MQ=255)
 tcatcatcaCCACCAGTCTGTATCGCCGCTCCACGTGCGATTTGATTTACTTTTTCTCTAAACACATCCAATACTTTACGCGTCATTGCCCCAGTATAAAAACGAGCATGTTTATATTCAGAAATAATATTAGCTGCTTTc                                                                                                                                        <  2:184669/141‑1 (MQ=255)
             ccAGTCTGTATCGCAGCTCCACGTGCGATTTGATTTACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGTTTTCCCGCCATCTaaa                                                                                                                            >  2:53753/1‑141 (MQ=255)
             ccAGTCTGTATCGCAGCTCCACGTGCGATTTGATTTACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGTTTTCCCGCCATCTaaa                                                                                                                            >  2:236900/1‑141 (MQ=255)
                           aGCTCCACGTGCGATTTGATTTACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGcac                                                                                                              <  1:111320/141‑1 (MQ=255)
                                            tatttACTTTTTCTGTTAATATATCCTATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTaa                                                                                             <  2:255909/140‑1 (MQ=255)
                                                tACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCAAACGTTGATCTTTTTTAatatt                                                                                         >  2:9171/1‑140 (MQ=255)
                                                   tttttCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGACCACGTTGATATTTTTTCATAtgttg                                                                                      >  1:12481/1‑141 (MQ=255)
                                                     tttCTGTTAATATATTCAATTCTTTACCCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTc                                                                                    <  1:43860/141‑1 (MQ=255)
                                                       tCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGt                                                                                  >  2:371772/1‑141 (MQ=255)
                                                        cTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTa                                                                                 <  1:299033/141‑1 (MQ=255)
                                                          gTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGt                                                                               >  1:56393/1‑141 (MQ=255)
                                                             aaTATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACg                                                                                                           <  1:195755/110‑1 (MQ=255)
                                                             aaTATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACg                                                                                                           >  2:195755/1‑110 (MQ=255)
                                                                  atCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTac                                                                       <  1:232472/141‑1 (MQ=255)
                                                                      aaTTCTTTACGCGTCATTTCTCTAGTATCAAAACGGTCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCa                                                                   <  2:307647/141‑1 (MQ=255)
                                                                          cTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTATCTGCTTTCCCGCCATCTAAAATCACACCTTGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTaa                                                               >  1:215664/1‑141 (MQ=255)
                                                                                                         aGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTACTGCATTTaa                                >  2:123903/1‑141 (MQ=255)
                                                                                                          gCATGTGTATAGTCTGGAATAATATTAGCTGTTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAGATAACTAATCATAGCGTCTAATGCATTTAAc                               <  1:236900/141‑1 (MQ=255)
                                                                                                          gCATGTGTATAGTCTGGAATAATATTAGCTGTTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAAc                               <  1:53753/141‑1 (MQ=255)
                                                                                                            aTGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAACGTCTAATGCATTTAACGc                             <  1:371772/141‑1 (MQ=255)
                                                                                                                    aGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATAAGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTg                     >  1:135875/1‑141 (MQ=255)
                                                                                                                     gTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGc                    >  2:75468/1‑141 (MQ=255)
                                                                                                                      tCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCa                   <  1:116360/141‑1 (MQ=255)
                                                                                                                         ggAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCAttt                >  2:167771/1‑141 (MQ=255)
                                                                                                                         ggAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCAttt                >  2:13853/1‑141 (MQ=255)
                                                                                                                          gAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTCTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCAtttt               <  1:214097/141‑1 (MQ=255)
                                                                                                                            ataataTTAGCTGCTTTCCCGCCATCTAAAATCACACCATTCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCa             <  1:8922/141‑1 (MQ=255)
                                                                                                                            ataataTTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAATTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCa             <  2:25522/141‑1 (MQ=255)
                                                                                                                             taataTTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCag            >  1:209027/1‑141 (MQ=255)
                                                                                                                               atatTAGCTGCTTTCCCGCCATCTAAAATCAAACCATGCACCCGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCagag          >  2:134300/1‑141 (MQ=255)
                                                                                                                                       tgctTACCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTGTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCAGAGGCATGAGc  <  1:313484/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
TTCATAATCACAACCAGTCTGTATCGCAGCTCCACGTGCGATTTGATTTACTTTTTCTGTTAATATATCCAATTCTTTACGCGTCATTGCTCTAGTATAAAAACGAGCATGTGTATAGTCTGGAATAATATTAGCTGCTTTCCCGCCATCTAAAATCACACCATGCACACGTTGATCTTTTTTAATATGTTGTCGTAGTTGTGCTACACCATTAAAATAACTAATCATAGCGTCTAATGCATTTAACGCTTCATCTGCATTTTCAGAGGCATGAGC  >  CP000730/2250965‑2251240

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