Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,539,250 C→G 100% intergenic (‑349/+37) aroC ← / ← USA300HOU_1405 chorismate synthase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,539,2500CG97.6% 47.0 / ‑6.3 17intergenic (‑349/+37)aroC/USA300HOU_1405chorismate synthase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (7/9);  minor base A (0/1);  total (7/10)
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.

CATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTCCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTGT  >  CP000730/1539118‑1539385
                                                                                                                                    |                                                                                                                                       
catcatATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTaaa                                                                                                                                 >  2:160323/1‑141 (MQ=255)
                  taaAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTaaaa                                                                                                               >  1:282044/1‑141 (MQ=255)
                     aaaTGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAatat                                                                                                            >  1:293640/1‑141 (MQ=255)
                                                                  taataTACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATaa                                                               <  1:59325/141‑1 (MQ=255)
                                                                   aataTACTGACCACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAAt                                                              <  1:294483/141‑1 (MQ=255)
                                                                             cTACATTTCTTAATACAAAAGAACGAACTTTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCtttt                                                    <  2:24803/141‑1 (MQ=255)
                                                                             cTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCtttt                                                    <  1:236001/141‑1 (MQ=255)
                                                                                  tttCTTAATACAAAAGAACAAACTGGTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGATAAGCGCTTTTGATATACTTTGACTTCAAAAATATAATCAATTCTTTTACTaa                                               >  1:180182/1‑141 (MQ=255)
                                                                                       tAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACAt                                          <  1:83522/141‑1 (MQ=255)
                                                                                          aaaaaaaGAACGAACTTTTAAATATTAATAAAGTTAATCAGAAAGTATAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATaaa                                       <  2:181081/138‑1 (MQ=255)
                                                                                            cAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATaaaaa                                     >  2:222833/1‑141 (MQ=255)
                                                                                                         tGTTAATTATTAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGagag                        <  2:24589/141‑1 (MQ=255)
                                                                                                              attattAATAAAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATgaacatcagt                   >  2:283866/1‑133 (MQ=255)
                                                                                                                    aataaAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTAAAAAAAAAAAAAAAATATAAATGAGAGTCAGTATaaaa             <  2:282044/141‑1 (MQ=255)
                                                                                                                      taaAGTTAATCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATTTAAATGAGAGTCAGTATAAAaat           >  1:197858/1‑141 (MQ=255)
                                                                                                                              aTCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTg   <  2:10871/141‑1 (MQ=255)
                                                                                                                               tCAGTGCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTGt  <  2:202515/141‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
CATCATATAATCATTAATTAAAAATGATTAGCCGTTTTAATATAATCGTATCGCAACTTTAGCAAATAATATACTGACTACATTTCTTAATACAAAAGAACGAACTGTTAATTATTAATAAAGTTAATCAGTCCGTTTAAAAAAATTATGATATTAAAATATACATATGTTAAGCGCTTTTGTTATACTTTGACTTCAAAAATATAATCAATTCTTTTACTAATACATAAAAAATATAAATGAGAGTCAGTATAAAAATAATTTTTGT  >  CP000730/1539118‑1539385

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