Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,678 T→C 100% H50R (CAC→CGC)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,6780TC100.0% 126.4 / NA 36H50R (CAC→CGC) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/15);  total (21/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAGTAC  >  CP006881/3724543‑3724811
                                                                                                                                       |                                                                                                                                     
gAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCa                                                                                                                                    <  2:123470/139‑1 (MQ=255)
       cgATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGTTTTGTTTGAACTTTTCATTATCGTCGCGCCCGACa                                                                                                                               >  3:176191/1‑137 (MQ=255)
               aTGTGATCGCTCAGCGGCGCCGCGCGCTCTCCTTCCGCATAGCTGATGATTTCCTCAGCTATTTCTATATGCTCCTCCGGCCGTGTTTGCAGGATTTGTTTGAAGTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                                                     >  7:231949/1‑139 (MQ=255)
               aTGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                                                     <  6:255465/139‑1 (MQ=255)
               aTGTGATCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                                                     <  4:310472/139‑1 (MQ=255)
                    aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                                                <  8:241739/139‑1 (MQ=255)
                    aTCGCTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                                                <  2:134787/139‑1 (MQ=255)
                    aTCGCTCAGCGGCGCCGCGAGCTCTCCGTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                                                >  7:294900/1‑139 (MQ=255)
                        cTCAGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATACTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGGTTTGTTTGAACTTTTCATTCTCGTCGCGCCCGACAAAGATGTTTTCCACTTTAt                                                                                                            >  5:217506/1‑139 (MQ=255)
                           aGCGGCGCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCa                                                                                                         <  2:37984/139‑1 (MQ=255)
                            gcggctCCGCGAGCTCTCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCAt                                                                                                        <  6:265616/139‑1 (MQ=255)
                                                      tAGCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCttttttt                                                                              >  6:232057/1‑139 (MQ=255)
                                                        gCTGATGATATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTg                                                                             >  3:6852/1‑138 (MQ=255)
                                                            atgatATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                                                <  6:272769/99‑1 (MQ=255)
                                                            atgatATCCTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACt                                                                                                                >  5:272769/1‑99 (MQ=255)
                                                                   ccTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCTTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATTAGATCATTTTTCTTTTTTTGAAAAGCGATTcc                                                                 >  7:60846/1‑139 (MQ=255)
                                                                   ccTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTcc                                                                 >  2:86213/1‑139 (MQ=255)
                                                                   ccTCAGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTcc                                                                 >  5:122128/1‑139 (MQ=255)
                                                                       aGCTATTTCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTATGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATAATTTTTCTTTTTTTGAAAAGAGATTcacggt                                                             >  1:137538/1‑134 (MQ=255)
                                                                            tttCTATATGCTCCTCCGGCAGTGTTTGCAGGATTTTTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCAATTTATTCATAGGGATGAGATAATTTTTCTTTTTTTGAAAAGCGATTCCAGGTCCCAt                                                        >  7:298792/1‑139 (MQ=255)
                                                                              tctatatGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGAGCTTTTCCACGTTATTCATAGGGATGAGATCATATTTCTTTTTTGGAAAAGCGATTCCCGGTCCCATCa                                                      >  1:308287/4‑139 (MQ=255)
                                                                               ctatatGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCa                                                      <  7:12885/136‑1 (MQ=255)
                                                                                         ctccGGCAGTGTTTTCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGAGATTCCCGGTCCCATCAACAttttttag                                           >  3:137764/1‑137 (MQ=255)
                                                                                         ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                                                     >  6:266990/1‑65 (MQ=255)
                                                                                         ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCtttt                                                                                                                     <  5:266990/65‑1 (MQ=255)
                                                                                         ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                                           >  3:133587/1‑139 (MQ=255)
                                                                                         ctccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                                           >  7:65058/1‑139 (MQ=255)
                                                                                         ctccGGCAGTGGTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTcc                                           <  4:263827/139‑1 (MQ=255)
                                                                                          tccGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCt                                          >  5:139300/1‑139 (MQ=255)
                                                                                            cGGCAGTGTTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGa                                        >  4:92246/1‑139 (MQ=255)
                                                                                                 gtgtTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGatcatc                                   <  4:6852/139‑1 (MQ=255)
                                                                                                 gtgtTTGCAGGATTTGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGatcatc                                   <  3:348836/139‑1 (MQ=255)
                                                                                                              ttGTTTGAACTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTaat                      >  6:183995/1‑138 (MQ=255)
                                                                                                                     aaCTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCACGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCAt               >  8:200576/1‑139 (MQ=255)
                                                                                                                       cTTTTCATTCTCGTCGCGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCAttg             <  6:123569/139‑1 (MQ=255)
                                                                                                                       cTTTTCATTCGCGTCACGCACGACAAAGATCGTTTCAACTTTATTCAACGCGATGAGATCATTTATCTTTTTTTTAAAAGCGATTCCCTGTCCCATCACGATTATTTCCTGATCATCCTCTTTTATTAAAGCCGCAttg             <  4:137764/139‑1 (MQ=255)
                                                                                                                                  cgtcgcGCACGACAAAGATCTTTTCCACTTTATTAATAGGGATGAGATCATTTTTCTTTTTTTGAAAACCGATTCCCGGACCCATCACCATTTTTTCCTGAACATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATaa  >  7:299902/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                     
GAAAGCGCGATATGAATGTGATCGCTCAGCGGCGCCGCGAGCTCCCCTTCCGCATAGCTGATGATATCCTCAGCTATTTCAATATGCTCCTCCGGCAGTGTTTGCAGGATTTGTTTAAACTTTTCATTCTCGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAGTAC  >  CP006881/3724543‑3724811

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