Predicted mutation
evidence seq id position mutation annotation gene description
RA ACB122 120,983 (T)9→8 intergenic (+280/+92) 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.95.5% 179.5 / 2.8 45intergenic (+272/+100)polA/engBDNA polymerase I/putative GTP‑binding protein EngB
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (28/14);  total (28/17)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTC  >  ACB122/120850‑121110
                                                                                                                             |                                                                                                                                       
aCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgcc                                                                                                                               <  8:300990/136‑4 (MQ=255)
aCTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttgc                                                                                                                                <  8:1717/135‑3 (MQ=255)
    aaTGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGttttttttggacgg                                                                                                                            >  3:199496/1‑129 (MQ=255)
               aaGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAACTCCGGGGGGTTTTTTGCCCGGCATTTGCCgg                                                                                                                 >  7:231666/1‑134 (MQ=255)
               aaGCCCTTCGATGCTCCTATGTTGAGTGAAGTGTTGGCAGATATCTGGACCCCGCCCTATCGGTCCGGACTTGGACCCCGAACTTCCCACTCCCCATACGAAGTCCGGGGGTTTTTTTGCCCGGAATTTGCCgg                                                                                                                 >  7:58789/1‑134 (MQ=255)
                       cGATGCTCCTATGTTGGGTTAAGTGTTGGCAGATATCGGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGt                                                                                                         <  7:275762/134‑1 (MQ=255)
                        gaagCTCCTATCTTCTGTTAAGTCTATCCAGAACTCTGGACCCCGCCCTAGCGGTCCGGACGTGGAGCCCGCACTTCCGCCTCCCCCTACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGt                                                                                                         <  2:155413/130‑1 (MQ=255)
                         aTGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGtttt                                                                                                      >  1:49610/1‑135 (MQ=255)
                           gCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTtct                                                                                                    <  1:223357/135‑1 (MQ=255)
                               cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGATCGGGTTTTCTcgct                                                                                                >  7:245709/1‑134 (MQ=255)
                               cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTcgagt                                                                                               >  5:330974/1‑133 (MQ=255)
                               cTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGa                                                                                               <  4:212493/136‑1 (MQ=255)
                                  tgttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             <  8:81395/135‑1 (MQ=255)
                                    ttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                           <  7:201951/135‑1 (MQ=255)
                                    ttgtGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACCTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                           <  3:170686/135‑1 (MQ=255)
                                     tgtgTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGc                                                                                           <  1:189300/134‑1 (MQ=255)
                                         ttAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATAAGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                     >  5:145854/1‑136 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             >  4:314396/1‑127 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             <  3:314396/127‑1 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGt                                                                                      >  2:97258/1‑134 (MQ=255)
                                          tAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCACGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTc                                                                                     >  4:269660/1‑135 (MQ=255)
                                            aGTGTTGGCAGATATCTGGCCCCCGCCCTAGCGGTCCGGACTTGGACCCCGCACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                   >  1:183094/1‑135 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                                  >  8:203894/1‑136 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGaa                                                                                  >  3:291366/1‑136 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGa                                                                                   >  1:298023/1‑135 (MQ=255)
                                            aGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCACCCCAAACGAAGTACGGGGTTTTTTTTGCCCGGAATTAGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCg                                                                                    >  2:101819/1‑134 (MQ=255)
                                                    cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACAccgc                                                                          >  5:158556/1‑136 (MQ=255)
                                                    cAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTGTTTTTTGCCCGGAATTTGCCGGGAGCGGGTGTTCTCGCGATTGCGCTGTCGAATACCccgc                                                                          >  7:126810/1‑136 (MQ=255)
                                                      gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             <  6:311700/115‑1 (MQ=255)
                                                      gATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAAt                                                                                             >  5:311700/1‑115 (MQ=255)
                                                           cTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGATGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATCTCGAATACACCGCCGCATTc                                                                   >  7:221552/1‑136 (MQ=255)
                                                             ggACCCCCCCCTCGCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACTAAGTCCGGGTTTTTTTTTGCCCGGAATTTGCCGGTAGCGGGTTTTCTCGCGAATGCGATGTCTAAAACACCGCCGCATTcgc                                                                 <  4:212791/136‑1 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGCGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGGAAAcc                                                        >  5:229695/1‑136 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTaagc                                                         >  1:231617/1‑133 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAcc                                                        >  7:185674/1‑136 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCAATCGCGGGTAAAcc                                                        >  5:301418/1‑136 (MQ=255)
                                                                      ccTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGG‑‑TTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAAccc                                                       >  5:141749/1‑136 (MQ=255)
                                                                          gCGGGCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCgg                                                     <  3:304719/135‑1 (MQ=255)
                                                                                 ggACTTGGACCCCGAACTTCGCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                              <  2:23491/135‑1 (MQ=255)
                                                                                 ggACTTGGACCCCGAACTTCCCCCTCCCCATCCGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGgc                                              <  1:24446/135‑1 (MQ=255)
                                                                                   aCTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCa                                           <  6:158556/136‑1 (MQ=255)
                                                                                               aaCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCTATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAACCCACTAAAc                                >  8:116435/1‑135 (MQ=255)
                                                                                                aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGa                              >  6:232575/1‑136 (MQ=255)
                                                                                                aCTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGa                              >  4:292168/1‑136 (MQ=255)
                                                                                                 cTTCCCCCTCCCCATACGAAGTCGGGGGTTTTTTTTGCCGGTATTTTGCCGGGAGGGGGTTTTCTCGCTAATGCGATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACAACTAAACGa                              <  2:231617/135‑1 (MQ=255)
                                                                                                       ccTCCCCATACGAAGTCCGGGGTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCACACCACTAAACGACCGGCtt                       >  3:93587/1‑136 (MQ=255)
                                                                                                            ccATACTACGTCGGGGTTTTTTTTTCCCCGTATTTTGCGGGGAGCGGTTTTGCCCGCGATTGGGATGTCTAATACACCGCCGCATTCGGGTTTAACCCCGGGTCAGGCGCAAACCACAAAACGACCGGCTTGTCtt                  <  8:125024/136‑1 (MQ=255)
                                                                                                            ccATACGAAGTCCGGGGTTTTTTTTGCACGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCCTTCGCGGGTAACCCCGGGTCAGGCGCAAACCACTAAACGCCCGGCTTGTCtg                  >  6:269587/1‑135 (MQ=255)
                                                                                                                        ccggggTTTTTTTTGCCGGTAATTTGCCGGGAGCGGGTTTTCTCGCAAATGCTATGTCTAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc       <  7:203894/134‑1 (MQ=255)
                                                                                                                        ccggggTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTtccatc       >  8:259171/3‑136 (MQ=255)
                                                                                                                             tttttttttGCCGGGATTTTGCGGGGAGTGGGTTTTCTGGCGAATGCTATGTCGATTCCCCGCCCGAATTGGCGGTTAAACGGGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTc  <  8:317904/136‑1 (MQ=255)
                                                                                                                             |                                                                                                                                       
ACTGAATGGCTGGTGAAGCCCTTCGATGCTCCTATGTTGTGTTAAGTGTTGGCAGATATCTGGACCCCGCCCTAGCGGTCCGGACTTGGACCCCGAACTTCCCCCTCCCCATACGAAGTCCGGGGTTTTTTTTTGCCCGGAATTTGCCGGGAGCGGGTTTTCTCGCGAATGCGATGTCGAATACACCGCCGCATTCGCGGGTAAACCCGGGTCAGGCGCAAACCACTAAACGACCGGCTTGTCTTCCAGTTCCATCCAGTC  >  ACB122/120850‑121110

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: