Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,264,192 C→T E14K (GAA→AAA)  qorA ← putative quinone oxidoreductase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,264,1920CT96.3% 90.3 / NA 27E14K (GAA→AAA) qorAputative quinone oxidoreductase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (16/10);  total (17/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 = 9.75e-01

ACACGATCGCCTGCCTTAATATGCTTTACACCACTGCCGACTTTACTCACGATGCCTGCTGCTTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCCTCCGCACTCGATCGAAGCGGAGCAGTCGGTGTTGGGCGGTTTA  >  NC_000913/4263978‑4264436
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
acacGATCGCCTGCCTTAATATGCTTTACACCACTGCCGACTTTACTCACGATGCCTGCTGCTTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTg                                                                                                                                                                                                                    <  1:254863/250‑1 (MQ=255)
         cctgcctTAATATGCTTTACACCACTGCCGACTTTACTCACGATGCCTGCTGCTTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGGTATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAAATCTACGGCTTTAAGTAATTTTGGGCCAACGTGCTTGTGAAATTCAATTCGGGTTGCCATGTGtcc                                                                                                                                                                                                          >  1:356022/1‑251 (MQ=255)
                     tgcattACACCACTGCCGACTTTACTCACGATGCCTGCTGCTTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGGTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATAGGGAGGAGTGAACTCTAACGCTTGTAAGAATTTTGGGGCCACCGGGTTGTGAAATTCAATTTGTGTTTTCATGTGTTCTCCCGAACGTAt                                                                                                                                                                                              >  1:204017/5‑251 (MQ=255)
                                           tACTCACGATGCCTGCTGCTTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATg                                                                                                                                                                         <  1:286689/250‑1 (MQ=255)
                                                         gctgctGCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCTCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTTGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGGATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGAcc                                                                                                                                                           <  1:179620/250‑1 (MQ=255)
                                                            gctTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCAt                                                                                                                                                        <  1:164484/250‑1 (MQ=255)
                                                                  gTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTTTTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGGGAAATCTACGGCTTTAAGGACTTTCTGGGCCCCCCGCTTGTGGAATTTAATTTGTGTTGCCCTGTGTCCTCCCGAACGTATCGTCCGGGGCTGGTTTATATTATAAAATTTTGAcccctttttttt                                                                                                                                                 >  2:196311/1‑242 (MQ=255)
                                                                                             ggcggGTAAAGGCCGCTGCGGATATATGTGTCGCTAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGGAGTACTTTCTGGCCACCGTGCTTGTGAAATTTCATTCGTGTTTCCCTTTGTCCTCCCGAACCTATAGGCAGGGGCTGCTTTATATGATCAAATTTCGACCCCTTTTTTTTCTCGGGTACTCCATTTACTATTGCAgg                                                                                                                      >  2:97803/1‑251 (MQ=255)
                                                                                                gggTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGaaa                                                                                                                   <  1:196311/251‑1 (MQ=255)
                                                                                                      aGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAAcc                                                                                                             <  1:271122/251‑1 (MQ=255)
                                                                                                       ggCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCCCCGTGCTTGTGAAATTCAATTCGTGTTTCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATCTGATAAAGTTTCGACCCCTTCTTTATCTCGGTAACTCCCTTCACTATGGCAGGGAATAAAccc                                                                                                            >  2:140929/1‑251 (MQ=255)
                                                                                                                 gATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGTTAACTCCATTCACTATGGCAGGAAATAAACCCTTCaacaaa                                                                                                   >  1:361411/1‑250 (MQ=255)
                                                                                                                                 aaaTTGATGCCGCTGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCcg                                                                                  >  1:87816/1‑251 (MQ=255)
                                                                                                                                   aTTGATGCCGATGGCTGTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTGTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCcgcg                                                                               <  2:58077/251‑1 (MQ=255)
                                                                                                                                           cGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGAt                                                                       >  1:191683/1‑251 (MQ=255)
                                                                                                                                              tGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTAATTTCGGGCCACCGTGCTTGTGAAAATCAATTCGTGTTGCCATGTGTCCTCCAGAACGAATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCAGAAGCGATcc                                                                     >  1:59194/1‑250 (MQ=255)
                                                                                                                                               ggCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCcac                                                                   >  2:51811/1‑251 (MQ=255)
                                                                                                                                               ggCTTTATTTTCGACCTGGATTTCATTCTCCACCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACCTTCGGACCAACATGCTTGTAAAATTCAATTCGTGATGCCCTATGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAATTTCCACCCATTCTTTATCTCGGTAACTCCATTCACTATCGACGGCAATAAACCCTTCAACAAACAGCCAGCTCAACCCCCCGAACGCGATCcac                                                                   >  1:26095/1‑251 (MQ=255)
                                                                                                                                                   ttATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGCACCCCGCGAACGCGATCCACAAGt                                                               >  1:117959/1‑251 (MQ=255)
                                                                                                                                                            aCCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCt                                                      >  1:23000/1‑251 (MQ=255)
                                                                                                                                                                            ccgccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCGGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCCTCCGCAc                                      <  1:150577/251‑1 (MQ=255)
                                                                                                                                                                              gccgGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCCTTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGTGATCCACAAGTTGCCGGGCTGAAAGTGCCCCCGCAccc                                    >  2:288025/1‑249 (MQ=255)
                                                                                                                                                                                                  ctctACGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCCTCCGCACTCGATCGAAGCGGAGCAGTCgg                >  1:366687/1‑251 (MQ=255)
                                                                                                                                                                                                      aCGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGa                                                                                                                                                                        >  1:97401/1‑96 (MQ=255)
                                                                                                                                                                                                      aCGGCTTGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGa                                                                                                                                                                        <  2:97401/96‑1 (MQ=255)
                                                                                                                                                                                                        ggcgtGAAGTACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGTAGGAAATAAACCCTGCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCCTCCGCACTCGATCGAAGCGGAGCAGTCGGTGTTgg          <  2:233794/247‑1 (MQ=255)
                                                                                                                                                                                                                  tACTTTCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGCGTCCTACAGAACGTATCGTCAGGGTCTGCTTCCTATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCCTCCGCACTCGATCGAAGCGGAGCAGTCGGCGGTGGGCGGTTTa  >  1:40997/1‑250 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                     
ACACGATCGCCTGCCTTAATATGCTTTACACCACTGCCGACTTTACTCACGATGCCTGCTGCTTCGGTGCCTAATCCGCTGGGTAGCGATGGCGGCGGGTAAAGGCCGCTGCGGATATATGTGTCGATAAAATTGATGCCGATGGCTTTATTTTCGACCTGGATTTCATTCTCCGCCGGATCGGCAGGAGTGAACTCTACGGCTTGAAGTACTTCCGGGCCACCGTGCTTGTGAAATTCAATTCGTGTTGCCATGTGTCCTCCAGAACGTATCGTCAGGGTCTGCTTCATATGATAAAGTTTCGACCCATTCTTTATCTCGGTAACTCCATTCACTATGGCAGGAAATAAACCCTTCAACAAACAGCAGGCTGAACCCCGCGAACGCGATCCACAAGTTGCCGGGCTGAAAGTGCCTCCGCACTCGATCGAAGCGGAGCAGTCGGTGTTGGGCGGTTTA  >  NC_000913/4263978‑4264436

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: