Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,587 C→T 100% G80G (GGG→GGA U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,5870CT99.1% 164.0 / ‑6.6 53G80G (GGG→GGAU712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/33);  minor base G (2/0);  total (20/33)
Fisher's exact test for biased strand distribution p-value = 1.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGAT  >  CP006881/3724454‑3724715
                                                                                                                                     |                                                                                                                                
ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                                                                                <  5:166190/136‑1 (MQ=255)
ttttATAGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGGCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTctcc                                                                                                                                <  1:201485/136‑1 (MQ=255)
   tataGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCg                                                                                                                           <  6:267585/138‑1 (MQ=255)
     ttgagCGCCTTTATCTCATGCAGCAATTTATTCTACACAAGCAGCCCATTTTGAATCCTTTCGATCACAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                         <  3:42095/136‑1 (MQ=255)
     taGCGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                         <  8:5153/138‑1 (MQ=255)
     taGAGCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                         <  5:63263/138‑1 (MQ=255)
        agCGCCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                     >  5:261075/1‑139 (MQ=255)
          cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                     >  7:284052/1‑137 (MQ=255)
          cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                     <  8:284052/137‑1 (MQ=255)
          cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCt                                                                                                                    >  4:127815/1‑138 (MQ=255)
          cgcCTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCt                                                                                                                    >  5:138192/1‑138 (MQ=255)
             ctttATCTAATGCAGAAATTTCTTCTGCACTAGCAGCGCATTTTGAATCATTTCTATCGCATATGAGAAATGGTCGGAAAGCGCGATATGATTGTGATCGCTCAGCGGCGCGGCGAGCTCTCCTTCCGCATAGCtgatg                                                                                                                <  8:40069/138‑1 (MQ=255)
                tATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTgatgat                                                                                                              <  3:157897/138‑1 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                         <  4:32006/120‑1 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCa                                                                                                                         >  3:32006/1‑120 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                      <  8:142696/139‑1 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                      >  7:142696/1‑139 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCGTCCGCATAGCTGATGATATCCTCAg                                                                                                      >  3:35446/1‑139 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATCTGATCGCTCAGCGGCTCCGCGAGCTCGCCTTCCGCATAGCTGATGAAATCcacac                                                                                                      >  3:130842/1‑138 (MQ=255)
                       tGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGGCGGAAAGCGCGATATGAATGTGAGCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAg                                                                                                      <  1:274032/139‑1 (MQ=255)
                          atcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                                                                                   <  5:306465/137‑1 (MQ=255)
                          agcaATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTa                                                                                                   <  4:131945/139‑1 (MQ=255)
                            caATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTAtt                                                                                                 >  2:140697/1‑139 (MQ=255)
                                   ttCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGGCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTatat                                                                                          <  7:82770/139‑1 (MQ=255)
                                        cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                     >  1:248580/1‑107 (MQ=255)
                                        cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGc                                                                                                                     <  2:248580/107‑1 (MQ=255)
                                        cacaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGctc                                                                                      <  3:127815/138‑1 (MQ=255)
                                         acaAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGCGCTCTCCTTCCGCATAGGTGATGATTTCCTCATCTATTTCTATATGctcct                                                                                    >  5:190974/1‑139 (MQ=255)
                                             gcagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCgg                                                                                >  8:290668/1‑139 (MQ=255)
                                              cagcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGc                                                                               <  3:54807/139‑1 (MQ=255)
                                                gcCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCCGCGGCGCCGCGCGGTCTCCTTCCGCATAGCTGATGATATCCTCAGCTAGTTCTATATGCTCCTCCGGCAg                                                                             >  3:90042/1‑139 (MQ=255)
                                                       ttGAATCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTg                                                                       <  7:290668/138‑1 (MQ=255)
                                                           aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGa                                                                  >  4:39706/1‑139 (MQ=255)
                                                           aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGa                                                                  >  3:225906/1‑139 (MQ=255)
                                                           aTCCTTTCGATCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGAGCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGa                                                                  <  3:71225/139‑1 (MQ=255)
                                                                    tgcGCAAAGGTGAAATGGTCGGAAAGCGCGATATGAATGTGTTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGtttga                                                         <  6:179737/137‑3 (MQ=255)
                                                                    aTCGCAAAGGACAAATGGTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGtttga                                                         <  6:261075/139‑3 (MQ=255)
                                                                                     gTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGc                                                                                        <  5:179441/91‑1 (MQ=255)
                                                                                     gTCGGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGc                                                                                        >  6:179441/1‑91 (MQ=255)
                                                                                       cggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgtcgcg                                      <  1:258783/137‑3 (MQ=255)
                                                                                        ggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgtcgcg                                      >  1:104214/2‑136 (MQ=255)
                                                                                         gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                    <  7:125128/139‑1 (MQ=255)
                                                                                         gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                    <  3:221099/139‑1 (MQ=255)
                                                                                         gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                    <  3:39706/139‑1 (MQ=255)
                                                                                          aaaGCGCGATATGAATGTGATCGCGCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCAc                                   <  5:188906/139‑1 (MQ=255)
                                                                                          aaaGCGCGATATGAATGCGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCAc                                   <  5:184905/139‑1 (MQ=255)
                                                                                                 gagatGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATCTGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAg                            <  8:82457/136‑1 (MQ=255)
                                                                                                           tgATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTcc                   <  1:140697/138‑1 (MQ=255)
                                                                                                              tCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                <  4:35446/138‑1 (MQ=255)
                                                                                                               cGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACtt               <  6:190974/138‑1 (MQ=255)
                                                                                                                 cTCAGCGGCGCCGCGCGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGGAGTGTTTGGAGGGTTTGTTTGTACTTTTCATTATCGTCGCGCACGACAAAGATGGTTTGCACttttt            >  1:153007/1‑137 (MQ=255)
                                                                                                                 cTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTAt            <  7:43069/139‑1 (MQ=255)
                                                                                                                 cTCAGCGGCGCCGCGAGCTCGCCTTCCGCATAGCGGCTGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTGGTTTGTAGTTTTGATTCTCGTCGCGCCCCACCAAGATCTTTTCCACgtttt            >  4:178105/1‑137 (MQ=255)
                                                                                                                     gcggcgCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTTTATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGATCTTTTCATTCTCGTCTCGCACGACAAAGATCTTTACCACTTTATTCAt        <  1:111819/139‑1 (MQ=255)
                                                                                                                           ccgcGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGAt  >  2:65234/1‑139 (MQ=255)
                                                                                                                                     |                                                                                                                                
TTTTATAGAGCGCTTTTATCTCATGCAGCAATTTATTCTGCACAAGCAGCCCATTTTGAATCCTTTCGATCGCAAAGGACAAATGGTCAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGAT  >  CP006881/3724454‑3724715

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: