Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,659 A→G 100% F56F (TTT→TTC)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,6590AG99.6% 192.2 / ‑6.9 52F56F (TTT→TTC) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (27/24);  minor base C (0/1);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.81e-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.

CAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATT  >  CP006881/3724541‑3724788
                                                                                                                      |                                                                                                                                 
cggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgtcgcg                                                                                                               <  1:258783/137‑3 (MQ=255)
 ggAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCtcgtcgcg                                                                                                               >  1:104214/2‑136 (MQ=255)
  gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                                                                                             <  7:125128/139‑1 (MQ=255)
  gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                                                                                             <  3:39706/139‑1 (MQ=255)
  gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                                                                                             <  3:221099/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)
                                                      cATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACa                                                                                                        <  7:37597/92‑1 (MQ=255)
                                                      cATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACa                                                                                                        >  8:37597/1‑92 (MQ=255)
                                                      cATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCttttt                                                         >  3:181170/1‑139 (MQ=255)
                                                      cATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCttttt                                                         >  4:32604/1‑139 (MQ=255)
                                                      cATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACATTATTAATAGGGATGAGATCATTTTTCttttt                                                         >  1:251001/1‑139 (MQ=255)
                                                        tAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCttttttt                                                       >  6:174838/1‑139 (MQ=255)
                                                           cTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTGTGAACTTTTCATTCTCGTCGCGCACGCCAAAGATCTTTTCCACTTTATTCATTGGGATGAGATCATTTTTCTTTTTTTGa                                                     >  8:89578/1‑138 (MQ=255)
                                                              atgatATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                              >  5:72411/1‑94 (MQ=255)
                                                              atgatATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                              <  6:72411/94‑1 (MQ=255)
                                                              atgatATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                         <  7:95281/99‑1 (MQ=255)
                                                              atgatATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                         >  8:95281/1‑99 (MQ=255)
                                                                     ccTCAGCTATTTCTATATGCTCCTCCGGCAGTGGTGGCAGGATTTGTTCGAACTTTTCATTCTCGGCGCCCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTcc                                          <  6:293725/139‑1 (MQ=255)
                                                                         aGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGt                                      <  2:33674/139‑1 (MQ=255)
                                                                         aGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGt                                      <  2:38384/139‑1 (MQ=255)
                                                                         aGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGt                                      >  1:33674/1‑139 (MQ=255)
                                                                         aGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGt                                      >  1:38384/1‑139 (MQ=255)
                                                                           cTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTcc                                    >  5:253093/1‑139 (MQ=255)
                                                                              tttCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCAt                                 >  1:96944/1‑139 (MQ=255)
                                                                                 ctatatGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCa                               <  4:181170/136‑1 (MQ=255)
                                                                                         tcctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGAttttt                       >  3:217238/1‑138 (MQ=255)
                                                                                           ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCa                                                                                  <  6:95006/77‑1 (MQ=255)
                                                                                           ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCa                                                                                  >  5:95006/1‑77 (MQ=255)
                                                                                           ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGCTTCCCGCTCCCATCACTCTTTTTTcc                    >  1:246375/1‑139 (MQ=255)
                                                                                           ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                    >  6:148056/1‑139 (MQ=255)
                                                                                           ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                    <  6:138192/139‑1 (MQ=255)
                                                                                           ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                    >  8:176229/1‑139 (MQ=255)
                                                                                           cgccGGCAGTGTTTGCAGGATTTGTTTCAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGATATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                    <  5:24065/137‑1 (MQ=255)
                                                                                            tccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGaa                                                    <  6:92763/106‑1 (MQ=255)
                                                                                            tccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGaa                                                    >  5:92763/1‑106 (MQ=255)
                                                                                              cGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTGTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCCTCACGCTTTTTTCCt                   >  6:122989/1‑137 (MQ=255)
                                                                                                  aGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGatcat             <  4:217238/139‑1 (MQ=255)
                                                                                                     gtTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCct          >  7:211853/1‑139 (MQ=255)
                                                                                                     gtTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCct          >  3:118658/1‑139 (MQ=255)
                                                                                                        tgaaGGATTTGTTTGTACTTTTCATTCTCTTCGCGAACGACAAAGATCTTTTCCACTTTATTAATAGGGTTGTGATCTATTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCctct        <  2:246375/135‑1 (MQ=255)
                                                                                                              aTTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTttatt  >  6:265640/1‑138 (MQ=255)
                                                                                                              aTTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTttatt  >  4:292873/1‑138 (MQ=255)
                                                                                                                      |                                                                                                                                 
CAGAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATT  >  CP006881/3724541‑3724788

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: ‑