Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 120,975 Δ1 bp 100% intergenic (+272/+100) polA → / ← engB DNA polymerase I/putative GTP‑binding protein EngB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122120,9750T.93.2% 235.2 / 8.9 59intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base . (41/14);  minor base G (0/4);  total (41/18)
Fisher's exact test for biased strand distribution p-value = 6.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

GAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTC  >  ACB122/120853‑121092
                                                                                                                          |                                                                                                                     
gAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgccc                                                                                                            >  4:102780/1‑130 (MQ=255)
 aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGGACTTCCCCCTCCCCATACGAAGTGCGGGGTTTTTTTTGCCCg                                                                                                          >  8:277997/1‑134 (MQ=255)
 aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCg                                                                                                          >  3:97578/1‑134 (MQ=255)
            aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCg                                                                                                >  3:117858/1‑133 (MQ=255)
            aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCg                                                                                                >  6:115040/1‑133 (MQ=255)
                      aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGtttt                                                                                    >  8:78029/1‑135 (MQ=255)
                      aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTAGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑‑TATTTTTGCCCGGAAATAGCCGGGTGCGGGTTTtc                                                                                   >  6:193653/1‑135 (MQ=255)
                            cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgc                                                                               >  7:205306/1‑134 (MQ=255)
                             tATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTACCCCTCCCCATACGAAGTCCGGGGGGTTTTTTGCCCGGATTTTGCCGGGAGCGGGTTTTCTCGCGa                                                                             >  6:181569/1‑135 (MQ=255)
                             tATGTTGTGTTAAGTGTTGGCAGATATCTCGACCCCGCCCTAGCGGTCCGGACTTGGACCCAGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                             <  5:95937/135‑1 (MQ=255)
                                gttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                           <  8:205306/134‑1 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTTCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                    >  5:248593/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                    >  5:297126/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                    >  8:255638/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                    >  6:15984/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                    >  7:90346/1‑134 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGAAGTc                                                                   >  6:201914/1‑135 (MQ=255)
                                       tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGAGATGt                                                                    >  2:187830/1‑134 (MQ=255)
                                         aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                 >  7:115214/1‑135 (MQ=255)
                                         aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                 >  8:45772/1‑135 (MQ=255)
                                         aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                 >  8:290372/1‑135 (MQ=255)
                                            gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                               >  3:212387/1‑134 (MQ=255)
                                            gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                               >  4:225197/1‑134 (MQ=255)
                                            gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAAt                                                               >  6:247298/1‑134 (MQ=255)
                                            gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATa                                                              <  8:115214/135‑1 (MQ=255)
                                            gtTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATa                                                              >  8:81754/1‑135 (MQ=255)
                                              tGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGCACTTCCCCCTCCCCATACGAAGTCCGGGGGTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATaca                                                            >  7:284695/1‑135 (MQ=255)
                                                 cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGCATGCGATGTCGAATACAccg                                                         >  3:178830/1‑135 (MQ=255)
                                                 cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccg                                                         >  7:303823/1‑135 (MQ=255)
                                                 cAGATATCTGGACCCCGCCCTAGCGGACCGGACTTGGACCCCGAACTTCCCCCTCCCCATACCAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCCGGATTTCTCGCGAATGCGATGTCGAATACAccg                                                         >  3:85021/1‑135 (MQ=255)
                                                      atCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGTTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCa                                                    <  8:221902/135‑1 (MQ=255)
                                                        cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                   <  7:250353/134‑1 (MQ=255)
                                                        cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAt                                                   <  3:102780/134‑1 (MQ=255)
                                                            cccccGCCCTCGCGGTCCGTACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGGGGGTTTTCTCGCGATTGCGATGTCGAATACACCGCCGCATTcgcg                                              <  7:146480/134‑1 (MQ=255)
                                                             ccccGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGACGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCgg                                             <  8:160723/135‑1 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCCTACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                       >  3:34872/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgc                                               <  4:125675/127‑1 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTcgc                                               >  3:125675/1‑127 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTGAAc                                       >  1:120522/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                       >  6:224469/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                       >  5:186073/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                       >  2:32539/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                       >  7:30956/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAc                                       >  1:17269/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCAGCATACGCGGGTAAAc                                       >  8:172013/1‑135 (MQ=255)
                                                                   ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTGTTATCGCGAATGCGCTGTCGAATACACCGCCGCCTTCGCGGGGAAAc                                       >  5:183024/1‑135 (MQ=255)
                                                                    cTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGTATGCGATGGCGAATACACCGCCGCATTCGCGGGTAAAcc                                      >  1:299546/1‑135 (MQ=255)
                                                                       gCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                   <  6:186073/135‑1 (MQ=255)
                                                                        cGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCg                                    <  4:117858/133‑1 (MQ=255)
                                                                        cGGGCCGGACTTGCACCCCCCACATCCCCCCCCCCATCCGACGTCGGGGGTTTTTTTTGCCCGGAATTTGCGGGGAGCGGGTTTTCTCGCGAATGCGTTGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                   <  3:159467/134‑1 (MQ=37)
                                                                          gTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCggg                                  >  5:62592/1‑133 (MQ=255)
                                                                              ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                            <  7:290372/135‑1 (MQ=255)
                                                                              ggACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                            <  7:255638/135‑1 (MQ=255)
                                                                                aCTTGGACCCCGAACTTCCCCCTCCCCAAACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCACATTCGCGGGTAAACCCGGGTCAGgcg                           >  6:22697/1‑134 (MQ=255)
                                                                                       ccccGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCa                    >  4:95154/1‑134 (MQ=255)
                                                                                                accccccccccATACTATGTCGGGGGTTTTTTTTGCCCGGATTTTGCGGGGAGGGGTTTTTCTCGCTATTGGTTTGTCTAATACCCCGCCGCATTCGGGGGTAAACCCGGGTCAGGCGAAAACCACTACACGAcc          <  8:90346/134‑1 (MQ=255)
                                                                                                acccccTCCCCATACGAAGTCGGGGGTTTTTTTTGCCGGTATTTTGCGGGGAGGGGGTTTTCTCGGGAATGCGTTGTCGAAAACACCGCCGCATTCGCGGGTAAACCCGGGTCTGGCGCAAACCACTAAACGAcc          <  8:303823/134‑1 (MQ=255)
                                                                                                  ccccTCCCCATCCGATGGCGGGGGTTTTTTTTGCGCGTTATTGGCGGGTTGGGGTTTTTCTCGCGATGCCGATGTCTAATCGCCGGCGGCATTCGCGGGTAAACCCGGGTCAGGCGAAAACCACTACACGAcc          <  4:85021/133‑1 (MQ=255)
                                                                                                    ccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAACCGCCCGGCt      >  8:26674/1‑135 (MQ=255)
                                                                                                        cccATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCTAATGCGATGTCGAATACCCCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTc  <  5:115040/135‑1 (MQ=255)
                                                                                                                          |                                                                                                                     
GAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTC  >  ACB122/120853‑121092

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