Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,808 G→A 100% L7L (CTA→TTA)  U712_19200 ← SacPA operon antiterminator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,8080GA97.1% 143.5 / ‑5.8 45L7L (CTA→TTA) U712_19200SacPA operon antiterminator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/28);  minor base T (2/0);  total (17/28)
Fisher's exact test for biased strand distribution p-value = 1.37e-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.

CGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAGTACTTTATAGATTTTCAAGAGGCTATCTCTCCCGCCAATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTG  >  CP006881/3724673‑3724929
                                                                                                                                       |                                                                                                                         
cgtcgcGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCAGCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTTTTTATTAc                                                                                                                        >  5:225187/1‑139 (MQ=255)
 gtcgcGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATa                                                                                                                         <  3:106863/137‑1 (MQ=255)
  tcgcGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTATTAAtt                                                                                                                      >  4:41236/1‑139 (MQ=255)
     cgCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGAGCATTGTTCTTTTTTTGAAAAGCGAAGCCCGGGCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTtat                                                                                                                   <  4:360847/138‑1 (MQ=255)
        aCGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGa                                                                                                                <  4:68278/139‑1 (MQ=255)
           aCAAAGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGAttt                                                                                                             >  7:27744/1‑139 (MQ=255)
               agagCTTTTCCACTTTATTCATAGTGATGCGAGCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCaa                                                                                                         <  8:270492/139‑1 (MQ=255)
               aGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCaa                                                                                                         <  6:392281/139‑1 (MQ=255)
               aGATCTTTTCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCaa                                                                                                         <  3:230782/139‑1 (MQ=255)
               aGATCTTTTCCACTTAATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCaa                                                                                                         <  4:209764/139‑1 (MQ=255)
                gATATTTTCCAGTTAATTCATCGCTCTCATATCTTTTTTATTTTTTTCAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAag                                                                                                        <  2:199160/139‑1 (MQ=255)
                    ttttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGAATCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTATAATTAAAGCCGCATTGTAGTTTAATACTTTATAGATTTTCAAGAggt                                                                                                    >  1:317493/1‑138 (MQ=255)
                      ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTtat                                                                                                                   >  3:170183/1‑122 (MQ=255)
                      ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTtat                                                                                                                   <  4:170183/122‑1 (MQ=255)
                      ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTa                                                                                                  >  3:56306/1‑139 (MQ=255)
                      ttCCACTTTATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTa                                                                                                  >  7:17684/1‑139 (MQ=255)
                              tATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTccc                                                                                          >  5:246789/1‑139 (MQ=255)
                              tATTCATAGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCATGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTacc                                                                                          >  8:148203/1‑139 (MQ=255)
                                ttCATAGGGATGATATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCGCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGc                                                                                        <  4:134040/139‑1 (MQ=255)
                                     aGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGa                                                                                                                >  7:177731/1‑110 (MQ=255)
                                     aGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGa                                                                                                                <  8:177731/110‑1 (MQ=255)
                                     aGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAAtt                                                                                   >  6:380110/1‑139 (MQ=255)
                                     aGGGATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGACAAtt                                                                                   >  6:114388/1‑139 (MQ=255)
                                        gATGAGATCATTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGc                                                                                >  6:228842/1‑139 (MQ=255)
                                          tGAGATCATTTTTCTTTTTTTGAAAATAGATTCACGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCTCATTGTTGTTTAATACTTTATATATTTTCAAGAGGTTATCTCTCACGCCAATTTGCtt                                                                              >  1:222304/1‑139 (MQ=255)
                                                 aTTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACt                                                                        <  8:17684/138‑1 (MQ=255)
                                                 aTTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACt                                                                        <  1:165116/138‑1 (MQ=255)
                                                 aTTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACt                                                                        <  8:52820/138‑1 (MQ=255)
                                                 aTTTTTCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACt                                                                        <  1:165150/138‑1 (MQ=255)
                                                   ttttCTTTTTTTGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGCATggg                                                                     >  7:193164/1‑139 (MQ=255)
                                                          tttttGAAAAGCGATTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGAtttt                                                                <  6:113630/137‑1 (MQ=255)
                                                           ttttGAAAAGCGATTCCCGGTCCCCTCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTAc                                                              <  6:323032/138‑1 (MQ=255)
                                                                       aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTATAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                 <  2:317493/139‑1 (MQ=255)
                                                                       aTTCCCGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCg                                                 <  8:281883/139‑1 (MQ=255)
                                                                          cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                               <  5:228842/135‑1 (MQ=255)
                                                                          cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                               <  4:56306/135‑1 (MQ=255)
                                                                          cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                               <  4:179842/135‑1 (MQ=255)
                                                                          cccggTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                               <  3:317955/135‑1 (MQ=255)
                                                                               taccATCACGATTTTTGACTGATCATCCTCTTTTATTAAAGACGCAATGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTaa                                         <  5:33696/137‑1 (MQ=255)
                                                                                          ttttttCATGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGc                               <  6:246789/138‑1 (MQ=255)
                                                                                               tCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGct                          <  3:41236/138‑1 (MQ=255)
                                                                                                               ttattaAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                               <  2:233077/101‑1 (MQ=255)
                                                                                                               ttattaAAGCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGcc                                               >  1:233077/1‑101 (MQ=255)
                                                                                                                 attaAAGCCGCATTGTTGGTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGAt        <  5:82780/138‑1 (MQ=255)
                                                                                                                       gCCGCATTGTTGTTTAATACTTTATAGATTTTCAAGAGGTTATCTCTCCCGCCAATTTGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTg  >  4:156360/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                         
CGTCGTGCACGACAAAGATCTTTTCCACTTTATTCATAGGGATGAGGTCATTTTTCTTTTTTTGAAAAGCGATTCCTGGTCCCATCACGATTTTTTCCTGATCATCCTCTTTTATTAAAGCCGCATTGTTGTTTAGTACTTTATAGATTTTCAAGAGGCTATCTCTCCCGCCAATTAGCTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTG  >  CP006881/3724673‑3724929

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