Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 511,836 A→C 13.9% E55D (GAA→GAC U712_02450 → Protein sprT‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881511,8360AC13.9% 96.6 / 2.1 52E55D (GAA→GACU712_02450Protein sprT‑like protein
Reads supporting (aligned to +/- strand):  ref base A (27/17);  new base C (2/6);  total (29/23)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACAGAGGATATTTCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTGCAGCAGGTCAACGCGCCGAGGT  >  CP006881/511698‑511969
                                                                                                                                           |                                                                                                                                     
gACAGAGGATATTTCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGaa                                                                                                                                       >  2:228751/1‑139 (MQ=255)
 aacgagTATATTTATAAAAAATATTTAAAAAAGCCGTGCGCCCAGCAGGGGTTTTTAAATGACCGGAATAAAACCACGGGCGGCAGGTACCTGTTAACCTCCCATAACTTTTATTTCAACATAAAATATTTGATTGaa                                                                                                                                       <  4:235133/135‑1 (MQ=255)
  cAGAGGATTTTTCTCAAACATATTTTAAAACGCCGTTCCGCCATCTGGCGCTTTTTAATGACCGGATTAAAACCAGGGGCGGCAGGTACCGGTTAACCTCCAATAACATTGAATTTAACAGAAAATATTTCATTCAaca                                                                                                                                     <  8:14645/139‑1 (MQ=255)
     agGATATTTCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACgg                                                                                                                                  >  8:105426/1‑139 (MQ=255)
       gATATTTCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACgg                                                                                                                                  >  1:84906/1‑137 (MQ=255)
        ttatttattAAACATTTTTTAAAAAGCCGTTCCGCCATCCGGCGCTTTTTATTCACCGGCTTAAAACCACTGGCGGCAGGTACCTGTAAACCTCCAATAACATTGAATTCAACATAAAATATTTGATTGACCACGGCa                                                                                                                                <  2:61601/137‑1 (MQ=255)
           tttCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAgg                                                                                                                              >  6:195987/1‑137 (MQ=255)
             tCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGggag                                                                                                                           >  1:239204/1‑138 (MQ=255)
                 aaaCATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACt                                                                                                                      >  1:111018/1‑139 (MQ=255)
                 aaaCATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACt                                                                                                                      >  6:222085/1‑139 (MQ=255)
                   aCATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCCCTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGa                                                                                                                    >  2:183825/1‑139 (MQ=255)
                      tatTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTTACCTCCCATAACATAGAATTGAACAGCAAATATTTGATTGACCACGGCAGGGAGGAACTGATTg                                                                                                                 >  7:143266/1‑139 (MQ=255)
                       atTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTgg                                                                                                                >  2:290302/1‑139 (MQ=255)
                          taaaaacGCCGTTCCGCCTTCCGGCGTTTTTTATTGACCGGATAAAACCCTGTGGCGGGTGTTAGCTGTAAACCCCCAAAAACTTTCAATTCAACATAAAATTTTTTATTCACCACGGCAGGGAGGAACTGATTGGca                                                                                                              <  4:345028/136‑1 (MQ=255)
                             aaacGCCGTTCCCCAACCAGGCTTTTTTTATTGACCGGCAAAAAAACACTGGCGGCCGGTTCCTGTTAACCTCCCAAAACATTTAATTGAAAATAAAATTTTTTATTTACCACGGCAGGGAGGAACTGATTGGcatcat                                                                                                          <  4:316689/139‑1 (MQ=255)
                                   cGTACCGCCACCAGGCGAGTTTTAATGACCGGATTAAACCCCCGGGCGGCAGGCACCTGTAAACCTCCCATACCTTTGAATGGAACAGAAAATATTTGATTGACCACGGCAGGGAGGAACTGATTGGCATCATCAAGCa                                                                                                    <  3:319585/139‑1 (MQ=255)
                                     ttCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACg                                                                                                  >  3:205673/1‑139 (MQ=255)
                                     ttCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACg                                                                                                  >  3:210892/1‑139 (MQ=255)
                                          ccATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGaa                                                                                                                                       <  6:360684/97‑1 (MQ=255)
                                          ccATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGaa                                                                                                                                       >  5:360684/1‑97 (MQ=255)
                                          ccATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCAGAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGctc                                                                                             >  8:367045/1‑139 (MQ=255)
                                             tCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGc                                                                                          >  8:301889/1‑139 (MQ=255)
                                                                cttattAAAACCAGGGGCGGGAGTTACGTTTTAACATCCAAAACATTTTATTTGAAAATAAAATATTTTATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTg                                                                       <  4:210892/135‑1 (MQ=255)
                                                                cgtcTTAATCCCAGTGGCCGGAGGTAGCTTTTAACACCCAAAAATATTTATTTGAACAAAAAATATTTCATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTg                                                                       <  4:205673/136‑1 (MQ=255)
                                                                      aaaaCCACTGGCGGCAGGTACCTGTAAACCTCCAATAACATTTAATTGAACATAAAATATTTGATTTAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGaa                                                                 <  7:105426/139‑1 (MQ=255)
                                                                       aacccccTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTCAACATAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGaaa                                                                <  2:289078/133‑1 (MQ=255)
                                                                       aaaCCCCGGGCGGCAGGTACCTTTAAACCTCCAATAACATTTAATTTAAAATAAAATTTTTGATTCAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGaaa                                                                <  1:228751/139‑1 (MQ=255)
                                                                       aaaCCACTGGCGGCAGGTACGTGTTAACCTCCCATAACTTTGAATTGAACAGAAAATATTTGATTTAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGaaa                                                                <  4:86914/139‑1 (MQ=255)
                                                                             cTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAATCACGAGCTCTGCAATTATCATCTTCATCTTGAAGGAAAAGTGTa                                                          >  3:269068/1‑139 (MQ=255)
                                                                             cTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTa                                                          >  1:157147/1‑139 (MQ=255)
                                                                             cTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGCGAGGAAATCATTGGCATCATCAAGCACGAGATCTGCCATAATCATCTACATCTGGAAGGAAAAGGGTa                                                          >  2:299460/1‑139 (MQ=255)
                                                                             cTGGCGGCAGGTACCTGTTAACCTCCCAAAACATAGAATTGAACAGAAAATATTTGATTGAACACGGCAGTGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGt                                                           >  2:20560/1‑138 (MQ=255)
                                                                                gcggcAGGTACCTGTTAACCTCCCATAACTTTGAATTGAACATAAAATATTTTATTTAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACa                                                        <  7:299950/138‑1 (MQ=255)
                                                                                   gctgttACCTGTTTACAACCAAAAACTTTTATTTTAAAAGAAAAATTTTTAATAACCTCGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGca                                                    <  7:353111/134‑1 (MQ=255)
                                                                                   gcAGGAACCTGTAAACCACCCATAACATTCAATTGAAAACAAAATTTTTGATTCAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGca                                                    <  4:333501/139‑1 (MQ=255)
                                                                                     aGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGca                                                    >  5:58647/1‑137 (MQ=255)
                                                                                        tACCTTTTAACCAACCAAAAAATTAAATATAAAATAAAATATTTTAAATAACAAGGCAGGGAGGAACTGAGTGGCATCATCCAGCACGAGCTCAGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCacag                                                 <  1:346408/137‑1 (MQ=255)
                                                                                          ccTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCacagaga                                              >  3:279874/1‑138 (MQ=255)
                                                                                              ttAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCCTCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACACGCacagagaca                                            >  5:268190/1‑136 (MQ=255)
                                                                                                    ccccAAAAAATTTAAATCAACACAAAATATTTCATGCACCACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTa                                   <  2:262465/138‑1 (MQ=255)
                                                                                                     cccATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCAc                                                                                                   >  8:380628/1‑74 (MQ=255)
                                                                                                     cccATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCAc                                                                                                   <  7:380628/74‑1 (MQ=255)
                                                                                                     cccAAAACATTGAATTCAAAACAAAATATTTGATTCAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAg                                  <  1:183825/139‑1 (MQ=255)
                                                                                                       cATAACATTGAATTGAACATAAAATTTTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGaa                                <  8:233229/139‑1 (MQ=255)
                                                                                                       cATAACATTGAATTCAACATAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGaa                                <  4:133440/139‑1 (MQ=255)
                                                                                                                tgaattgaaCAGAAAATATTTGATTGACCACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTg                        >  6:61184/1‑139 (MQ=255)
                                                                                                                 gaattgaaCAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCAc                                                                                                   <  2:272299/63‑1 (MQ=255)
                                                                                                                 gaattgaaCAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCAc                                                                                                   >  1:272299/1‑63 (MQ=255)
                                                                                                                   attgaaCAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTgcag                     <  7:301889/139‑1 (MQ=255)
                                                                                                                       aaCAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGaa                                                                 >  4:139540/1‑91 (MQ=255)
                                                                                                                       aaCAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGaa                                                                 <  3:139540/91‑1 (MQ=255)
                                                                                                                            aaaaTATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTGCAGCAGGTCAAc            >  4:273854/1‑139 (MQ=255)
                                                                                                                            aaaaTATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGAACAGAGACAGAGATTTTAGAATGCTATTGCAGCAGGTCAAc            >  7:247474/1‑139 (MQ=255)
                                                                                                                                      attgaACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTgcag                     >  7:253949/1‑120 (MQ=255)
                                                                                                                                      attgaACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTgcag                     <  8:253949/120‑1 (MQ=255)
                                                                                                                                      attgaACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTGCAGCAGGTCAACGCGCCGAGGt  >  1:169663/1‑139 (MQ=255)
                                                                                                                                      attgaACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTGCAGCAGGTCAACGCGCCGAGGt  >  6:153402/1‑139 (MQ=255)
                                                                                                                                           |                                                                                                                                     
GACAGAGGATATTTCTGAAACATATTTTAAAAAGCCGTTCCGCCATCAGGCGCTTTTTAATGACCGGCTTAAAACCACTGGCGGCAGGTACCTGTTAACCTCCCATAACATTGAATTGAACAGAAAATATTTGATTGAACACGGCAGGGAGGAACTGATTGGCATCATCAAGCACGAGCTCTGCCATTATCATCTTCATCTTGAAGGAAAAGGGTACAAGCACAGAGACAGAGATTTTAGAATGCTATTGCAGCAGGTCAACGCGCCGAGGT  >  CP006881/511698‑511969

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