Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 712,137 C→T 9.0% R154H (CGC→CAC)  U712_03445 ← Uncharacterized protein yezG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP006881712,1370CT 9.0% 127.2 / 2.3 44R154H (CGC→CAC) U712_03445Uncharacterized protein yezG
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base T (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

GCGTTTTTGCCTTTACCCGTTGCTTCTTGGAGCTTTTGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCTTTTTTCTGTTCGGACATGATGGATTTAGATTTTTTAT  >  CP006881/712002‑712275
                                                                                                                                       |                                                                                                                                          
gCGTTTTTGCCTTTACCCGTTGCTTCTTGGCGCTTTTGGCAGTCGGATTTGATGGATTGTTGGTTTGGTTCAGCCAGGGCGTATTCTGTTTTCAGGTGTTCGTACAAGGCGTGCAATTTTTGCAGGGTTTTTTTGCggt                                                                                                                                         >  3:156853/1‑137 (MQ=255)
       tGCCTTTACCCGTTGCTTCTTGGAGCTTTTGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGt                                                                                                                                   <  2:263849/138‑1 (MQ=255)
                   ttGCTTCTTGGAGCTTTTGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTTTCGTTGGCATCGGCCAg                                                                                                                      >  3:101434/1‑139 (MQ=255)
                      cttcttGGAGCTTTTGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTc                                                                                                                   <  3:209154/139‑1 (MQ=255)
                         cttGGAGCTTTTGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTc                                                                                                                >  7:127021/1‑139 (MQ=255)
                                   ttGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTAGTTGGCATCGGCCAGCTCGTCTTTGAAGGtt                                                                                                      >  1:172160/1‑139 (MQ=255)
                                     ggAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGttt                                                                                                     <  1:366800/138‑1 (MQ=255)
                                     ggAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGttt                                                                                                     <  4:101434/138‑1 (MQ=255)
                                              aTTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGaaaa                                                                                           <  1:188843/139‑1 (MQ=255)
                                                        ttGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGc                                                                                 >  2:258018/1‑139 (MQ=255)
                                                                         aCAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGat                                                                <  6:322092/139‑1 (MQ=255)
                                                                         aCAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGat                                                                <  5:94706/139‑1 (MQ=255)
                                                                         aCAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGat                                                                <  8:176114/139‑1 (MQ=255)
                                                                             ggCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTTCGTTTGTCGTTGGCATCGGCCAGCTAGTCTTTGAAGGTTTCTGTTGAAAACAGGTCTAGCGGGAGGATGTCGTCGATATCg                                                            >  1:316894/1‑139 (MQ=255)
                                                                             ggCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTATGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGatttag                                                            >  6:314388/1‑135 (MQ=255)
                                                                              gtgTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGt                                                           >  4:243625/3‑139 (MQ=255)
                                                                              gCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGata                                                               <  2:6002/135‑1 (MQ=255)
                                                                              gCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGata                                                               >  1:6002/1‑135 (MQ=255)
                                                                              gCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGt                                                           <  1:74728/139‑1 (MQ=255)
                                                                                gTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGttt                                                         <  4:156853/139‑1 (MQ=255)
                                                                                   ttCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCa                                                                                                                              >  2:212825/1‑67 (MQ=255)
                                                                                   ttCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCa                                                                                                                              <  1:212825/67‑1 (MQ=255)
                                                                                   ttCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAgg                                                      <  3:255209/139‑1 (MQ=255)
                                                                                   ttCTGGTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGTGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAgg                                                      <  5:223774/139‑1 (MQ=255)
                                                                                       gTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTACAGTGTTTTTTTGCGTTTGTCGTTGGCCTCGGCCAGCTCCGGTTTGTAGGTTTGTGGTGAAAACAGGTCGAGCGGGAGGTTGTGGGGGATATCGTTTAtgttgt                                                  >  2:62723/1‑139 (MQ=255)
                                                                                        ttttCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGTGTTTTTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTATGTTGAAAAAATGTCGAGCGTGAGGATGTCGTCGATATCGTTTAGGATGTc                                                 >  6:223774/1‑139 (MQ=255)
                                                                                        ttttCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGTGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAgg                                                      <  5:218366/134‑1 (MQ=255)
                                                                                        ttttCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGTGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAgg                                                      >  6:218366/1‑134 (MQ=255)
                                                                                             aGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCa                                            <  2:172160/139‑1 (MQ=255)
                                                                                                  ttCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCt                                       >  4:43537/1‑139 (MQ=255)
                                                                                                  ttCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCt                                       >  5:95456/1‑139 (MQ=255)
                                                                                                  ttCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCt                                       >  6:361374/1‑139 (MQ=255)
                                                                                                    cGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAgg                                                      >  8:54376/1‑122 (MQ=255)
                                                                                                    cGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAgg                                                      <  7:54376/122‑1 (MQ=255)
                                                                                                          aGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCa                                            >  6:169452/1‑126 (MQ=255)
                                                                                                          aGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCa                                            <  5:169452/126‑1 (MQ=255)
                                                                                                           ggCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGata                                                               >  7:109577/1‑106 (MQ=255)
                                                                                                           ggCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGata                                                               <  8:109577/106‑1 (MQ=255)
                                                                                                               tCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGAt                                                                         <  7:26198/92‑1 (MQ=255)
                                                                                                               tCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGAt                                                                         >  8:26198/1‑92 (MQ=255)
                                                                                                               tCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCTTTTTTCTGTTCgg                          >  5:17131/1‑139 (MQ=255)
                                                                                                               tCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATAGTTTATGATGTCGTTCATCGGGTTTTTGTGTTCgg                          >  5:51535/1‑139 (MQ=255)
                                                                                                                              gtTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCTTTTTTCTGTTCGGACATGATGGAtt              >  2:336842/1‑136 (MQ=255)
                                                                                                                                  ttttGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCTTTTTTCTGTTCGGACATGATGGATTTAGAttt       >  7:61322/1‑139 (MQ=255)
                                                                                                                                       cGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCTTTTTTCTGTTCGGACATGATGGATTTAGATTTTTTAt  <  7:201754/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                          
GCGTTTTTGCCTTTACCCGTTGCTTCTTGGAGCTTTTGGAAGTCGGATTTGATGAATTGTTCGTTTGGTTCAGACAGGGCGTATTCTGTTTTCAGGTCTTCGTCCAAGGCGTCCAATTTTTCCAGTGTTTTTTTGCGTTTGTCGTTGGCATCGGCCAGCTCGTCTTTGAAGGTTTCTGTTGAAAACAGGTCGAGCGGGAGGATGTCGTCGATATCGTTTAGGATGTCTTTCATCGCTTTTTTCTGTTCGGACATGATGGATTTAGATTTTTTAT  >  CP006881/712002‑712275

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: