Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,078 C→T 100% S250N (AGT→AAT)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,0780CT100.0% 190.5 / NA 56S250N (AGT→AAT) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/28);  total (28/28)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCA  >  CP006881/3723943‑3724210
                                                                                                                                       |                                                                                                                                    
aaTCCTGATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCAtttt                                                                                                                                   >  2:181220/1‑139 (MQ=255)
      gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa                                                                                                                              >  4:258408/1‑138 (MQ=255)
      gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaa                                                                                                                              >  3:321618/1‑138 (MQ=255)
      gATCAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATACTGGCAAATGTAATTGATAATCATTTTTCaaa                                                                                                                              >  5:103574/1‑138 (MQ=255)
       atcaataaaAACACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaaaa                                                                                                                            <  8:390974/136‑1 (MQ=255)
       atcaataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCaaaaa                                                                                                                            <  4:268857/136‑1 (MQ=255)
          aataaaAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAActc                                                                                                                         >  6:245256/1‑139 (MQ=255)
              aaaaaaCTTCCAGAGATCAATCCATGAAAGTCAAAGGGGTTATTCGGAGTGATCTTGAAGACGCTGGACTTGCAGCGTGATTTAGCCGGCCTCGGTTTCCTGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGc                                                                                                                      <  7:350242/138‑1 (MQ=255)
                       ccAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgc                                                                                                            >  1:324082/1‑139 (MQ=255)
                        cAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgcg                                                                                                           <  4:88907/139‑1 (MQ=255)
                          gagaTCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGca                                                                                                         <  5:7679/139‑1 (MQ=255)
                                aCTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGa                                                                                                   <  8:220292/139‑1 (MQ=255)
                                  tCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGa                                                                                                   >  8:111863/1‑137 (MQ=255)
                                    ccTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAgg                                                                                               <  2:324082/139‑1 (MQ=255)
                                       gCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGcg                                                                                             >  3:285202/1‑137 (MQ=255)
                                            tCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgcgc                                                                                                          <  6:36504/120‑1 (MQ=255)
                                            tCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAgcgc                                                                                                          >  5:36504/1‑120 (MQ=255)
                                            tCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGaa                                                                                       <  3:258408/139‑1 (MQ=255)
                                            tCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGaa                                                                                       >  4:149227/1‑139 (MQ=255)
                                                     ttATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                     <  3:34391/100‑1 (MQ=255)
                                                     ttATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGcc                                                                                                                     >  4:34391/1‑100 (MQ=255)
                                                     ttATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAAGACtttttt                                                                              >  8:357628/1‑139 (MQ=255)
                                                          ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGTATACTTTTTTTGGa                                                                          <  4:366600/138‑1 (MQ=255)
                                                          ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                          <  5:233140/138‑1 (MQ=255)
                                                          ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                          <  8:32231/138‑1 (MQ=255)
                                                          ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                          <  1:364101/138‑1 (MQ=255)
                                                          ccGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGa                                                                          <  2:178490/138‑1 (MQ=255)
                                                             agagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATTTTCATTTTTCAAAAACTCCGCCAGCTCCTGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                       >  8:331869/1‑138 (MQ=255)
                                                             agagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGa                                                                       >  6:315609/1‑138 (MQ=255)
                                                              gagaTCTTGAAGCCGCTGGACATGCAGCGTGATATAGCCGGCCTCGGTTTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGCTCGGCGAATGAATACTTTTTTTGGATGaaa                                                                     <  4:231135/139‑1 (MQ=255)
                                                              gagaTCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGaaa                                                                     <  2:79273/139‑1 (MQ=255)
                                                                  tCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATa                                                                   >  5:108316/1‑137 (MQ=255)
                                                                      gAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGc                                                               <  7:374147/137‑1 (MQ=255)
                                                                            gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGca                                                                                                         <  8:203161/89‑1 (MQ=255)
                                                                            gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGca                                                                                                         >  7:203161/1‑89 (MQ=255)
                                                                            gCTGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATctcctc                                                       >  1:245374/1‑139 (MQ=255)
                                                                              tGGACATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGt                                                     >  1:116008/1‑139 (MQ=255)
                                                                                 aCATGCAGCGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCa                                                  <  1:292834/139‑1 (MQ=255)
                                                                                         cGTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATgg                                          <  8:108979/139‑1 (MQ=255)
                                                                                          gTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtctc                                                          <  4:77852/122‑1 (MQ=255)
                                                                                          gTGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCAtctc                                                          >  3:77852/1‑122 (MQ=255)
                                                                                           tGATATAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAg                                        <  5:245256/139‑1 (MQ=255)
                                                                                             atatAGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGTATGAAATAAAGCATCTCCTAGTCCCTACGATGGcgcg                                      >  6:264699/1‑139 (MQ=255)
                                                                                                taGCCGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCt                                    <  5:389797/138‑1 (MQ=255)
                                                                                                   ccGGCCTCGGATTCCGGCAAATGTTATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATAGTTTTTTTGGATGAAATAAAGCAGCTCATCGTCCATACTATGTCGCGCTTCAt                                >  2:254777/1‑139 (MQ=255)
                                                                                                   ccGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAt                                >  7:115540/1‑139 (MQ=255)
                                                                                                   ccGGCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAt                                <  3:235552/139‑1 (MQ=255)
                                                                                                      gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCa                                                  >  6:199210/1‑118 (MQ=255)
                                                                                                      gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCa                                                  <  5:199210/118‑1 (MQ=255)
                                                                                                      gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAttt                              >  2:211421/1‑138 (MQ=255)
                                                                                                      gCCTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTg                             >  4:273325/1‑139 (MQ=255)
                                                                                                      gCCTCGGATTCCGGCAAATGTAATTGATAGTCATTTTTCAAAAACTCCTCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTGTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCAttt                              >  5:187186/1‑138 (MQ=255)
                                                                                                       ccTCGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGa                            >  6:374278/1‑139 (MQ=255)
                                                                                                          cGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCt                                    <  2:240967/128‑1 (MQ=255)
                                                                                                          cGGATTCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCt                                    >  1:240967/1‑128 (MQ=255)
                                                                                                              ttCCGGCAAATGTAATTGATATTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAAc                     <  6:17697/138‑1 (MQ=255)
                                                                                                                                  atTCATTTTTCAAAAACTCCGCCAGCTCCAGCGCACATTGATAGGCGAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTGCATTTGAGTCCAACCGGCTGACCGCGTATCGCa  >  4:60299/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                    
AATCCCGATAAATAAAAAAACTTCCAGAGATCACTCCCTGCAAGTCAAAGCGGTTATTCCGAGAGATCTTGAAGACGCTGGACATGCAGCGTGATATAGCCGGCCTCGGACTCCGGCAAATGTAATTGATATTCACTTTTCAAAAACTCTGCCAGCTCAAGCGCACATTGATAGGCAAATGAATACTTTTTTTGGATGAAATAAAGCATCTCCTCGTCCATACGATGGAGCGCTTCATTTGATTCCAACCGACTGACTGCGTATCTCA  >  CP006881/3723943‑3724210

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