Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,894,839 T→C L12P (CTC→CCC)  pabB → aminodeoxychorismate synthase subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,894,8390TC84.2% 38.4 / ‑2.4 19L12P (CTC→CCC) pabBaminodeoxychorismate synthase subunit 1
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (9/7);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 5.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCC  >  NC_000913/1894611‑1895084
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
gCTCTTCTGCCACCAGCGCGATGGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTACATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTTCCCTGGCGTCAGGACGCCGCTGaa                                                                                                                                                                                                                                 <  1:234316/251‑1 (MQ=255)
                  cgATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCATGGATGCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTa                                                                                                                                                                                                               <  2:37724/251‑1 (MQ=255)
                   gATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTa                                                                                                                                                                                                               >  1:218767/1‑250 (MQ=255)
                                                cctgcAGTTCCGGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCGGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTGCATGCTCCCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGTACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTacac                                                                                                                                                                                 <  2:63223/247‑1 (MQ=255)
                                                                   cTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCATGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTACCGCTGTGATTACTTTAACCTGGCGACAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCTGGCTATGCCGatcatc                                                                                                                                                              >  1:201058/1‑251 (MQ=255)
                                                                        catcgttttGCTGCTGTTCATGGTGGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATGGTTGGCGCTAATTATGTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGCTTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCCCCCTCTTAAGCCTCCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGtat                                                                                                                                                         <  2:218767/246‑1 (MQ=255)
                                                                                        ttCATGGGAGAGTGAAGGTAAACCTACAAACATTGTTAAATCCTGCTAAAGTGTTGGCGCTAATTATTTCATGCTACCCGACACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGCTTACTTTGCCCTGGCGTCAGGACGCCGCTGACTTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTCCAATCCGTGTCTACCTATCATCCTTATAGCCGCTTTGTTATtg                                                                                                                                         >  2:75023/1‑251 (MQ=255)
                                                                                                  agtgaAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgccg                                                                                                                                >  1:71933/1‑250 (MQ=255)
                                                                                                       aGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCgccgagccgg                                                                                                                          >  1:166820/1‑250 (MQ=255)
                                                                                                        ggTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCTATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGAt                                                                                                                         >  1:137472/1‑251 (MQ=255)
                                                                                                                                             gtggcgcaaattatttaagccacaccGGCACATTGCCAGTAGATGCAGGACTGATGAAGCCTTTACCCCCCGCTGTGATTACTTTCCCCTGGCGTCAGGACGCCTCTGAATTTTATTTCCCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATGGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCgttgtt                                                                                    <  2:166820/227‑1 (MQ=255)
                                                                                                                                                         attTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTACGCCACCTGCCGTGGGCGCTGCTTTTACACTCCGGCTATGCCGCTCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCAATTTACACACTTTCGGTAAAGAAACCGTTGTTAGttgaagcgc                                                                         >  2:24007/1‑242 (MQ=255)
                                                                                                                                                                     ccGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACg                                                            >  1:123218/1‑251 (MQ=255)
                                                                                                                                                                           caTAGCCAGTAGAGTCAGGACTGATGAAACCGTTATCTCCCGCTGTGATTGCTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACt                                                      <  1:107718/251‑1 (MQ=255)
                                                                                                                                                                            aTAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACtg                                                     >  1:88230/1‑251 (MQ=255)
                                                                                                                                                                                    aaaatctacgACTAATGAAAACGTTATCTCCCGCTTTTATTACTTTCCTCTGGCGCCGTGACGCCGCTGAATTTTATTTCCCCCACTTAACCCCCCTGCCGTGGGCGATTCTTTTACACTCCGGCTATGCCTATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCAGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACCCACAACGACCACTGATGACCCg                                             <  2:137472/242‑1 (MQ=255)
                                                                                                                                                                                                         cgtcaTCTCCCGCTTTGTTTACTTGCCCCTGGCGTCAGGACGCCGCGGAATTTTATTTCTCTCGCTTAAGCCACCTGCCGTGGGCGATGTTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGTACTTTAACCACTTTCTGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGTCCCGCTACAGGTGCTCCAGCAGGTg                        <  1:157431/247‑1 (MQ=255)
                                                                                                                                                                                                          gTTATCTCCCGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGc                       <  1:145932/251‑1 (MQ=255)
                                                                                                                                                                                                            tATCTCACGCTGTGATTACTTTACCCTGGCGTCAGGACGCCGCTGACTTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGCTCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACAGTTGTTAGTGAAAGCGACAAACGCAAACCGACCACTGATGACCCGCTACAGGTGCTACATCAggggatg                     >  2:261848/1‑246 (MQ=255)
                                                                                                                                                                                                                                ttaccCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGcc  <  1:74163/250‑1 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     
GCTCTTCTGCCACCAGCGCGATTGCCTGGCCGCTGCTCATCCCTTGTGCCATCAGTTCCTGGATCCGCTCGACAGCTTTTTGCTGCTGTTCATGGGTGAGTGAAGGTAAACCTGCAAACATTGTTAACTCCTGCTAAATTGTTGGCGCTAATTATTTCATGCTACCCGGCACATAGCCAGTAGAGTCAGGACTGATGAAGACGTTATCTCCCGCTGTGATTACTTTACTCTGGCGTCAGGACGCCGCTGAATTTTATTTCTCCCGCTTAAGCCACCTGCCGTGGGCGATGCTTTTACACTCCGGCTATGCCGATCATCCGTATAGCCGCTTTGATATTGTGGTCGCCGAGCCGATTTGCACTTTAACCACTTTCGGTAAAGAAACCGTTGTTAGTGAAAGCGAAAAACGCACAACGACCACTGATGACCCGCTACAGGTGCTCCAGCAGGTGCTGGATCGCGCAGACATTCGCC  >  NC_000913/1894611‑1895084

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