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,0850GT84.6% 83.6 / 10.6 41intergenic (‑45/‑891)cspA/xpaCcold shock protein CspA/5‑bromo‑4‑chloroindolyl phosphate hydrolysis protein XpaC
Reads supporting (aligned to +/- strand):  ref base G (4/2);  new base T (14/19);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATT  >  CP000730/1439953‑1440207
                                                                                                                                    |                                                                                                                          
aaTGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAAGCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTAc                                                                                                                    <  1:319819/141‑1 (MQ=255)
   gTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaa                                                                                                                 <  1:14568/141‑1 (MQ=255)
    tACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACaaaaa                                                                                                               >  2:156160/1‑141 (MQ=255)
     cgtattcgtacATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACaaaaa                                                                                                               >  1:305124/11‑141 (MQ=255)
      cGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAAt                                                                                                             >  1:174373/1‑141 (MQ=255)
            cGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGc                                                                                                       <  2:315088/141‑1 (MQ=255)
                      tcCTTCAACTTCGATAAAGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACt                                                                                             >  2:196186/1‑141 (MQ=255)
                        cTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTc                                                                                            >  2:368768/1‑141 (MQ=255)
                        cTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTc                                                                                            >  2:368734/1‑141 (MQ=255)
                                   aTAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCttt                                                                                 <  1:222306/141‑1 (MQ=255)
                                      aaGCCGAATCCTTTTTAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCAAATAAAAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCaa                                                                             <  2:136393/141‑1 (MQ=255)
                                       aGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCaa                                                                             >  2:242578/1‑141 (MQ=255)
                                        gCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAat                                                                            >  1:251963/1‑141 (MQ=255)
                                            aaTCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATgg                                                                        <  2:305124/141‑1 (MQ=255)
                                             aTCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGa                                                                       >  1:184375/1‑141 (MQ=255)
                                               cgtttttCAGCGGTAAACCATTTAACTGTAGCTTGTTTCATTATCTTAATCCTGCATGACGAAGTTTCATTCAATATGCTTATTCTCGTATTACAAAAAATACATGCAATTTGCGCTCAGATATTCTTTGCAATATGGaat                                                                     <  2:384963/139‑1 (MQ=255)
                                               ccTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaat                                                                     <  2:384680/141‑1 (MQ=255)
                                                cTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCCCATATTACAAAAAATACATGCAAATTTCACTCACATACTCTTTACAATATGGaata                                                                    >  2:117841/1‑141 (MQ=255)
                                                cTTTTNCAGCGTTAANCCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAATACTAAAATTCATTCAATATGCTTNTTCTNATATTACNAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaata                                                                    <  2:340432/141‑1 (MQ=255)
                                                  ttttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaataat                                                                  <  1:959/141‑1 (MQ=255)
                                                  ttttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGaataat                                                                  <  2:410039/141‑1 (MQ=255)
                                                   tttCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAatc                                                                 >  2:386550/1‑141 (MQ=255)
                                                                   tttAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGta                                                 >  2:68999/1‑141 (MQ=255)
                                                                          gTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGt                                          >  2:364573/1‑141 (MQ=255)
                                                                                 gTTTCATAATCTGAACCCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGtta                                   >  1:416366/1‑141 (MQ=255)
                                                                                    tCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGttatta                                <  2:305659/141‑1 (MQ=255)
                                                                                     cATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTAc                               <  1:132738/141‑1 (MQ=255)
                                                                                      aTAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTa                                                        >  2:152860/1‑115 (MQ=255)
                                                                                      aTAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTa                                                        <  1:152860/115‑1 (MQ=255)
                                                                                       tAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCa                             >  1:276590/1‑141 (MQ=255)
                                                                                           cTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATAGTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAAtt                         <  1:166978/141‑1 (MQ=255)
                                                                                             gAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGc                       <  1:103143/141‑1 (MQ=255)
                                                                                             gAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGc                       <  1:191469/141‑1 (MQ=255)
                                                                                             gAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGc                       <  1:195726/141‑1 (MQ=255)
                                                                                              aaaCCTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGca                      >  2:264183/1‑141 (MQ=255)
                                                                                              aaaCCTCCAAGACTAAAATTCATTCAATATGCTTATTCCCAGATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGca                      >  1:412775/1‑141 (MQ=255)
                                                                                                  cTCCAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatac                  <  1:38585/141‑1 (MQ=255)
                                                                                                    ccAAGACTAAAATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATACTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGcaatacaa                >  1:148817/1‑141 (MQ=255)
                                                                                                        gACTAATATTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACaataat            <  2:228599/141‑1 (MQ=255)
                                                                                                            aaaaTTCATTCAATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATaaaa        <  2:206731/141‑1 (MQ=255)
                                                                                                                  cattcaATATGCTTATTCTCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATAtt  >  1:54877/1‑141 (MQ=255)
                                                                                                                                    |                                                                                                                          
AATGTACGAATACGTCATTTTCTCCTTCAACTTCGATAAAGCCGAATCCTTTTTCAGCGTTAAACCATTTAACTGTACCTTGTTTCATAATCTGAAACCTCCAAGACTAAAATTCATTCAATATGCTTATTCGCATATTACAAAAAATACATGCAAATTTCACTCAAATATTCTTTGCAATATGGAATAATCATCTGCTTAACTTGTATTATAGTGCAGTTATTACCAAATTGCAATACAATAATAAAAAATATT  >  CP000730/1439953‑1440207

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: