Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,810,588 T→C 10.2% S266G (AGC→GGC)  secY2 ← Sec family Type I general secretory pathway protein SecY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,810,5880TC10.2% 106.9 / 3.0 39S266G (AGC→GGC) secY2Sec family Type I general secretory pathway protein SecY
Reads supporting (aligned to +/- strand):  ref base T (12/23);  new base C (2/2);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CTTTAAGAAATCTTTGGATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCATATTTGTTGCGGAAACGTTC  >  CP000730/2810448‑2810727
                                                                                                                                            |                                                                                                                                           
ctTTAAGAAATCTTTGGATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCt                                                                                                                                             >  1:456859/1‑141 (MQ=255)
ctTTAAGAAATCTTTGGATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCt                                                                                                                                             >  1:41536/1‑141 (MQ=255)
ctTTAAGAAATCTTTGGATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCt                                                                                                                                             <  1:528957/141‑1 (MQ=255)
 tttAAGAAATCTTTGGATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCtt                                                                                                                                            >  1:231148/1‑141 (MQ=255)
               ggATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGa                                                                                                                              >  2:328961/1‑141 (MQ=255)
                gATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGat                                                                                                                             <  1:285261/141‑1 (MQ=255)
                          tAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATaaaataaaat                                                                                                                   >  1:76757/1‑141 (MQ=255)
                              aTTGATTAAAAATCTCGATAACAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAaatgaat                                                                                                               <  1:88806/141‑1 (MQ=255)
                               ttGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAaatgaatg                                                                                                              <  1:334510/141‑1 (MQ=255)
                               ttGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAaatgaatg                                                                                                              <  2:547991/141‑1 (MQ=255)
                                          tctcGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTaa                                                                                                   >  2:130673/1‑141 (MQ=255)
                                           ctcGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAg                                                                                                  <  1:63446/141‑1 (MQ=255)
                                           ctcGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAg                                                                                                  <  2:20214/141‑1 (MQ=255)
                                                tAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAAc                                                                                             >  1:499434/1‑141 (MQ=255)
                                                       tAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCt                                                                                      >  1:67206/1‑141 (MQ=255)
                                                       tAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCt                                                                                      <  1:350060/141‑1 (MQ=255)
                                                         aCCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGa                                                                                    >  2:367963/1‑141 (MQ=255)
                                                                aaCATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACtc                                                                             <  1:15186/141‑1 (MQ=255)
                                                                    tttGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACtcatca                                                                         <  2:127122/141‑1 (MQ=255)
                                                                          acACTAAATAGACTGAACTACCTACGTGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAg                                                                   <  1:243595/140‑1 (MQ=255)
                                                                           cACTACATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGt                                                                  <  1:265666/141‑1 (MQ=255)
                                                                                aaTATACTCAAATACCGACTGCACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATAc                                                             <  1:473298/141‑1 (MQ=255)
                                                                                      cTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCCTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATActgcctg                                                       >  2:283421/1‑141 (MQ=255)
                                                                                      cTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCCTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAAAAAAAGCTGAAATACTCATCATTAAAGTAATActgcctg                                                       >  2:150197/1‑141 (MQ=255)
                                                                                             aCCTACGGGACTATCAAATGTCAGCCTTGGCATCTCCTCGTATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGtt                                                <  1:156223/141‑1 (MQ=255)
                                                                                              ccTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTa                                               <  1:83409/141‑1 (MQ=255)
                                                                                              ccTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTCTATTCAATATAGATAAAATAAAATGAATGCCAGTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATTCTGCCTGCAGGGTTa                                               <  1:421852/141‑1 (MQ=255)
                                                                                                         aTCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGa                                    <  1:328961/141‑1 (MQ=255)
                                                                                                          tCAAACGTCAGCATAGGCAGCTTATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGat                                   <  1:215312/141‑1 (MQ=255)
                                                                                                            aaaTGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTGATACTGCCTGCAGGGTTAACTTTCCAAGataa                                 <  2:262753/141‑1 (MQ=255)
                                                                                                            aaaTGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGttaa                                 <  2:208840/141‑1 (MQ=255)
                                                                                                            aaaTGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGataa                                 <  1:84188/141‑1 (MQ=255)
                                                                                                              aTGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGataaat                               >  2:28936/1‑141 (MQ=255)
                                                                                                                   aGCATTGGCATGTCATCTGATATGCCTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGAt                          <  1:150197/141‑1 (MQ=255)
                                                                                                                   aGCATTGGCATGTCATCTGATATGCCTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTACCAAGATAAATAAGAt                          <  1:283421/141‑1 (MQ=255)
                                                                                                                     cATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGAttt                        >  1:172296/1‑141 (MQ=255)
                                                                                                                       ttGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAATATTTCa                      >  2:82683/1‑141 (MQ=255)
                                                                                                                         ggCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCata                    >  2:95873/1‑141 (MQ=255)
                                                                                                                                   cTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCATATTTGTTGCgg          >  1:413572/1‑141 (MQ=255)
                                                                                                                                       tatGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCATATTTGTTGCGGAAAc      >  2:384278/1‑141 (MQ=255)
                                                                                                                                         tGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCATATTTGTTGCGGAAACGt    <  2:456859/141‑1 (MQ=255)
                                                                                                                                         tGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTCTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCATATTTGTTGCGGAAACGt    <  2:41536/141‑1 (MQ=255)
                                                                                                                                           cTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAATTTTCCAAGATAAATAAGATTTCATATTTGTTGCGGAAACGTTc  <  1:548438/141‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
CTTTAAGAAATCTTTGGATTTTTGTTTAGTATTGATTAAAAATCTCGATAAAAAATAACCCAATAACATTTGAATCACTAAATATACTGAAATACCTACTGGACTATCAAATGTCAGCATTGGCATGTCATCTGATATGCTTTTATTAAACATAGATAAAATAAAATGAATGCCACTTTTTAAGAAAACAAAAGCTGAAATACTCATCATTAAAGTAATACTGCCTGCAGGGTTAACTTTCCAAGATAAATAAGATTTCATATTTGTTGCGGAAACGTTC  >  CP000730/2810448‑2810727

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: