Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,440,085 G→T intergenic (‑45/‑891) cspA ← / → xpaC cold shock protein CspA/5‑bromo‑4‑chloroindolyl phosphate hydrolysis protein XpaC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,440,0850GT95.7% 66.0 / ‑6.6 24intergenic (‑45/‑891)cspA/xpaCcold shock protein CspA/5‑bromo‑4‑chloroindolyl phosphate hydrolysis protein XpaC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (10/12);  minor base C (1/0);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATCAATTATTGGTAG  >  CP000730/1439940‑1440228
                                                                                                                                                 |                                                                                                                                               
ttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCa                                                                                                                                               >  2:216611/1‑148 (MQ=255)
     tGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAAACTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTa                                                                                                                                          >  2:6862/1‑148 (MQ=255)
         gAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaa                                                                                                                                      >  1:89538/1‑148 (MQ=255)
            aaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaaaaa                                                                                                                                   <  2:62511/148‑1 (MQ=255)
                  aCGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATg                                                                                                                             <  1:175702/148‑1 (MQ=255)
                         cGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCCTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAAttt                                                                                                                      <  2:209095/148‑1 (MQ=255)
                           tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCa                                                                                                                    <  2:301789/148‑1 (MQ=255)
                            cATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCAc                                                                                                                   >  1:34671/1‑148 (MQ=255)
                            cATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGAACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCCCATATTACAAAAAATACATGCAAATTTCAc                                                                                                                   >  1:211701/1‑148 (MQ=255)
                                                aTAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCACGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat                                                                                               <  2:305585/148‑1 (MQ=255)
                                                aTAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAATAAAATTCATTCAATATGCTTATAAAAAAAAAAAAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAatat                                                                                               <  2:146575/148‑1 (MQ=255)
                                                                   cAGCGTTAAACCATTTAACTGTACCATATTTCATAATCTGAAACCTCCAAGACTAAAATTCATGCAATATGCTTATTATAATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTaa                                                                            <  2:113389/148‑1 (MQ=255)
                                                                      cGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACtt                                                                         >  2:10042/1‑148 (MQ=255)
                                                                        ttAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGt                                                                       >  2:357101/1‑148 (MQ=255)
                                                                               aTTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAg                                                                >  1:170754/1‑148 (MQ=255)
                                                                                        tACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGttat                                                       <  1:216611/148‑1 (MQ=255)
                                                                                             tGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCa                                                  >  1:42183/1‑148 (MQ=255)
                                                                                             tGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTGATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCa                                                  >  1:11469/1‑148 (MQ=255)
                                                                                                 tCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAAtt                                              <  1:121695/148‑1 (MQ=255)
                                                                                                      aTCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGTAGTTATTACCAAATTGcaat                                         <  1:179491/148‑1 (MQ=255)
                                                                                                      aTCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaat                                         <  1:370350/148‑1 (MQ=255)
                                                                                                           aaaCCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatacaat                                    >  2:147577/1‑148 (MQ=255)
                                                                                                                               cattcaATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTaa                <  2:15915/148‑1 (MQ=255)
                                                                                                                                             aTTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATCAATTATTGGTAg  <  1:313149/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATTCATTTAATCAATTATTGGTAG  >  CP000730/1439940‑1440228

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: