Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,423,338 G→C intergenic (+26/‑62) ydaW → / → rzoR putative uncharacterized protein YdaW/putative prophage outer membrane lipoprotein RzoR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,423,3380GC93.8% 18.0 / ‑6.1 16intergenic (+26/‑62)ydaW/rzoRputative uncharacterized protein YdaW/putative prophage outer membrane lipoprotein RzoR
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (6/9);  minor base A (0/1);  total (6/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

AGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTT  >  NC_000913/1423189‑1423562
                                                                                                                                                     |                                                                                                                                                                                                                                
gtcaccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGTTGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGc                                                                                                                             >  1:210515/12‑251 (MQ=11)
   accgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGt                                                                                                                                >  1:202739/9‑245 (MQ=11)
   accgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGt                                                                                                                                <  2:202739/237‑1 (MQ=11)
           tATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgt                                                                                                                   <  1:202587/250‑1 (MQ=11)
                ccGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcg                                                                                                                 <  2:111449/247‑1 (MQ=11)
                                 atcatcgtcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGATGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGt                                                                                            >  1:383440/5‑251 (MQ=11)
                                                          ttAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGCAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGccacc                                                                   >  1:330743/1‑251 (MQ=11)
                                                                atgaCCGTGATAACCCCATTACCATCAAAGCCCAGCGCGACAAAAATGCCAGAAAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATTTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTATAGCGTAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCTCACCAAAGCAGTCTGGCAGTCAGTGCGTGAAGCCCCCACCGcc                                                             <  2:330743/248‑1 (MQ=11)
                                                                           aaCGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACGGAAGCTGGCGAACGCGGCAATTACGGACATGCAGATGCTTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTGGCCGCTGGTCGTCGTCGGTTGCACATCAAAGGGGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTgg                                                   <  2:132001/250‑1 (MQ=11)
                                                                              ccaattacccacaaagcccagtgctaCAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGAGGTTGCTGCGGTCGATGCAAAATACACGATGGAGTTAGCTTAAGCTAAAGCTGAAAATGATGCCCTGTGTGATGATTTTGCCGCTGTTCGTCTTCGGTTGCACATCAAAGCAGTCTTTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATaa                                               <  1:38158/226‑1 (MQ=11)
                                                                                             ccccAGCGCGCCAAAAATGCCAGAGAACTGAAGCGGGCTAACGCGGTAATTACTGAAATGCAGAGGCGTAAGCGTGATGTTGCTGCGCTCGATGAAAAATACACTAAGGAGTTAGCTGAGTCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGGCGCCGGTTGCCCATCAAAGCAGACGGTCAGCCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCg                                <  2:383440/248‑1 (MQ=11)
                                                                                                      taCAAAAATGCCAGAGAGCTGAAGCTGGCGAACGCGGCAATTCCTGACATGCAGATGCGTCAGCGTGAGGTTACTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAACCTGAAAATTATGCTCTGCGTGATGATGTTGGCGCTGATCGTCGTCGGTTGCACATCAAAGCAGTCTGCCAGTCAGTGCGGGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCgcagtga                          >  2:273839/2‑242 (MQ=11)
                                                                                                        cAAAAATGCCAGAGAACGTAAGCGGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGGTGCCGCTGGTCGTCGTCGGTTTCTCATCAAAGCAGTCGGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGaca                     <  1:398690/251‑1 (MQ=11)
                                                                                                        aaaaaagtccAGAGAACTGAAGCTGGCGAACGCTGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCGGCGTTCGATGCAAAATACACGAAGGATTTAGCTGATGCTAAAGCTGAAAATGAGGCTCTGCGTGATGATGTTGCCGCTGGTCGTGCTCGGTTGCACATCAAAGCAGTCTGTCAGTCATGGCGTGAAGCCCCAACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGaca                     <  2:95236/243‑1 (MQ=11)
                                                                                                             aTGCCAGAGAACTGAAGCGGGCGAACGCGGCAATTACTTTCATGCAGATGCTTCAGCGTGATGTTGCGGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCTCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCt                <  1:394977/251‑1 (MQ=11)
                                                                                                                           aaGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTAttt  >  2:350847/1‑251 (MQ=14)
                                                                                                                                                     |                                                                                                                                                                                                                                
AGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTT  >  NC_000913/1423189‑1423562

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