Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,536,804 T→C 13.7% F93F (TTT→TTC yhgA → transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,536,8040TC13.7% 65.5 / 3.1 29F93F (TTT→TTCyhgAtransposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base T (16/9);  new base C (2/2);  total (18/11)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TCAAGCTGGAGTCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCG  >  CP009273/3536659‑3536946
                                                                                                                                                 |                                                                                                                                              
tCAAGCTGGAGTCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTTCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAACGCACCCCAAACAAACTGATCGCATTTCGc                                                                                                                                             >  1:45468/1‑149 (MQ=255)
   aGCTGGAGTCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCatg                                                                                                                                          >  1:384870/1‑149 (MQ=255)
           tCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTCCGCa                                                                                                                                            <  1:200915/139‑1 (MQ=255)
           tCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTCCGCa                                                                                                                                            <  1:201208/139‑1 (MQ=255)
           tCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTCCGCa                                                                                                                                            >  2:201208/1‑139 (MQ=255)
           tCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTCCGCa                                                                                                                                            >  2:200915/1‑139 (MQ=255)
                   cTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTg                                                                                                                          >  1:182377/1‑149 (MQ=255)
                          gatgaCGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAAt                                                                                                                   <  2:282710/149‑1 (MQ=255)
                                         gAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGa                                                                                                    <  1:53745/149‑1 (MQ=255)
                                              cTATTCCGAGGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTg                                                                                               <  2:292552/149‑1 (MQ=255)
                                                         tgCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACg                                                                                    <  1:115421/149‑1 (MQ=255)
                                                                   ggTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGAtgg                                                                          >  2:253276/1‑149 (MQ=255)
                                                                              agcAAGGCCCGGAATACGTCTAGTGTCTGTTTGAACATCAAAGCACTCCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGCGGCCGATGGTGGTGCCAtgg                                                               <  2:439987/147‑3 (MQ=255)
                                                                                   ggACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAGTCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTgttgtt                                                          >  1:103209/1‑149 (MQ=255)
                                                                                     aCCCGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGtttt                                                        <  1:263736/149‑1 (MQ=255)
                                                                                       cAGGATACATCTATTGTCTGATTGAACACCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTAAGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTAc                                                      >  1:350170/1‑149 (MQ=255)
                                                                                        aGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTGGGATGATGCGTTACGCTATTGCCGCAATGCAAAATCCCCGTGATGCAGGATACACACAGTTGCCGATGCTGGTGACATTGATGTTTTAcc                                                     >  2:342983/1‑149 (MQ=255)
                                                                                             aCATCTATTGTCTGATTGAACATCATAGCTACTTAAGCGCGCGTATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGt                                                <  1:437109/149‑1 (MQ=255)
                                                                                             aCATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGt                                                >  2:304169/1‑149 (MQ=255)
                                                                                                         tGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTACCGCAATGCAAAATCACCTTGATGCTGGATACAAAACATTGCTGATGCGGGTGCCATGGTTGTTTTACCACGGTATTGATAGcccc                                    >  2:281204/1‑149 (MQ=255)
                                                                                                          gATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTGGCCGATGGTGGTGCCATTGTTGTTTTACCCCGGTATTGAAAGCccct                                   >  1:90905/1‑149 (MQ=255)
                                                                                                               aaCATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACACAAAGTTGCCGATGGTGGTGCCATTGTTGTTTTTCCACGGTGTTGAAAGCccctatccc                              >  2:109329/1‑149 (MQ=255)
                                                                                                               aaCATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCccctatccc                              >  1:390797/1‑149 (MQ=255)
                                                                                                               aaCATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCccctatccc                              >  1:355845/1‑149 (MQ=255)
                                                                                                                  aTCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCccctatccctat                           >  1:105196/1‑149 (MQ=255)
                                                                                                                  aTCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTCCCGCTATCCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGTTGCCATTGTTGTTTTACCACGGTATTGAAGACcccgttccctat                           >  2:194469/1‑149 (MQ=255)
                                                                                                                      aaGCACCTCAAACAAACTGATCGCATATCGCATGATGCGATACGCTATTGCCGCAATGCACAATCACCTTGGTGCATGATACGAAACGTTGCCGAGGGTGGTGCGATGGATGTTTTACCACGGTATTGAAAGTCCCACTCCCTTTTCGc                       >  2:36290/1‑149 (MQ=255)
                                                                                                                       aGCACCTCAAACAACCTGAACGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCCCCTATCCCTATTCGCt                      <  2:292976/149‑1 (MQ=255)
                                                                                                                                           cGCAGTTCGCATGATGCGTTACGCTATTGCCGCTCAACAACAACACCGTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCCCCTATTCCTATTCGCTGTGTTGGCTGGATTGTTTCg  <  1:243703/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCAAGCTGGAGTCGAGCAGCTTTATTGATGACGATCTGCGTGAAAGCTATTCCGATGTGCTGTGGTCGGTGAAAACGGAACAAGGACCAGGATACATCTATTGTCTGATTGAACATCAAAGCACCTCAAACAAACTGATCGCATTTCGCATGATGCGTTACGCTATTGCCGCAATGCAAAATCACCTTGATGCTGGATACAAAACGTTGCCGATGGTGGTGCCATTGTTGTTTTACCACGGTATTGAAAGCCCCTATCCCTATTCGCTGTGTTGGCTGGATTGTTTCG  >  CP009273/3536659‑3536946

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