Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,028,542 G→C 22.2% L35V (CTC→GTC)  yvqF ← hypothetical membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,028,5420GC22.2% 80.3 / 22.2 36L35V (CTC→GTC) yvqFhypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base G (18/10);  new base C (6/2);  total (24/12)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

AGATGATGATTGTATACTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATAT  >  CP000730/2028404‑2028678
                                                                                                                                          |                                                                                                                                        
aGATGATGATTGTATACTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTa                                                                                                                                        >  2:123980/1‑141 (MQ=255)
 gatgatgatTGTATACTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTaa                                                                                                                                       >  2:278099/1‑141 (MQ=255)
      tgatTGTATACTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGc                                                                                                                                  >  1:285795/1‑141 (MQ=255)
         ttGTATACTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTAGGTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCaa                                                                                                                               >  2:175080/1‑141 (MQ=255)
                cTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTc                                                                                                                        <  1:101218/141‑1 (MQ=255)
                             aatGTGAATGCAATTAATAGCGCGCCTATACAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATgtag                                                                                                           <  2:251320/141‑1 (MQ=255)
                               tgtgAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGACCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATgtagta                                                                                                         <  2:123365/141‑1 (MQ=255)
                               tgtaAATACAATTAATAGCGCGACTATCCAAAACGCCAAAAGGACACATATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCAAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATgtagta                                                                                                         <  2:240670/141‑1 (MQ=255)
                                gtgAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTaa                                                                                                        >  1:28914/1‑141 (MQ=255)
                                   aaTGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTaaaaa                                                                                                     >  2:4455/1‑141 (MQ=255)
                                      gCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTg                                                                                                  <  2:134838/141‑1 (MQ=255)
                                         aTTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCa                                                                                               >  2:358115/1‑141 (MQ=255)
                                                 cgcgcCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCAt                                                                                       >  2:12286/1‑141 (MQ=255)
                                                     ccTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAAt                                                                                   >  2:184723/1‑141 (MQ=255)
                                                       tATCCAAAACGCCAAAAGGCAACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGACTAAAAAGCAAAAAAATTAAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGc                                                                                 <  2:208975/141‑1 (MQ=255)
                                                        aTCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCa                                                                                >  2:63097/1‑141 (MQ=255)
                                                                 gacAAAAAGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGACTAAAAAGACAATATTTTAAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATg                                                                       <  2:204392/139‑1 (MQ=255)
                                                                  ccAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGACTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGa                                                                      >  1:302845/1‑141 (MQ=255)
                                                                      aaGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGACTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCaa                                                                  >  2:153405/1‑141 (MQ=255)
                                                                                    tATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGata                                                    <  1:4455/141‑1 (MQ=255)
                                                                                                  tatCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTc                                      <  2:39036/141‑1 (MQ=255)
                                                                                                   atCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCa                                     <  1:301097/141‑1 (MQ=255)
                                                                                                     ccTACATAAACTCATACACATCCCCATAATAGAACGAGTAAAAAGCCAATTTTCTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTCATATATATTTGCGTGTCATg                                   >  2:303114/1‑141 (MQ=255)
                                                                                                       tACATAAACTAATACACATCCCAATAATAGACCGACTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGtt                                 >  2:44178/1‑141 (MQ=255)
                                                                                                        aCATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGttt                                >  2:20105/1‑141 (MQ=255)
                                                                                                               cTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCgccgcc                         >  1:227202/1‑141 (MQ=255)
                                                                                                                 aaTACACATCCCAATAATAGAACGACTAAAAAGCCAATTTTTTTTAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATCTATATTTCTGACTCATGTTACGCAGCCtt                       >  2:130269/1‑141 (MQ=255)
                                                                                                                 aaTACACATCCCAATAATAGAACGACTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCtt                       >  1:233613/1‑141 (MQ=255)
                                                                                                                           ccAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCcac             >  1:63221/1‑141 (MQ=255)
                                                                                                                            cAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATCTATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCcaca            <  2:319753/141‑1 (MQ=255)
                                                                                                                            cAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCcaca            >  2:370499/1‑141 (MQ=255)
                                                                                                                            cAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCcaca            >  1:178838/1‑141 (MQ=255)
                                                                                                                            cAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCcaca            >  1:178872/1‑141 (MQ=255)
                                                                                                                             aataatAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAg           <  1:92699/141‑1 (MQ=255)
                                                                                                                             aataatAGAACGACTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAg           >  1:89155/1‑141 (MQ=255)
                                                                                                                                      aCGAGTAAAAAGCCAATTTTTTTTTTTTTTTTTTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCTTCTGCGTTCatat  >  2:92876/1‑141 (MQ=255)
                                                                                                                                          |                                                                                                                                        
AGATGATGATTGTATACTTATTAGACAATAATGTGAATGCAATTAATAGCGCGCCTATCCAAAACGCCAAAAGGCCACGTATTTTATGAAAATAAAGATATCCTACATAAACTAATACACATCCCAATAATAGAACGAGTAAAAAGCCAATTTTTTCAAAAAATATGTAGTAAAAATTGGCAATAATCATTAATGCAGTAAAAATGATCAACATTTGCGTTGATATATATTTGTGTGTCATGTTTCGCCGCCTTTCTATGCCACAGCGTTCATAT  >  CP000730/2028404‑2028678

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