Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,991 G→T 100% Q80K (CAG→AAG)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,9910GT88.3% 42.8 / ‑1.4 17Q80K (CAG→AAG) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (5/10);  total (7/10)
Fisher's exact test for biased strand distribution p-value = 1.54e-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: Not supported by required number of reads on both strands.

ACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACA  >  CP009273/2305849‑2306128
                                                                                                                                              |                                                                                                                                         
aCGGCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTc                                                                                                                                     <  1:60104/149‑1 (MQ=255)
       tCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTgg                                                                                                                              <  1:58620/149‑1 (MQ=255)
            gAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATAATGGATTTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTa                                                                                                                         >  1:346018/1‑149 (MQ=255)
             aaGTTACGGCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTgg                                                                                                                              <  2:36713/143‑1 (MQ=255)
                ttACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCg                                                                                                                     <  1:71288/149‑1 (MQ=255)
                        tAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTg                                                                                                             >  1:71571/1‑149 (MQ=255)
                                  cGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCa                                                                                                         <  2:174406/143‑1 (MQ=255)
                                       gtACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTCAc                                                                                                    >  2:166874/2‑143 (MQ=255)
                                          cGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTc                                                                                           <  1:13731/149‑1 (MQ=255)
                                                   tCGAAAGAACCCACATCCTGGAATTTCATACCTGCTAATTCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGaa                                                                                  <  1:140964/149‑1 (MQ=255)
                                                             ccACATCCTGGAATTTCAGACCTGCGAAGGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGc                                                                        >  1:240866/1‑149 (MQ=255)
                                                                                          gCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTAc                                                 <  2:280909/143‑1 (MQ=255)
                                                                                                                tttcgttttcAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAgg                     <  1:143615/149‑1 (MQ=255)
                                                                                                                     ttttcAGCGCTGTTGCCCTGGATCTTCTATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTCCATCTTTGTTGTCCGACACATAGTGCAg                      >  2:75553/1‑143 (MQ=39)
                                                                                                                     ttttcAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAGGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCGGAGAAATAGGGCAg                      >  2:310147/1‑143 (MQ=255)
                                                                                                                     ttttcAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACAACTTGGTTGTCAGAGAAAAAGTGCAg                      >  2:102512/1‑143 (MQ=255)
                                                                                                                                   gCCATGGATCTTATATTCCCACTGGCCGTATCCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACa  <  1:75324/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                         
ACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACATCTTTGTTGTCAGAGAAATAGTGCAGGCCGTCTACTTTACCGTACA  >  CP009273/2305849‑2306128

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: