Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 868,330 T→G 43.9% intergenic (+19/+27) loiP → / ← mneA Metalloprotease LoiP/Putative manganese exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122868,3300TG43.9% ‑2.5 / 3.9 20intergenic (+19/+27)loiP/mneAMetalloprotease LoiP/Putative manganese exporter
Reads supporting (aligned to +/- strand):  ref base T (8/2);  new base G (7/2);  total (16/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CACCGCCTGAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCA  >  ACB122/868202‑868445
                                                                                                                                |                                                                                                                     
cACCGCCTGAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATgg                                                                                                                 >  8:296149/1‑135 (MQ=255)
 aCCGCCTGAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCtggctggc                                                                                                                >  2:196614/1‑135 (MQ=255)
 aCCGCCTGAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGACAAGCCGTAAATGCCTGAATCCTGGGGCTGCagtgc                                                                                                                >  7:176780/1‑131 (MQ=255)
           ttCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCCAGGCAGcccccccc                                                                                                      >  4:130743/1‑134 (MQ=255)
           ttCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGAGTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAACGGCCTGAATCCTGGGGCGGCCTGGCCGcccccccc                                                                                                      >  3:79613/1‑134 (MQ=255)
            tCATGAGTACTCCCCCGGCGGCCGCGAGCCGCATTGCTGCGCTGCAGGCGGCGATTCCGACGGTGATGCCGCGGTATCCGGCGGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCCgg                                                                                                    <  7:294916/135‑1 (MQ=255)
                gAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGCCTGGGCGccccccccgggtt                                                                                                 >  3:297403/1‑129 (MQ=255)
                  gTACTCACCCGGCGCCCGCGAGCCGCATGGCTTCGCGGCCGGGGGCTATTCCCAAGGTTAGGCCGCTGTATCAGGCGGCCGCCAAGCCGTAAATGCCTTAATCGTGGGGCTGCATGGCAGCCCCCCCGGTTTTc                                                                                               <  6:119185/134‑1 (MQ=255)
                        aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCCTGGGAGCCCCCCCGGTGTTGCCTTTAc                                                                                        >  8:209475/1‑135 (MQ=255)
                        aCCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGGCTGGCCGCCCCCCCGGGTTTTCCGtttg                                                                                         >  6:281871/1‑134 (MQ=255)
                         ccgggcggcCGCGTGCCGAATTGCTGCGCTGCTGGCGGGGATTCCGAAGGTGATGCCGCTGTATCTGGCTGCCGCCAAGCGTTAAATGCCTGAATCCTGGGGCGGAATGGCAGCCCCCCCGGTTTTCCCTTTa                                                                                         <  7:296149/130‑1 (MQ=255)
                                tCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGCATGGGCCCCCCCCCGGTGTTACACTTACCCaaaacc                                                                                >  3:148833/1‑131 (MQ=255)
                                     gAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGGCTGGGAGCCCCCCCGGGTTTGCCGTTTGCCCAAGccgcc                                                                              >  1:219801/1‑130 (MQ=255)
                                                      tGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCCGGGTTTGGCGTTACCCCAACCCGCCCCTGGCCGTGATTGccca                                                          >  7:93556/1‑132 (MQ=255)
                                                        cAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGCCTTGGAGCCCCACCGGGTTTAGCATTTCCCCATCCCGACACAGGGCGAGATTGCCGaaa                                                        >  7:209853/1‑133 (MQ=255)
                                                                   tCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCCTGGCCGCCCCCCCGGTTTTCCCTTTTCCCCATCCACCCCCTGGCCGTAATTGGCGAATAGACGGcccccc                                            >  2:222676/1‑131 (MQ=255)
                                                                               gCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCGGGATGGCCGCCCCCCCGGGTTTTCCGTTACCCCATCCAGCCCCTGGCCGTCATTGCCGCATAGACGGCCACCCTTGCCAggaaga                                >  7:85928/1‑130 (MQ=255)
                                                                                ccGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGGATGGCCGCCCCCCCGGGTTTGCCGTTTACCAATACCCCCACGGGGCGTCAATTGCCAATAGACGGCCACGCTTGACAGttaga                                >  8:225339/1‑130 (MQ=255)
                                                                                      aagccggcggcCGCCACGCCTAAAATCCCTTAAACCTGGGGGGGCATGGCCGCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGc                         <  3:48867/130‑1 (MQ=255)
                                                                                                         cGTAAATGCCTGAATCCTGGGGGAGAAAGGCACCCCCCCCGGTTTTCGCGTTTCCCAATCCATCCACTGGCCTTGAATTCGGGAGAATCCGCCCCGAATTCCACTACGAAGCACGCCGATGCAGCCCggggggg       >  8:226921/1‑132 (MQ=255)
                                                                                                             aaTGCCTGAATCCTGGGGCTGCCTGGCAGCCCCCCCGGGTTTTCCGTTTCCCAATCCAGCCACTTGCCGTCAATTGCGAATATACGGCCCCCATTGCCCGGACGAAGAACCCCCTTGCCCCCCCGCGGCGCATCa  >  8:4417/1‑135 (MQ=255)
                                                                                                                                |                                                                                                                     
CACCGCCTGAGTTCATGAGTACTCACCCGGCGTCCGCGAGCCGCATTGCTTCGCTGCAGGCGGCGATTCCGAAGGTGATGCCGCTGTATCAGGCTGCCGCCAAGCCGTAAATGCCTGAATCCTGGGGCTGCATGGCAGCCCCCCCGGTTTTCCCTTTACCCAATCCAGCCACTGGCCTTCATTGCCGAATAGACGGCCACGATTGCCAGTACGAAGAACGCCGTTGCCGCCAGGCGGCGAATCA  >  ACB122/868202‑868445

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