Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,725,043 C→A 100% A42A (GCG→GCT)  U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,725,0430CA99.6% 144.5 / ‑6.9 41A42A (GCG→GCT) U712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (6/34);  total (7/34)
Fisher's exact test for biased strand distribution p-value = 1.71e-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.

CTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAA  >  CP006881/3724911‑3725175
                                                                                                                                    |                                                                                                                                    
cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTATGCTTTTATTCCTCATTAACAGCCATCCGCCAATCATAAAACCAAACCAGTCCTGCCGCCTGCAAagccgcc                                                                                                                                >  8:28091/1‑138 (MQ=255)
     ttGGCGATTCCTTGAAGGCTTGCAAACATCTGCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCCCCATCCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                                                                                           <  4:108913/139‑1 (MQ=255)
       ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                                         <  8:263802/139‑1 (MQ=255)
       ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                                         <  3:214229/139‑1 (MQ=255)
        gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                                         <  8:98911/138‑1 (MQ=255)
        gCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                                                                                         <  6:301406/138‑1 (MQ=255)
                  gAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGACACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACg                                                                                                               <  8:91692/138‑1 (MQ=255)
                   agggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCt                                                                                                             <  7:213152/137‑1 (MQ=255)
                     ggCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGa                                                                                                           <  3:190765/139‑1 (MQ=255)
                      gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                          <  1:319108/139‑1 (MQ=255)
                      gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                          <  2:129576/139‑1 (MQ=255)
                      gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                          <  8:181386/139‑1 (MQ=255)
                      gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                          <  4:180928/139‑1 (MQ=255)
                      gCTTGCAAAAATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                                                                                          <  6:90332/139‑1 (MQ=255)
                           cAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGaa                                                                                                     <  5:84371/139‑1 (MQ=255)
                                     cttTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACCGCCGCCGTTGTTTTGTATTTTAGCAGAGCCGAATAGACCGccg                                                                                           >  5:173606/1‑137 (MQ=255)
                                          ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTgg                                                                                                                       <  4:345736/106‑1 (MQ=255)
                                          ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTgg                                                                                                                       >  3:345736/1‑106 (MQ=255)
                                          ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATACACCGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGc                                                                                       <  2:125917/138‑1 (MQ=255)
                                           taAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                                                     <  3:190632/139‑1 (MQ=255)
                                           taAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                                                     <  7:218522/139‑1 (MQ=255)
                                           aaaaGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAACCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCg                                                                                     <  1:138551/138‑1 (MQ=255)
                                                   ggcggcTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGttttt                                                                             <  8:247543/139‑1 (MQ=255)
                                                        cTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTc                                                                        <  5:175822/139‑1 (MQ=255)
                                                                  ctgcttttCTTCCTCATTCACCGCCAGCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCTCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCAc                                                               <  8:304392/134‑1 (MQ=255)
                                                                               tCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTGTTTGGTATTTACGCGGAGCCGCATAGACCGCGGCAGCCGCCGGTTTTTGGTCCGGCGTCACCCAGGTGATCaaa                                                  >  3:62099/1‑138 (MQ=255)
                                                                                        gCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATcc                                        <  3:17932/139‑1 (MQ=255)
                                                                                         ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                       <  1:121353/139‑1 (MQ=255)
                                                                                         ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCca                                       <  5:283664/139‑1 (MQ=255)
                                                                                         ccATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGGCACCCAGTTCATCAAAAAGACCATCca                                       <  7:28091/139‑1 (MQ=255)
                                                                                           aTCCGCCATTCATAACACGAATCGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGGTTTTCTGTCCGGCGGCACCCAGTTCATCAAAAAGACCATCcaca                                     <  8:351555/139‑1 (MQ=255)
                                                                                           aTCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCcaca                                     <  4:224426/139‑1 (MQ=255)
                                                                                                   ttCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGGTCATCACAAAGGCCATCCACACCCgga                               >  1:231886/2‑135 (MQ=255)
                                                                                                   ttCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTATGTAAGCGGCGCCGAATAGACCGCGGCCGCCCCAGTTTTTCTGTCCGGCGTCCCCCACTTCATCAAAAAGACCCACCACACCCCTaa                              >  5:132908/2‑138 (MQ=255)
                                                                                                       tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATAAAAAATACCATCCACACCCGTAAAGATa                         <  3:174108/139‑1 (MQ=255)
                                                                                                       tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGATGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATa                         <  4:126397/139‑1 (MQ=255)
                                                                                                       tAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCATTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATa                         <  7:162821/139‑1 (MQ=255)
                                                                                                            cGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGc                    <  6:173606/139‑1 (MQ=255)
                                                                                                               aaCGAGTCCTGCGGCCTGCACAGCCGCCGTTCTTTGGTATGTACCCTGAGCCGAATAGACACCGGCAGCCCCAGTTTTTATGTCAGTCGTCCCCCAGGTAATCAAAAATACCAACCACACCCGTAAAGATACAAGCaaa                 >  7:300107/1‑139 (MQ=255)
                                                                                                                              cTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGaaaaa  <  3:299626/139‑1 (MQ=255)
                                                                                                                              cTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGAATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGaaaaa  <  6:137721/139‑1 (MQ=255)
                                                                                                                               tGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGCATAGACCGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGGAAAGATAAAAGCCAACCTTTGATTGaaaaa  >  4:66078/1‑138 (MQ=255)
                                                                                                                                    |                                                                                                                                    
CTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTCCGGCGTCACCCAGTTCATCAAAAAGACCATCCACACCCGTAAAGATACAAGCAAACGTTTCATAGAAAAA  >  CP006881/3724911‑3725175

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