Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,305,945 C→G 69.3% R95P (CGT→CCT)  ompC ← outer membrane porin protein C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,305,9450CG69.3% 54.8 / 33.6 39R95P (CGT→CCT) ompCouter membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base C (3/9);  new base G (9/18);  total (12/27)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

TGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACAT  >  CP009273/2305805‑2306083
                                                                                                                                            |                                                                                                                                          
tGAAGTTGTCAGAACCGTAGCTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACgg                                                                                                                                          <  1:150057/143‑1 (MQ=255)
     ttGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAgtt                                                                                                                               >  2:437054/1‑149 (MQ=255)
             aCCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTTCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTgttttcgt                                                                                                                       <  2:14673/149‑1 (MQ=255)
             aCCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTgttttcgt                                                                                                                       <  2:393314/149‑1 (MQ=255)
                   ggTGTCACCACCGAATGCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTgttttcgt                                                                                                                       <  1:445713/143‑1 (MQ=255)
                    ctgtCACCACCGAATTCTGGCAGCACCTCGGTGCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCACGAGTTgttttcgtt                                                                                                                      <  1:207887/142‑1 (MQ=255)
                       tcgccaccGAAGTCGGGCAGTACGTCGGTCCAGGAAGTTCCGTCTTAAACATCGCCGGAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCACACCTGCTAATGCCACAGGCGTCCAGGAGTTgttttcgttttc                                                                                                                   <  1:376255/140‑1 (MQ=255)
                           caccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAgcg                                                                                                               >  1:544105/1‑143 (MQ=255)
                              cGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTg                                                                                                            >  1:381592/1‑143 (MQ=255)
                                    cTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATc                                                                                                <  2:215099/149‑1 (MQ=255)
                                     tGGCAGTACGTCGGTCCGGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCACACCTGCGTATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCt                                                                                                     <  1:355869/143‑1 (MQ=255)
                                     tGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCt                                                                                               <  2:486253/149‑1 (MQ=255)
                                      ggCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTg                                                                                              <  2:192302/149‑1 (MQ=255)
                                       gaaGTACGTCGGGCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTgg                                                                                                   <  1:13992/141‑1 (MQ=255)
                                       gCAGTACGGCGGTCCAGGAAGTTACGTCATAGACAACGCCGTAGTTACGACCGTAGTCGTAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTgg                                                                                                   <  1:427427/143‑1 (MQ=255)
                                         aGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGata                                                                                           >  2:605982/1‑149 (MQ=255)
                                         aGTACGCCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCCGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGAt                                                                                                 <  1:523812/143‑1 (MQ=255)
                                             cggcGGCCCAGGAAGTGACGTCGTAACCAACGCGGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGa                                                                                             <  1:23665/140‑1 (MQ=255)
                                             cGTCGGTCCAGGAAGTTACGGCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCtta                                                                                             <  1:70946/143‑3 (MQ=255)
                                              gTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCGGat                                                                                            >  1:434860/1‑143 (MQ=255)
                                                  gTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTg                                                                                  >  2:209537/1‑149 (MQ=255)
                                                       ggAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCGCTGGCCGt                                                                             >  2:34288/1‑149 (MQ=255)
                                                       ggAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACATGCGAATGCCACCGGGGACCAGGAGCTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATCTTCCCCCTGGCCgg                                                                             >  2:246102/1‑148 (MQ=255)
                                                         aaGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTgg                                                                                 <  1:485031/143‑1 (MQ=255)
                                                               cGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTc                                                                     >  2:638226/1‑149 (MQ=255)
                                                                   aTAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCt                                                                 <  2:390147/149‑1 (MQ=255)
                                                                       acaacGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTGCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAg                                                                   <  1:183825/143‑1 (MQ=255)
                                                                              cGTAGTTACGAGCGGAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCa                                                            <  1:396448/143‑1 (MQ=255)
                                                                                        aCCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCAc                                            <  2:229168/149‑1 (MQ=255)
                                                                                              gTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCAGCTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCAc                                            <  1:376215/143‑1 (MQ=255)
                                                                                               tCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCAcc                                           >  1:245609/1‑143 (MQ=255)
                                                                                                         ccACCTCCTGGACTTTCAGACGTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTTTTGCTCTGGTTCTTATATTCCCACCGGCCGTAACCGGTCAGCTGGCCAGTATCCTGAGtttttcggatttagcga                                 <  1:236006/143‑15 (MQ=255)
                                                                                                         ccACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTTATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCa                                 <  1:208654/143‑1 (MQ=255)
                                                                                                           acaTCCTGGAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCAt                         <  2:602011/149‑1 (MQ=255)
                                                                                                                  ggAATTTCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTAGGGAGCAACACTGATGACCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTc                  >  2:628482/1‑149 (MQ=255)
                                                                                                                       ttCAGACCTGCGAATGCCACAGGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTc             <  2:286640/149‑1 (MQ=255)
                                                                                                                       ttCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCTCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGt                   <  1:228064/143‑1 (MQ=255)
                                                                                                                               tGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGTTCAGCTGGTCAGTAACCTGAGTTCCGCAGCGCGAGCCAAGACGCATGTAGGTCTGGTCGc           <  1:405544/143‑1 (MQ=255)
                                                                                                                                        cacaCGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACAt  >  1:419917/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                          
TGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCTACAT  >  CP009273/2305805‑2306083

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