Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,389,706 C→A 20.8% M158I (ATG→ATT ynaI ← mechanosensitive channel protein, very small conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,389,7060CA20.8% 43.5 / 3.1 23M158I (ATG→ATTynaImechanosensitive channel protein, very small conductance
Reads supporting (aligned to +/- strand):  ref base C (9/9);  new base A (2/3);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

AACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTGGGTTTTCTACGCTGATCGACGAAAACAGCGAGTTCGGTACGTACAATGGACGATTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGAT  >  CP009273/1389453‑1389948
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
aaCCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTGGGTTTTCTACGCTGATCGACGAAAACAGCGAGTTCGGTACGTACAATGGACGATTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTAAGAATAACATTACCGGCCATACCGACAGCCAGaccacc                                                                                                                                                                                                                                >  2:172695/1‑274 (MQ=255)
                                     tgctcgattttATACGCGGATCGACGAAAACAGCGAGTTTGTGACGTACAATGGACTATTGTCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTGCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATg                                                                                                                                                                                           <  2:83067/267‑1 (MQ=255)
                                                    gCTGATCGACGAAAACAGCGAGTTCGGTACGTACAATGGACGATTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAATTTACTCAGCATATCTTTCCCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGACTCATGCCGAAATGATCGCCaa                                                                                                                                                                           >  1:273892/1‑274 (MQ=255)
                                                                                        tGGACGATTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAACCAGAAAAAGAAACACAATAACAGTAATATtga                                                                                                                                       >  2:90194/1‑275 (MQ=255)
                                                                                           aCGATTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACAAAAGGTCAGCAAGCCAGAAAGGCTCATGCCAAAAAGTTCGCCATAAAGCAGAACAAGAACCACAATAATGGTGAttttattt                                                                                                                                     >  2:121454/1‑270 (MQ=255)
                                                                                              aTTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATAc                                                                                                                                 <  1:172701/275‑1 (MQ=255)
                                                                                                         tCGTAATTTTGGTAATTCGCCAGCCAATTTCCTCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAAt                                                                                                                       <  1:90196/274‑1 (MQ=255)
                                                                                                                                         gCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAAt                                                                                      <  1:121460/275‑1 (MQ=255)
                                                                                                                                                                               ccAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAACTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCtttt                                                                                           <  1:274443/232‑1 (MQ=255)
                                                                                                                                                                                       cAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCAt                                                                                                                                                                                                                                                  >  2:466165/1‑73 (MQ=255)
                                                                                                                                                                                           aCTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCg                                                                                                                                                                                                                                              <  1:466178/73‑1 (MQ=255)
                                                                                                                                                                                               aaaGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCAATACCGACAGCCAGacca                                                                                                                                                                                                                                  >  2:401171/1‑81 (MQ=255)
                                                                                                                                                                                               aaaGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCAATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCATCAGAAATAGCTTCATCTTTATAATAAc                                >  2:241425/1‑275 (MQ=255)
                                                                                                                                                                                                   gACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCAATACCGACAGCCAGACCACCaa                                                                                                                                                                                                                              <  1:401184/81‑1 (MQ=255)
                                                                                                                                                                                                     cGGTCGAAATGGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGCTATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCACAAt                          >  1:484491/1‑275 (MQ=255)
                                                                                                                                                                                                         aaaaaGAGAACATACTCCCGGAAACGAATTTATTAAAAATATCTTAACCGGCATAAACGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAgg                      <  2:484478/273‑1 (MQ=255)
                                                                                                                                                                                                                        tCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAg                                                                                                                                                                                                                                          >  1:428947/1‑48 (MQ=255)
                                                                                                                                                                                                                        tCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAgg                                                                                                                                                                                                                  >  2:257854/1‑72 (MQ=255)
                                                                                                                                                                                                                            ggAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAg                                                                                                                                                                                                                                      <  2:428934/48‑1 (MQ=255)
                                                                                                                                                                                                                            ggAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAg                                                                                                                                                                                                              <  1:257862/72‑1 (MQ=255)
                                                                                                                                                                                                                             gAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATCGc                                                 >  2:188883/1‑228 (MQ=255)
                                                                                                                                                                                                                             gAAAAGAAGTTACTCAGAATATCTTTACCGGCAATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGAt  <  1:241432/275‑1 (MQ=255)
                                                                                                                                                                                                                                 agaagTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTaa                                             <  1:188889/228‑1 (MQ=255)
                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
AACCAATGGTCGTGGTAATGCGGCGGTTGGTCATTCGTCCTGGGTTTTCTACGCTGATCGACGAAAACAGCGAGTTCGGTACGTACAATGGACGATTATCAAAGGTCGTAATTTTGGTAATTCGCCAGCCAATTTCCGCTACTGTACCTTCGATATTTCTGTCCGGTGAACGGATCCAGTCGCCAATACTGAAAGGACGGTCGAAATAGAGCATAATCCCGGAAAAGAAGTTACTCAGAATATCTTTACCGGCCATACCGACAGCCAGACCACCAATACCACCAAAGGTCAGCAAGCCAGAAAGGCTCATGCCGAAATGTTCGCCATAAAGCAGAACAAGCACCACAATAATGGTGATTTTGATGATACGCGACATAATCCGCGCACTGGTGATATCGCGACCTTTTTTAATCTGCTGTTTTTCAAACTGATTAATCAGCAGAAATAGCTTAATCGTCAGAATAACCGCAATCAGGGACGTACAGATAAAATCGAT  >  CP009273/1389453‑1389948

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