Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,700,747 G→A L84L (CTC→CTT acpS ← holo‑[acyl‑carrier‑protein] synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,700,7470GA89.5% 44.3 / ‑0.9 21L84L (CTC→CTTacpSholo‑[acyl‑carrier‑protein] synthase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base A (8/9);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

AATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCGAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTCGTTATCGCTTAATACGCGGCGTGCCAGGCGATCACCGGATCGGGCGATCACCGCTTCGATGCGAGCGATCTCCACAATATCCGTGCCTAAACCTAATATTGCC  >  NC_000913/2700510‑2700996
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
aaTCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTg                                                                                                                                                                                                                                               >  1:18900/1‑250 (MQ=255)
aaTCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTg                                                                                                                                                                                                                                               <  2:330020/250‑1 (MQ=255)
   cGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATa                                                                                                                                                                                                                                           <  1:414260/251‑1 (MQ=255)
       gttgttGAGGGTGCATGCTGCACAAAATTAACGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACATGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTATTGCCTCGCCCCATAGCCGTACCCGTGGTTTGCCAGACTCATCATTGAATCCTTc                                                                                                                                                                                                                                       >  2:113519/1‑251 (MQ=255)
             gAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTg                                                                                                                                                                                                                                 <  1:232339/251‑1 (MQ=255)
                                 aTTAAAGTTAAAAAGTAAAACCCCCGTTCCTTCCCAGTTCGGGGGTTTTACTTTTTAAAGAAAACGGTATTATTTTTAACTGTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCCAGCTCATCATTGAATACTTCAAATTTATTAAACGCCAGACCATTGc                                                                                                                                                                                                             <  2:97551/251‑1 (MQ=255)
                                          aaaaaGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTGTTACTTTTTAAAGAGAACGTTATTATATTTAACTGTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTGTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCGCATCATTGAAGACTTCAAATTGATTAAACGCCAGACCATGGCGGATCCCg                                                                                                                                                                                                     <  2:240841/250‑1 (MQ=255)
                                                                  ccAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCGAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGctt                                                                                                                                                                            >  2:173397/1‑251 (MQ=255)
                                                                                                        attTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCCCCGGc                                                                                                                                      >  1:262221/1‑251 (MQ=255)
                                                                                                                     aataatTACCGTGGCACAAGCATAGTTCCATTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCATCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCTGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTTCGttt                                                                                                                         <  2:428163/251‑1 (MQ=255)
                                                                                                                         atTACCATGACACAAGCATAGTGCCGCTCATCTGCCCGCGTTTCATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCCGACCACTGCGGATCCCGGTGCAAAACACTTTTACTGCTGCTTCTTTCACAGCAAACTGCTTCGCCAGAAAACGAAACCGCTGGTGGTGCGGTTTccc                                                                                                                     >  2:153045/1‑250 (MQ=255)
                                                                                                                               gTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACCTGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCAGTGCGGGTCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGGGTTGCGTTTTCCAGATAGc                                                                                                               >  2:153885/1‑251 (MQ=255)
                                                                                                                                   cacaAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGTATCCCGGTGCCAAACGCTTTTACTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGTTGGTGCGTTTTCCAGATAGCCCAt                                                                                                           >  1:113/1‑251 (MQ=255)
                                                                                                                                     caAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTc                                                                                                         <  1:384916/251‑1 (MQ=255)
                                                                                                                                             gTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTCGTTATCGc                                                                                                 >  1:20762/1‑251 (MQ=255)
                                                                                                                                              tGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTTCGTTTTCCAGATAGCCCATTc                                                                                                         >  2:195869/1‑242 (MQ=255)
                                                                                                                                              tGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTc                                                                                                         <  1:195869/242‑1 (MQ=255)
                                                                                                                                                     cATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTCGTTATCGCTTAATAc                                                                                          <  2:99460/250‑1 (MQ=255)
                                                                                                                                                                   aCATGTATATGTGCAACGCCCAGCTTTTCCGCCAGTCTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCAAGCTCATCATTGAATACGTCAAATTGATTAAACGCCAGTCCATTGCGGATCCCGGTGCCAAACGCTTTTGCGGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTCGTTAGCGCTTAATACGCGGCGTGCCAGGc                                                                            <  2:167214/250‑1 (MQ=255)
                                                                                                                                                                                     accAGCTTTTCCGCCAGTTTTAATGCCTCGCACCATAGCCGTAGCCGTGTTTTGCCAAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGACAAACGCTTTTGATGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTGCAAGATAGCCCATTAGTTATCGCTTAATACACGGGATGCCAGGAGTTCACCGGATCgagggt                                                          >  2:297018/2‑245 (MQ=255)
                                                                                                                                                                                                                                            cGAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTCCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTCGTTATCGCTTAATACGCGGCGTGCCAGGCGATCACCGGATCGGGCGATCCCCGCTTCGATGCGGGCGATCTCCACAATATCCGTGCCTAAACCTAATATTGcc  >  2:155143/1‑251 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                          
AATCGCTGTTGTTGAGGGTGCATGCTGCACAAAATTAAAGTTAAAAAGTAAAACCCCCGTTCCTTACCAGTTCGGGGGTTTTACTTTTTAAAGAGAACGGTATTATTTTTAACTTTCAATAATTACCGTGGCACAAGCATAGTGCCGCTCATCTGCCAGCGTTACATGCATATTTGCAACGCCCAGCTTTTCCGCCAGTTTTAATGCCTCGCCCCATAGCCGTAGCCGTGGTTTGCCGAGCTCATCATTGAATACTTCAAATTGATTAAACGCCAGACCATTGCGGATCCCGGTGCCAAACGCTTTTGCTGCGGCTTCTTTCACAGCAAAACGCTTCGCCAGAAAACGCACCGGCTGGTGGTGCGTTTTCCAGATAGCCCATTCGTTATCGCTTAATACGCGGCGTGCCAGGCGATCACCGGATCGGGCGATCACCGCTTCGATGCGAGCGATCTCCACAATATCCGTGCCTAAACCTAATATTGCC  >  NC_000913/2700510‑2700996

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