Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,724,980 T→C 100% K63K (AAA→AAG U712_19205 ← Uncharacterized protein ywcI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,724,9800TC96.2% 168.7 / ‑5.5 51K63K (AAA→AAGU712_19205Uncharacterized protein ywcI
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (11/37);  minor base G (3/0);  total (14/37)
Fisher's exact test for biased strand distribution p-value = 1.75e-02
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.

CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724851‑3725105
                                                                                                                                 |                                                                                                                             
cTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTcc                                                                                                                      <  3:159417/139‑1 (MQ=255)
cTTTTGACTGGGATTATACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTttcttc                                                                                                                       <  4:291009/138‑1 (MQ=255)
       cTGGGATTTAACCTCCTTGTGACCGCCATGGAAAGAGGCGGGCCCGCGCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTca                                                                                                               <  5:219606/139‑1 (MQ=255)
          ggATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAg                                                                                                            >  7:98911/1‑139 (MQ=255)
          ggATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGGATTATTCTTTTCTTACACATTAACAg                                                                                                            >  2:207385/1‑139 (MQ=255)
               ttACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCAt                                                                                                        <  1:257353/138‑1 (MQ=255)
                     cctTTTGACCGCCATGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCATTCGTCGTTTCTTGCTGATTCACAGCCCACCGCCat                                                                                                 >  2:41979/1‑138 (MQ=255)
                      ctTTTGACCGCCATGTACAGATGCGGGCCCGCACTCTCCTGCTTTGCCGAGTCCTTGAAGGCTTGCAACCAGATTCTTATATATAGCGCGGCGGCTTTCGCTTTCTGCTTTTGGTCCTGATTCACCGCCAAACGcccta                                                                                                >  6:64229/1‑136 (MQ=255)
                                   tGTAAAGATGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                   <  6:254436/139‑1 (MQ=255)
                                          aggCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  1:207385/137‑2 (MQ=255)
                                          aggCGGGCCCGCGCTCTCCGGCGTGGGCGATTCCTTGCAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACCGCCATCCGCCATGCATAAAACGAAACGAGTcct                                                                            <  6:255654/137‑2 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCTGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  4:253458/139‑2 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  1:36770/139‑2 (MQ=255)
                                          aTGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  1:41979/139‑2 (MQ=255)
                                           tGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTGGCGTTGGGGTTTTGTTGCTCATTCACAGCCCTCCGCCATTCATAAAACGAAACGAGTcctg                                                                           >  2:7427/1‑137 (MQ=255)
                                           tGCGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGCGTTGTGGGTTTATTACTCCTTAACAGCCCTACGCCCGTCATAAAACGAAACCAGTcctg                                                                           >  1:125917/1‑137 (MQ=255)
                                            gCGGGCCCGATCTCTCCTGCTTTGGCGATTCCTGGAAGGCTTGCAAACATCTTCTTTTATAATGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCAGCCATTCATAAAACGAAACGAGTcct                                                                            <  4:82914/137‑2 (MQ=255)
                                             cGGGCCCGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTCGCTTTATGGTTTTATTCCACATTAACAGCCATCCCCCATTCATAAAACCAAACGAGTCCagcg                                                                         >  4:214229/1‑139 (MQ=255)
                                              gggCCCGCTCTCTCCTGCTTTGGCGATTCCTTGCAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCGTTGGCGTTGATCGTTTCTTGCGCATTCCACGCCAACCGCCCTTCAGATAACGCAACGATTCCTGCgg                                                                        >  7:351555/1‑139 (MQ=255)
                                               ggCCCGCTCGCGCCTGCTTTGGCCTTTCCTTGAAGGCTTGCAAACATCTTCTTTTATAACGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATACACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGc                                                                       <  6:24266/139‑1 (MQ=255)
                                                     ctctctCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCAc                                                                 <  3:49371/139‑1 (MQ=255)
                                                            cTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTAGCTTTATGCTTTTATTCCTCATTAACAGCCATCCGCCAATCATAAAACCAAACCAGTCCTGCCGCCTGCAAagccgcc                                                          >  8:28091/1‑138 (MQ=255)
                                                                 ttGGCGATTCCTTGAAGGCTTGCAAACATCTGCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCCCCATCCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCt                                                     <  4:108913/139‑1 (MQ=255)
                                                                  tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                   >  3:283581/1‑108 (MQ=255)
                                                                  tGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                   <  4:283581/108‑1 (MQ=255)
                                                                   ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  3:214229/139‑1 (MQ=255)
                                                                   ggCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCctttcttt                                                   <  8:263802/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                                    <  4:180928/139‑1 (MQ=255)
                                                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  1:319108/139‑1 (MQ=255)
                                                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  8:181386/139‑1 (MQ=255)
                                                                                  gCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  2:129576/139‑1 (MQ=255)
                                                                                  gCTTGCAAAAATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAg                                    <  6:90332/139‑1 (MQ=255)
                                                                                     tGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTGTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                   >  3:115820/1‑89 (MQ=255)
                                                                                     tGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAAc                                                                                   <  4:115820/89‑1 (MQ=255)
                                                                                       cAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTGGTATGTACGCTGAGCCGaa                               <  5:84371/139‑1 (MQ=255)
                                                                                             tcttctTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            <  3:143385/88‑2 (MQ=255)
                                                                                             tcttctTTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTcct                                                                            >  4:143385/1‑87 (MQ=255)
                                                                                                 cttTTATAAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACCGCCGCCGTTGTTTTGTATTTTAGCAGAGCCGAATAGACCGccg                     >  5:173606/1‑137 (MQ=255)
                                                                                                      ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTgg                                                 >  3:345736/1‑106 (MQ=255)
                                                                                                      ataAAGCGCGGCGGCTTTCGCTTTCTGCTTTTCTTCCTCATTCACAGCCATCCGCCATTCATAAAACGAAACGAGTCCTGCGGCCTGCACAGCCGCCTTTCTTTgg                                                 <  4:345736/106‑1 (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)
                                                                                                                                 |                                                                                                                             
CTTTTGACTGGGATTTTACCTCCTTTTGACCGCCATGTAAAGATGCGGGCCAGCTCTCTCCTGCTTTGGCGATTCCTTGAAGGCTTGCAAACATCTTCTTTTATAAAGCGCGGCGGCTTTCGCTTTTTGTTTTTCTTCCTCATTCACAGCCATCCGCCACTCATAAAACGAAACAAGTCCCGCGGCCTGCACCGCCGCTTTTCTTTGGTATGTACGCTGAGCCGAATAGACAGCGGCAGCCGCAGTTTTTCTGTC  >  CP006881/3724851‑3725105

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: