Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,647,572 Δ1 bp coding (1150/1335 nt) oprB‑II ← carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,5720C.97.2% 152.1 / ‑2.0 36coding (1150/1335 nt)oprB‑IIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base . (16/19);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

TCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCC  >  NC_002947/1647430‑1647705
                                                                                                                                              |                                                                                                                                     
tCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATagag                                                                                                                                 >  3:109618/1‑148 (MQ=255)
      tGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCcggc                                                                                                                           <  1:108195/148‑1 (MQ=255)
          acacCCCGCCTGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGcc                                                                                                                       >  4:37380/1‑148 (MQ=255)
               ccgcccgGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGt                                                                                                                   >  4:58844/1‑147 (MQ=255)
                    cgGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCa                                                                                                             <  1:20377/148‑1 (MQ=255)
                                gTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCt                                                                                                 <  7:61286/148‑1 (MQ=255)
                                   cTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTg                                                                                               >  2:39417/1‑147 (MQ=255)
                                    tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGc                                                                                              <  2:31424/147‑1 (MQ=255)
                                    tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCg                                                                                             >  3:131723/1‑148 (MQ=255)
                                    tGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCg                                                                                             >  5:35404/1‑148 (MQ=255)
                                       aGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATa                                                                                          <  1:84712/148‑1 (MQ=255)
                                        ggTTGGGGCGTACCGTGAGCCAATCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAg                                                                                         <  3:133646/148‑1 (MQ=255)
                                               gCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGAcc                                                                                                              >  4:77130/1‑120 (MQ=255)
                                               gCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGt                                                                                  >  5:47318/1‑148 (MQ=255)
                                                cGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCa                                                                                                             <  3:77130/120‑1 (MQ=255)
                                                     cGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGt                                                                            <  8:97897/148‑1 (MQ=255)
                                                       tGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGGTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGc                                                                          <  5:127790/148‑1 (MQ=255)
                                                                   ccAAGGGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGCCCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAgg                                                              >  7:138202/1‑148 (MQ=255)
                                                                            tGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCg                                                     <  6:47318/148‑1 (MQ=255)
                                                                             gCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGa                                                    >  8:120650/1‑148 (MQ=255)
                                                                             gCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGa                                                    <  7:133770/148‑1 (MQ=255)
                                                                               cGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGTGTTGACGTGCACGCTGGCCAGGGCGAAACCGAt                                                   <  8:101722/147‑1 (MQ=255)
                                                                               cGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATg                                                  <  6:35404/148‑1 (MQ=255)
                                                                                        gCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTg                                         <  8:2144/148‑1 (MQ=255)
                                                                                          tCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGCGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGgc                                       >  4:21378/1‑148 (MQ=255)
                                                                                                cTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCgg                                   >  4:36096/1‑147 (MQ=255)
                                                                                                    aCTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATc                              <  8:103835/147‑1 (MQ=255)
                                                                                                               tGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCttt                  <  1:31280/148‑1 (MQ=255)
                                                                                                               tGCACTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCtt                   >  3:151088/1‑147 (MQ=255)
                                                                                                                  aCTGGCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGaa               <  5:117697/148‑1 (MQ=255)
                                                                                                                      gCAGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATa             >  1:141493/1‑146 (MQ=255)
                                                                                                                        aGGAAGCCCGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGACACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAg         >  6:116862/1‑148 (MQ=255)
                                                                                                                              cccGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgc   <  4:27760/148‑1 (MQ=255)
                                                                                                                               ccGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc  >  8:104155/1‑148 (MQ=255)
                                                                                                                               ccGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc  <  4:151088/148‑1 (MQ=255)
                                                                                                                               ccGGGTTGTCGTAGTATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTcctgcc  >  1:136521/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                     
TCGACCTGCGACACCCCGCCCGGGTGGCGGATGTACTGCAGGTTGGGGCGTACCGTGAGCCAGTCGGCCAAGTGAATGCCGTAATACAGCTCGGCGCTGTACTCGGTGTCCTGCACTGGCAGGAAGCCCGGGTTGTCGTAGTCATAGAGGCCGGCGGCCTGGTTGACCAGGCGGGCGTTCTTGCGATAGGCAGGGTTGACGTGCACGCGGGCCAGGGCGAAACCGATGTCGTCCTTGGCGCGGGCATCGAAAGGCCCTTTGAATACCAGTCCTGCC  >  NC_002947/1647430‑1647705

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