Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 2,801,996 G→A 100% A327T (GCC→ACC)  U712_14440 → UPF0118 membrane protein ytvI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068812,801,9960GA99.9% 120.2 / ‑7.7 42A327T (GCC→ACC) U712_14440UPF0118 membrane protein ytvI
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (16/24);  minor base C (1/0);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 4.15e-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.

AAGCAATCGGAATTGGTATCCTGTACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCGCCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAATATGGTCCTATATCACCGTTCAGC  >  CP006881/2801886‑2802125
                                                                                                              |                                                                                                                                 
gacagaTCGGAATTGGTATCCTGTACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCgg                                                                                                       <  3:22209/134‑1 (MQ=255)
aaGCAATCGGAATTGGTATCCTGTACCTCGTCGTTCTTATTCAGCGGCATCTGACAGAGCCAAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCgg                                                                                                       <  4:209031/139‑1 (MQ=255)
aaGCAATCGGAATTGGTATCCTGTACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCgg                                                                                                       <  1:340265/139‑1 (MQ=255)
            ttGGTATCCTGTACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACt                                                                                           <  5:331946/139‑1 (MQ=255)
                tATCCTGTACCTCGTCGTGAGTATTTAGCCGGGGGTGACAGAGCCCCACATACTAAGTAAAACTAAGGGGTACGATCCGCTCGGCGCACTCCTCTCCCATTTTTCCGGCGTTAAATTAATACGCGTTTTAGGAcctgag                                                                                       >  4:22209/1‑134 (MQ=255)
                       tACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCggg                                                                                  <  3:117018/137‑1 (MQ=255)
                       tACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCggg                                                                                  <  2:186548/137‑1 (MQ=255)
                       tACCTCGTCGGTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGAGCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGGTTTTAGGACTGATTGCggg                                                                                  <  1:252578/137‑1 (MQ=255)
                            cgtcgtTCTTATTCAGCGGCAGCTGCCAGCGCCCAAGATACTAAGTAAATCTATCGGGATCGCTCCGCAGGCAACCCGCATAACCCTGTTTGCCCGCTGTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACAGGcg                                                                           >  6:185906/1‑138 (MQ=255)
                              tcgtTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGt                                                                         >  5:68732/1‑139 (MQ=255)
                              tcgtTCTTATTCAGCGGCAGCTGACAGACCCCAAGAAACTAATTAAATCTATAGGGATCGATCCGCTGGCAACCCTCATACCCCAGTTTGCCGGCTTTTAATTATTAGGGTTTTTATGACTGATTGCGGGACCGGCTg                                                                          >  4:117018/1‑138 (MQ=255)
                              tcgtTATTATTAAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATAGATCCGCTGGCCACACTCCTCACACTGTTTGCCTGCGTTAAATTTTTAGGGTTTTTAGGCAGTATTGGGGGACCGGCTg                                                                          >  2:89872/1‑138 (MQ=255)
                                gtTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTAc                                                                       <  2:146827/139‑1 (MQ=255)
                                gtTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTAc                                                                       <  4:55364/139‑1 (MQ=255)
                                      aTTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGtc                                                                  <  1:356835/138‑1 (MQ=255)
                                                    gACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTAAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCttt                                                    >  4:12966/1‑138 (MQ=255)
                                                           cccAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGc                                                                                                                         >  8:155114/1‑62 (MQ=255)
                                                           cccAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGc                                                                                                                         <  7:155114/62‑1 (MQ=255)
                                                                 aTACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGgcg                                      <  2:307561/139‑1 (MQ=255)
                                                                     tAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGTGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCacaac                                           >  5:206861/1‑130 (MQ=255)
                                                                     tAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGTGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCacaac                                           <  6:206861/130‑1 (MQ=255)
                                                                     tAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCt                                                      >  1:298250/1‑119 (MQ=255)
                                                                     tAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCt                                                      <  2:298250/119‑1 (MQ=255)
                                                                      aaGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGCCCGGCTGTACTCGTCCTCATCCACGCTTTCATCCCAACAGGCGCAtt                                  >  8:339813/1‑138 (MQ=255)
                                                                       aGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTaa                                <  1:69408/139‑1 (MQ=255)
                                                                           aaTCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAg                             >  4:175170/1‑138 (MQ=255)
                                                                           aaTCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAg                             >  2:158962/1‑138 (MQ=255)
                                                                           aaTCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAg                             <  6:68732/138‑1 (MQ=255)
                                                                            aTCTATCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGa                            >  4:324825/1‑138 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCgg                                                                                                       >  5:205655/1‑58 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCgg                                                                                                       <  6:205655/58‑1 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGtc                                                                  <  1:134938/95‑1 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGtc                                                                  >  2:134938/1‑95 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAata                        <  3:3123/137‑1 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAata                        >  6:183191/1‑137 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAata                        <  5:183191/137‑1 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAata                        <  3:12966/137‑1 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAata                        <  1:90982/137‑1 (MQ=255)
                                                                                 tCGGGATCGATCCGCTGGCAACACTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAATATg                      >  1:243483/1‑139 (MQ=255)
                                                                                                     acacTCATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAGAAGAAATATGGTCCTATATCACCGTTCAGc  <  1:89872/139‑1 (MQ=255)
                                                                                                          cATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTaa                                >  2:280984/1‑104 (MQ=255)
                                                                                                          cATCACCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTaa                                <  1:280984/104‑1 (MQ=255)
                                                                                                              |                                                                                                                                 
AAGCAATCGGAATTGGTATCCTGTACCTCGTCGTTCTTATTCAGCGGCAGCTGACAGAGCCCAAGATACTAAGTAAATCTATCGGGATCGATCCGCTGGCAACACTCATCGCCCTGTTTGCCGGCTTTAAATTATTCGGGTTTTTAGGACTGATTGCGGGACCGGCTGTACTCGTCATCATCCAAGCTTTCATCACAACAGGCGCATTAAAAGAAATATGGTCCTATATCACCGTTCAGC  >  CP006881/2801886‑2802125

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: