Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,172,775 A→G 8.3% V5A (GTT→GCT)  U712_16295 ← Uncharacterized protein yurT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,172,7750AG 8.3% 136.0 / 4.3 49V5A (GTT→GCT) U712_16295Uncharacterized protein yurT
Reads supporting (aligned to +/- strand):  ref base A (20/24);  new base G (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

ATCATCTGGAGAAACAAGCGTTATCCACCTGAATTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTATCCATCC  >  CP006881/3172648‑3172913
                                                                                                                               |                                                                                                                                          
atcatcTGGAGAAACAAGCGTTATCCACCTGAATTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGAtttt                                                                                                                                  >  7:236110/1‑138 (MQ=255)
     cTGGAGAAACAAGCGTTATCCACCTGAATTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTTTAAAACTCCAGTGCCATGTCTTGATCTTGTACATATATACTGGTAACAACGATTTTCATACt                                                                                                                            >  3:117/1‑139 (MQ=255)
             tcAAGCGGTATCCACCTGAATTTCCCCATGGGAACTTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATAGTGTTTTTcc                                                                                                                    <  4:102724/138‑1 (MQ=255)
             aCAAGCGTTATCCACCTGAATTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTcc                                                                                                                    <  8:153957/139‑1 (MQ=255)
                            cTGAATTTCCCCCTGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATatgaa                                                                                                     <  8:97275/139‑1 (MQ=255)
                            cTGAATTTCCCCATGGGAACGTCTTCCTTTTTTATAAAACCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATatgaa                                                                                                     >  6:293905/1‑139 (MQ=255)
                               aaTTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATatgaatg                                                                                                   <  4:282710/138‑1 (MQ=255)
                               aaTTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATatgaatg                                                                                                   <  6:38720/138‑1 (MQ=255)
                               aaTTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATatgaatg                                                                                                   <  3:184345/138‑1 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  1:48980/138‑1 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  1:67535/138‑1 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  4:117/138‑1 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  8:179130/138‑1 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  5:287547/138‑1 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCAt                                                                                             >  3:140111/1‑139 (MQ=255)
                                    ccccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGGATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  6:108124/138‑1 (MQ=255)
                                    acccATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  2:201817/137‑1 (MQ=255)
                                         tGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGat                                                                                                          <  4:288593/121‑1 (MQ=255)
                                         tGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGat                                                                                                          >  3:288593/1‑121 (MQ=255)
                                         tGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              >  2:94541/1‑133 (MQ=255)
                                         tGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCa                                                                                              <  1:94541/133‑1 (MQ=255)
                                           ggAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACt                                                                                       <  1:272260/138‑1 (MQ=255)
                                           ggAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTa                                                                                      <  8:236110/139‑1 (MQ=255)
                                                 tCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTaa                                                                                >  1:104919/1‑139 (MQ=255)
                                                                             cTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAaac                                                           >  3:199541/1‑132 (MQ=255)
                                                                             cTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAaac                                                           <  4:199541/132‑1 (MQ=255)
                                                                                tataAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCcaca                                                   <  4:140111/137‑1 (MQ=255)
                                                                                          cAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAGCAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCcac                                                    >  4:199046/1‑126 (MQ=255)
                                                                                          cAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAGCAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCcac                                                    <  3:203104/126‑1 (MQ=255)
                                                                                          cAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAGCAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCcac                                                    <  3:199046/126‑1 (MQ=255)
                                                                                          cAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAGCAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCcac                                                    >  4:203104/1‑126 (MQ=255)
                                                                                          cAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCagag                                       >  1:223654/1‑139 (MQ=255)
                                                                                             tGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGcc                                    <  5:293905/139‑1 (MQ=255)
                                                                                             tGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGcc                                    <  2:223654/139‑1 (MQ=255)
                                                                                             tGCCATGGCTTGCTCTTGTACAACTCTACTGGTACCACCGATTGTCCTACTGTTTCTCCCTCCTTCGCTCTGCCTGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACCCCCCCCACATTTAGCCGAGGcc                                    <  5:35465/139‑1 (MQ=255)
                                                                                               ccATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCt                                  >  8:264531/1‑139 (MQ=255)
                                                                                                    tCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCt                             >  3:272687/1‑139 (MQ=255)
                                                                                                         aTCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGAttt                        >  2:19190/1‑139 (MQ=255)
                                                                                                         aTCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGAttt                        >  5:74560/1‑139 (MQ=255)
                                                                                                              gTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTtata                    >  6:147840/1‑138 (MQ=255)
                                                                                                              gTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTtata                    >  1:180404/1‑138 (MQ=255)
                                                                                                                   aaTATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAatta              <  5:147840/139‑1 (MQ=255)
                                                                                                                        aCTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCCGTTCTGATTTTATAAAattatt            <  2:258678/136‑1 (MQ=255)
                                                                                                                        aCTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCCGTTCTGATTTTATAAAattatt            >  1:258678/1‑136 (MQ=255)
                                                                                                                         cTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAattatt            >  6:181293/1‑135 (MQ=255)
                                                                                                                         cTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAattatt            <  5:181293/135‑1 (MQ=255)
                                                                                                                            gTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTatcca     >  8:21509/1‑139 (MQ=255)
                                                                                                                            gTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTatcca     >  4:219457/1‑139 (MQ=255)
                                                                                                                            gTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTatcc      >  6:198287/1‑138 (MQ=255)
                                                                                                                              aacaacGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTatccatc   >  5:16287/1‑139 (MQ=255)
                                                                                                                               acaacGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTatccatcc  >  1:189859/1‑139 (MQ=255)
                                                                                                                               acaacGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTatcaatcc  >  3:169596/1‑139 (MQ=255)
                                                                                                                               |                                                                                                                                          
ATCATCTGGAGAAACAAGCGTTATCCACCTGAATTTCCCCATGGGAACGTCTTCCTTTTTTATAAACCCAAGTTTTTCTGTATAAAACTCCAGTGCCATGTCTTGATCTTGTACAAATATACTGGTAACAACGATTTTCATACTGTTTTTCCCTCCTTCGATATGAATGAAGCATGGTACTATGCTAACAATCTGGCTACAAGCAAAACAACCCACACATTTATCAGAGGCCCTGTTCTGATTTTATAAAATTATTCTATCCATCC  >  CP006881/3172648‑3172913

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