Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,440,007 Δ1 bp 100% coding (34/201 nt) cspA ← cold shock protein CspA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,440,0070C.97.8% 205.4 / ‑3.3 46coding (34/201 nt)cspAcold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (28/17);  total (29/17)
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
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  CP000730/1439875‑1440140
                                                                                                                                    |                                                                                                                                     
ggTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTaaa                                                                                                                              <  1:611391/141‑1 (MQ=255)
             tactTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTa                                                                                                                 >  1:403008/1‑141 (MQ=255)
                  cAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTg                                                                                                            >  1:91474/1‑141 (MQ=255)
                   aaaCTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGt                                                                                                           >  1:287199/1‑141 (MQ=255)
                       tCAACCGCCTGTCCTCCTTTCAATGATTTTCAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCTATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                                                                       <  1:399830/141‑1 (MQ=255)
                       tCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                                                                       <  1:436237/141‑1 (MQ=255)
                          aCAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATaa                                                                                                    >  1:104477/1‑141 (MQ=255)
                           cAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAAt                                                                                                   >  1:41231/1‑141 (MQ=255)
                                  aCCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTCAACCACTTAACTGTACCTTGTTTCATAATCTGAAAc                                                                                            <  1:405263/141‑1 (MQ=255)
                                          ctAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTTTCCTTCAACTTCGTTAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATTTGTATCCTCCAAGa                                                                                    >  2:611120/1‑141 (MQ=255)
                                           tAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGaa                                                                                   >  2:328716/1‑140 (MQ=255)
                                           tAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGAc                                                                                   <  1:756185/141‑1 (MQ=255)
                                            aaTGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACTTCCAAGACt                                                                                  >  1:644834/1‑141 (MQ=255)
                                             aTGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCAATAAAACGTATAACTGGACCTTGTTTCATAATCTGAAACCTCCGAGACTa                                                                                 >  1:525397/1‑141 (MQ=255)
                                                       aCCATCTTGGTTGATTGCTGAAAAATGTACGACTACGTCATTTTCTCCTTCAACTTCGATAAAACCGAATCCCTTTTAGCGTTAAACCATTCAACTGTCCCTTGTTTCATCATATGAAACCCACAAGACTAAAcatcattc                                                                       >  2:651651/1‑141 (MQ=255)
                                                       aCCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCAGAATCCTTTATAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAattcattc                                                                       >  2:83070/1‑141 (MQ=255)
                                                            cTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAatat                                                                  <  2:104477/141‑1 (MQ=255)
                                                              tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGCCTAAAATTCATTCAATATGc                                                                >  1:192466/1‑141 (MQ=255)
                                                                 ttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTa                                                             >  1:719108/1‑141 (MQ=255)
                                                                   aaTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAtt                                                           >  1:354463/1‑141 (MQ=255)
                                                                   aaTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAtt                                                           >  2:527266/1‑141 (MQ=255)
                                                                    aTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTc                                                          >  2:302421/1‑141 (MQ=255)
                                                                       gCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCa                                                       >  1:748123/1‑141 (MQ=255)
                                                                          gAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCatat                                                    >  1:216279/1‑141 (MQ=255)
                                                                             aaaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTAc                                                 >  2:145648/1‑141 (MQ=255)
                                                                              attGTACGACTACGTCAGTGTCTCATTCAACTTCGTTAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACa                                                <  1:420215/139‑1 (MQ=255)
                                                                               aTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa                                               <  2:663140/141‑1 (MQ=255)
                                                                                 gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGCATCCTTTTCAGCGTTAAACCATTTAACTGTACCTTGCTTCATAATCTTAAACCTCCAAGACTAAAATTCACTCAATATTCTTATTCGCATATCACataa                                             >  1:407798/1‑141 (MQ=255)
                                                                                    cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAAt                                          >  1:497923/1‑141 (MQ=255)
                                                                                         tcGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAATCCATTTAACTGTACCTTGTTTCATAATCTGAAACCTTCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACGAAAAATACATg                                     <  2:231209/140‑1 (MQ=255)
                                                                                         atctCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGCTAAACCATTTAACTGTACCTTGGTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCGTATTCGCATATTACAAAAAATACATg                                     <  2:78349/138‑1 (MQ=255)
                                                                                         aCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATg                                     <  2:265631/141‑1 (MQ=255)
                                                                                            tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCaa                                  <  2:231590/141‑1 (MQ=255)
                                                                                            tCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCaa                                  <  1:710969/141‑1 (MQ=255)
                                                                                               tgttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAAtt                               <  1:437307/139‑1 (MQ=255)
                                                                                                 ttCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTc                             >  2:585452/1‑141 (MQ=255)
                                                                                                    tcCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACt                          >  1:689868/1‑141 (MQ=255)
                                                                                                     ccTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACAATTTAACTGTACCTTGTTTCATAATCTGAAACCTCGAAGACTAAAATTCATTCAATATGCTTATTCGCATATTGAAAAAAATACATGCAAAtttacatt                         >  1:111482/1‑136 (MQ=255)
                                                                                                           aCTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAatat                   <  1:607515/141‑1 (MQ=255)
                                                                                                              tCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCt                <  1:741135/141‑1 (MQ=255)
                                                                                                                    aaGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAGTTTCACTCATATATTCTTTGCaa          >  1:57478/1‑141 (MQ=255)
                                                                                                                     aGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAat         >  2:629592/1‑141 (MQ=255)
                                                                                                                     aGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAat         >  2:338824/1‑141 (MQ=255)
                                                                                                                      gCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCCTAATCTGAAACGTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACATCAAATACATGCAAATTTCACTCAAATATTTTTTCCAatt        >  2:607248/1‑140 (MQ=255)
                                                                                                                          aaTCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGa    >  1:116818/1‑141 (MQ=255)
                                                                                                                            tCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat  <  2:341537/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
GGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  CP000730/1439875‑1440140

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: