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.95.5% 90.9 / ‑1.5 22coding (34/201 nt)cspAcold shock protein CspA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (12/9);  minor base T (1/0);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  CP000730/1439872‑1440140
                                                                                                                                       |                                                                                                                                     
cgcgGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGt                                                                                                                                  >  1:61654/1‑140 (MQ=255)
                tactTCAAACTCAACAGCTTTACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCGTTTT‑AGCGTTAAACCATTTAACTga                                                                                                                  >  1:215706/1‑139 (MQ=255)
                 actTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTa                                                                                                                 <  1:71080/140‑1 (MQ=255)
                         cTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTGGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGttt                                                                                                         >  2:178019/1‑140 (MQ=255)
                           cAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCa                                                                                                       <  1:192195/140‑1 (MQ=255)
                                                aTGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACt                                                                                  <  1:8868/140‑1 (MQ=255)
                                                            cATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTCACTGTACCTTTTTTCATAATCTGAAACCTACAATACTAAAAttcattca                                                                      >  2:313774/1‑140 (MQ=255)
                                                                 tgtttAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                                                                 <  2:100563/137‑1 (MQ=255)
                                                                 tGGTTAATTGCTGAAAAATGTACGAATACGTCCTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                                                                 >  2:59144/1‑140 (MQ=255)
                                                                 tGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATg                                                                 <  2:14892/140‑1 (MQ=255)
                                                                      aaTTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAGAGCCGAATCCTCTTT‑AGCGTTAAACCATTTAGCTGTACCTTTTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTAt                                                            <  2:2470/140‑1 (MQ=255)
                                                                             gAAAAATGTACGAATACGTCTTTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAGTTCATTCAATATGCTTATTCGCata                                                     >  1:215364/1‑140 (MQ=255)
                                                                             gAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCata                                                     >  1:341/1‑140 (MQ=255)
                                                                                  aTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCACTCAATATGCCTATTCGCATATTACa                                                <  1:178019/140‑1 (MQ=255)
                                                                                   tGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaa                                               >  1:280825/1‑140 (MQ=255)
                                                                                   tGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCAAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGGCTAAAATTCATTCAATATGCTTATTGGCATATTACaa                                               >  1:143341/1‑140 (MQ=255)
                                                                                     tACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaaaa                                             >  2:181169/1‑140 (MQ=255)
                                                                                            aCGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACAt                                      >  2:207367/1‑140 (MQ=255)
                                                                                                           tCAACTTCGATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAGAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCaa                       >  2:176857/1‑140 (MQ=255)
                                                                                                                   gATAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCtt               <  2:306401/140‑1 (MQ=255)
                                                                                                                     tAAAGCCGAATCCTTTTT‑AGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAATCCTCCAAGACTAAAATTCATTTAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTg             <  1:203642/140‑1 (MQ=255)
                                                                                                                                 ccttttTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat  >  2:354772/3‑140 (MQ=255)
                                                                                                                                       |                                                                                                                                     
CGCGGTCGCCTTCAACTACTTCAAACTCAACAGCTTGACCTTCTTCTAATGATTTGTAACCATCTTGGTTAATTGCTGAAAAATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAAT  >  CP000730/1439872‑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: ‑