Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 232,840 T→A V96E (GTG→GAG)  dnaQ → DNA polymerase III subunit epsilon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273232,8400TA95.1% 147.0 / ‑0.3 41V96E (GTG→GAG) dnaQDNA polymerase III subunit epsilon
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (22/17);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.66e-01

GTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAA  >  NZ_CP009273/232753‑232918
                                                                                       |                                                                              
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  2:940255/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  2:203364/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              >  2:713059/1‑90 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  2:1161847/90‑1 (MQ=255)
gTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGa                                                                              <  1:558168/90‑1 (MQ=255)
   ttGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  1:144707/87‑1 (MQ=255)
   ttGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              >  2:144707/1‑87 (MQ=255)
     gCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCAt                                                                         >  2:649014/1‑90 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  2:205636/82‑1 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              >  1:205636/1‑82 (MQ=255)
        gATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAAc                                                                      >  2:40288/1‑90 (MQ=255)
          tGACTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGACGATGAGTTCATGGACTATGTTCGCGGGCCGGAGGTGGAGATCAATAACGc                                                                    >  1:309970/1‑90 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  2:890921/80‑1 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              >  1:890921/1‑80 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGc                                                                    <  1:415444/90‑1 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGc                                                                    <  2:579073/90‑1 (MQ=255)
          tGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGc                                                                    <  1:469081/90‑1 (MQ=255)
            aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAg                                                                  >  2:1339463/1‑90 (MQ=255)
            aaTTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAg                                                                  >  2:804846/1‑90 (MQ=255)
                 ttGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTc                                                             <  2:149445/90‑1 (MQ=255)
                     tCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGata                                                         <  1:773676/90‑1 (MQ=255)
                       gATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATc                                                       <  1:649014/90‑1 (MQ=255)
                        aTAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCg                                                      >  2:391982/1‑90 (MQ=255)
                               cACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTAt                                               <  2:731934/90‑1 (MQ=255)
                                  gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              <  1:1370508/56‑1 (MQ=255)
                                  gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGgaga                                                                              >  2:1370508/1‑56 (MQ=255)
                                  gTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGa                                            >  1:584862/1‑90 (MQ=255)
                                      gCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTAc                                        >  1:213819/1‑90 (MQ=255)
                                           aGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGtt                                   >  2:1053953/1‑90 (MQ=255)
                                                ccGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGt                              <  2:653272/90‑1 (MQ=255)
                                                ccGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGt                              <  2:213819/90‑1 (MQ=255)
                                                ccGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGt                              <  1:586330/90‑1 (MQ=255)
                                                   aTGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGc                           >  2:1346904/1‑90 (MQ=255)
                                                        ttCATGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTAAAg                      <  2:474831/90‑1 (MQ=255)
                                                           aTGGACTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAgcgc                   <  1:674848/90‑1 (MQ=255)
                                                               aCTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGata               >  1:553505/1‑90 (MQ=255)
                                                                cTATATTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGatat              >  2:694737/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATAtt             >  1:241209/1‑90 (MQ=255)
                                                                 tataTTCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTAAGGAGTACGAGTTTTCGTTGTTTAAGCGCGATAtt             >  1:773000/1‑90 (MQ=255)
                                                                     ttCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa         <  1:1339463/90‑1 (MQ=255)
                                                                     ttCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa         <  1:804846/90‑1 (MQ=255)
                                                                     ttCGCGGCGCGGAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGa         <  1:382415/90‑1 (MQ=255)
                                                                       cgcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg       <  2:1037749/90‑1 (MQ=255)
                                                                       cgcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAg       <  1:694737/90‑1 (MQ=255)
                                                                        gcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa      >  2:1318380/1‑90 (MQ=255)
                                                                        gcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa      >  2:1287447/1‑90 (MQ=255)
                                                                        gcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa      >  2:786668/1‑90 (MQ=255)
                                                                        gcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa      >  1:1573908/1‑90 (MQ=255)
                                                                        gcggcgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGa      >  1:1000854/1‑90 (MQ=255)
                                                                            cgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa  >  2:471824/1‑90 (MQ=255)
                                                                            cgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa  >  1:656306/1‑90 (MQ=255)
                                                                            cgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa  >  1:655997/1‑90 (MQ=255)
                                                                            cgcggAGTTGGAGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCaa  >  1:35407/1‑90 (MQ=255)
                                                                                       |                                                                              
GTATTGCCGATGAATTTTTGCTCGATAAGCCCACGTTTGCCGAAGTAGCCGATGAGTTCATGGACTATATTCGCGGCGCGGAGTTGGTGATCCATAACGCAGCGTTCGATATCGGCTTTATGGACTACGAGTTTTCGTTGCTTAAGCGCGATATTCCGAAGACCAA  >  NZ_CP009273/232753‑232918

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: