Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 1,250,361 C→T 13.9% A44V (GCA→GTA)  ftsY → Sec family Type I general secretory pathway protein signal recognition particle protein FtsY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,250,3610CT13.9% 82.7 / 5.2 36A44V (GCA→GTA) ftsYSec family Type I general secretory pathway protein signal recognition particle protein FtsY
Reads supporting (aligned to +/- strand):  ref base C (18/13);  new base T (3/2);  total (21/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GGAGCAATAATGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATGGGTGCTAAATTTAAA  >  CP000730/1250222‑1250497
                                                                                                                                           |                                                                                                                                        
ggagCAATAATGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCa                                                                                                                                         <  2:201906/141‑1 (MQ=255)
ggagCAATAATGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCa                                                                                                                                         <  1:218103/141‑1 (MQ=255)
ggagCAATAATGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCa                                                                                                                                         <  1:392303/141‑1 (MQ=255)
 gagCAATAATGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCaa                                                                                                                                        >  2:287613/1‑141 (MQ=255)
          tGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAg                                                                                                                               >  2:243661/1‑141 (MQ=255)
          tGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAg                                                                                                                               <  1:344285/141‑1 (MQ=255)
             gCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTATGAAGGTTCAACGCAGGACGCAAATGATTTAGCAg                                                                                                                            >  2:382902/1‑141 (MQ=255)
               ttttttAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGaa                                                                                                                          <  2:313086/141‑1 (MQ=255)
                   ttAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATg                                                                                                                      <  2:163416/141‑1 (MQ=255)
                             aaaaGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGaa                                                                                                            >  2:321759/1‑141 (MQ=255)
                                 gATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAaag                                                                                                        <  2:30901/141‑1 (MQ=255)
                                    aaGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAAGAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGTAAATGATTTAGCAGAAAATGCTGAAGTGAAAaagaag                                                                                                     >  2:427034/1‑141 (MQ=255)
                                          gCAACAAATAAAGGAAATGAAGAAGTTAAATGCTTAACAGAAGAACAAGGTAAAGAAAAATTAGAAGATACACATTCTGACGGCTCCACGCAGGATGCAAATGATTTAGCAGAAACTGCTGAAGCGCAAAAGAAGCCACgg                                                                                               >  2:278979/1‑140 (MQ=255)
                                               aaaTAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGtt                                                                                          <  2:384052/141‑1 (MQ=255)
                                                aaTAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTg                                                                                         >  2:394908/1‑141 (MQ=255)
                                                  tAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACGAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGaaaa                                                                                                          <  1:4285/122‑1 (MQ=255)
                                                  tAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACGAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGaaaa                                                                                                          >  2:4285/1‑122 (MQ=255)
                                                              agaagTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGAttt                                                                           >  1:55953/1‑141 (MQ=255)
                                                                 agTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTtga                                                                        <  1:90504/141‑1 (MQ=255)
                                                                    tAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTtgatga                                                                     >  1:230534/1‑141 (MQ=255)
                                                                         ccTTAAGAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGTAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATg                                                                >  1:373106/1‑141 (MQ=255)
                                                                             aaCAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCtt                                                            <  1:156088/141‑1 (MQ=255)
                                                                                 gaagaaCAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGTAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAata                                                        <  1:156480/141‑1 (MQ=255)
                                                                                     aaCAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCaa                                                    >  1:369141/1‑141 (MQ=255)
                                                                                      aCAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAAt                                                   >  1:57978/1‑141 (MQ=255)
                                                                                                 aCAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTg                                        >  2:37410/1‑141 (MQ=255)
                                                                                                    aaTTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTgaag                                     >  2:188404/1‑141 (MQ=255)
                                                                                                       tAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGaaa                                  >  1:263567/1‑141 (MQ=255)
                                                                                                           agaTACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGa                              >  1:427596/1‑141 (MQ=255)
                                                                                                                    ttCTGAAGGTTCAACGCAGGACGTAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCaaaa                     <  1:31111/141‑1 (MQ=255)
                                                                                                                        gAAGGTTCAACGCAGGACGTAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATg                 >  2:268782/1‑141 (MQ=255)
                                                                                                                          aGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATggg               >  2:308848/1‑141 (MQ=255)
                                                                                                                            gTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATGGGTg             >  1:424451/1‑141 (MQ=255)
                                                                                                                              tCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATGGGTGCt           <  1:329534/141‑1 (MQ=255)
                                                                                                                                  cGCAGCACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCGCAAAAAATGGGTGCTAAAt       <  1:278979/141‑1 (MQ=255)
                                                                                                                                       gACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATGGGTGCTAAATTTaaa  >  2:61433/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GGAGCAATAATGAGCTTTTTTAAACGCTTAAAAGATAAGTTTGCAACAAATAAAGAAAATGAAGAAGTTAAATCCTTAACAGAAGAACAAGGTCAAGACAAATTAGAAGATACACATTCTGAAGGTTCAACGCAGGACGCAAATGATTTAGCAGAAAATGCTGAAGTGAAAAAGAAGCCACGCAAGTTGAGTGAAGCGGATTTTGATGACGATGGCTTAATATCAATTGAAGATTTTGAAGAAATTGAAGCTCAAAAAATGGGTGCTAAATTTAAA  >  CP000730/1250222‑1250497

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