Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 910,666 A→G 38.5% intergenic (‑50/+149) CJ019_00872 ← / ← oprD_4 hypothetical protein/Porin D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143910,6660AG38.5% 25.5 / 43.5 39intergenic (‑50/+149)CJ019_00872/oprD_4hypothetical protein/Porin D
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base G (10/5);  total (23/16)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.49e-01

AGGCAACTGCCCAGGTTGAGGATCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCGAAGCCTTCTTCGCGGGTAACC  >  ACB143/910543‑910791
                                                                                                                           |                                                                                                                             
aggcaaCTGCCCAGGTTGAGGATCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTc                                                                                                                   <  7:242304/136‑1 (MQ=255)
aggcaaCTGCCCAGGTTGAGGATCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTAGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTc                                                                                                                   <  3:119977/136‑1 (MQ=255)
     aCTGCCCAGGTTGAGGATCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGTGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACa                                                                                                              <  8:44543/136‑1 (MQ=255)
                   ggAGCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGGGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACg                                                                                                <  3:131927/136‑1 (MQ=255)
                      tCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             <  4:287152/136‑1 (MQ=255)
                        aGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGc                                                                                           <  2:87422/136‑1 (MQ=255)
                        aGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGc                                                                                           <  1:103830/136‑1 (MQ=255)
                        aGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGc                                                                                           <  6:111374/136‑1 (MQ=255)
                        aGCTCGGGGAGCGACCTGTAGACTGCGAAACGCGGGGATTGGGGGGGCATGCGCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGc                                                                                           <  3:40197/136‑1 (MQ=255)
                           tCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             >  6:257905/1‑131 (MQ=255)
                           tCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             <  5:257905/131‑1 (MQ=255)
                           tCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACg                                                                                        >  7:35285/1‑136 (MQ=255)
                                 tGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTAc                                                                                  >  3:135430/1‑136 (MQ=255)
                                  gCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTAcc                                                                                 >  6:73441/1‑136 (MQ=255)
                                      ccTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTAc                                                                             >  6:278834/1‑136 (MQ=255)
                                      ccTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTAc                                                                             >  7:242357/1‑136 (MQ=255)
                                      ccTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTAc                                                                             >  8:113178/1‑136 (MQ=255)
                                       cTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTTTACCGTACt                                                                            >  8:261465/1‑136 (MQ=255)
                                            gACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCAGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTc                                                                        >  7:180337/1‑135 (MQ=255)
                                                  gAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGAGATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGc                                                                 >  1:17188/1‑136 (MQ=255)
                                                   aaaCGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGcc                                                                <  2:145543/136‑1 (MQ=255)
                                                    aaCGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCg                                                               >  8:225984/1‑136 (MQ=255)
                                                    aaCGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCg                                                               >  2:175087/1‑136 (MQ=255)
                                                    aaCGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCg                                                               <  3:167192/136‑1 (MQ=255)
                                                    aaCGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTTTCCGACAACATTTGTCCGACCGAAGCACCGTGGACCGGACTTTTCGATTTGCCg                                                               >  8:105318/1‑136 (MQ=255)
                                                              ttGGGTGGGCATGCTCGCTACTGCGTGGGCAGGTCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTg                                                     >  5:84040/1‑136 (MQ=255)
                                                                        aTGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCggg                                           >  1:189066/1‑136 (MQ=255)
                                                                          gCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGtt                                         >  6:44991/1‑136 (MQ=255)
                                                                                        gggCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             <  5:244058/70‑1 (MQ=255)
                                                                                        gggCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             >  6:244058/1‑70 (MQ=255)
                                                                                         ggCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGCGCGGGTTTACCCGCGAACACCg                          >  1:42104/1‑136 (MQ=255)
                                                                                                 cATTCGCGGCGGCGATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGCGCGGGTTTTCCCGCGAACACCGGCGCAGcc                  >  4:279841/1‑136 (MQ=255)
                                                                                                  aTTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGGACACCGGCGCAGccg                 >  3:134279/1‑135 (MQ=255)
                                                                                                      gcggcggcgATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCGAAGCcttct              >  7:52504/1‑135 (MQ=255)
                                                                                                             cgATATTCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCCGGTTTAACCGCGAACACCCGCGAAGCCTTATTCGCggg       >  6:93141/1‑135 (MQ=255)
                                                                                                             cgATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCGAAGCCTTCTTCGCGGGt      <  7:284847/136‑1 (MQ=255)
                                                                                                             cgATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCGAAGCCTTCTTCGCGGGt      <  1:234026/136‑1 (MQ=255)
                                                                                                             cgATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCCCGAACACCCCCCAAGACGTCTTAGCggg       >  3:74708/1‑135 (MQ=255)
                                                                                                                  ttCTAGGGAGGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCGAAGCCTTCTTCGCGGGTAAcc  <  3:274759/135‑1 (MQ=255)
                                                                                                                          aaGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             <  7:232642/36‑1 (MQ=255)
                                                                                                                          aaGCGGGCATTGTCGGACAAGATTTGTCCGACGGaa                                                                                             >  8:232642/1‑36 (MQ=255)
                                                                                                                           |                                                                                                                             
AGGCAACTGCCCAGGTTGAGGATCAGCTCGGGGTGCGACCTGTAGACTGCGAAACGCGGGGATTGGGTGGGCATGCTCGCTACTGCGTGGGCAGGGCCATTCGCGGCGGCGATATTCTAGGGAAGCGGGCATTGTCGGACAAGATTTGTCCGACGGAAGCACGGTGTACCGTACTTTTCGATTTGCCGTGGACCCTTGTAGGAGCGGGTTTACCCGCGAACACCGGCGAAGCCTTCTTCGCGGGTAACC  >  ACB143/910543‑910791

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