Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,573,300 T→G D69A (GAT→GCT)  eutD ← phosphate acetyltransferase EutD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,573,3000TG97.5% 141.3 / ‑3.5 40D69A (GAT→GCT) eutDphosphate acetyltransferase EutD
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (20/19);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01

TACCGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGATCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGCTTTCAGCACACG  >  NC_000913/2573163‑2573430
                                                                                                                                         |                                                                                                                                  
tACCGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGCTGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGagca                                                                                                                                  >  2:796311/1‑140 (MQ=255)
tACCGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGagca                                                                                                                                  >  2:960634/1‑140 (MQ=255)
 aCCGGCGCTGACCATTGCTGCGGCGAAC‑TCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGagcaat                                                                                                                                <  1:213593/140‑1 (MQ=255)
  ccGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGatcaat                                                                                                                                >  2:797703/1‑140 (MQ=255)
  ccGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGagcaat                                                                                                                                <  2:926520/140‑1 (MQ=255)
  ccGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGagcaat                                                                                                                                <  2:496343/140‑1 (MQ=255)
  ccGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGagcaat                                                                                                                                <  2:1175058/140‑1 (MQ=255)
     gcgcTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCAc                                                                                                                             >  2:975475/1‑140 (MQ=255)
     gcgcTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCAc                                                                                                                             >  2:935200/1‑140 (MQ=255)
       gcTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGATCAATCACCt                                                                                                                           <  1:797703/140‑1 (MQ=255)
        cTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCACATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTg                                                                                                                          >  1:61987/1‑140 (MQ=255)
            ccATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGc                                                                                                                      <  1:333597/140‑1 (MQ=255)
                tgctgcGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGt                                                                                                                  >  1:1035648/1‑140 (MQ=255)
                     cggcgAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATc                                                                                                             >  2:867423/1‑140 (MQ=255)
                     cggcgAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATc                                                                                                             >  2:1020048/1‑140 (MQ=255)
                       gcgAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGc                                                                                                           <  1:712531/140‑1 (MQ=255)
                       gcgAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGc                                                                                                           <  2:1175550/140‑1 (MQ=255)
                          aaCATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCAc                                                                                                        >  1:1152032/1‑140 (MQ=255)
                          aaCATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCAc                                                                                                        >  2:183462/1‑140 (MQ=255)
                                    ggTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTg                                                                                              >  1:123598/1‑140 (MQ=255)
                                      tCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAg                                                                                            >  1:88038/1‑140 (MQ=255)
                                         gTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAgcgc                                                                                         >  1:334979/1‑140 (MQ=255)
                                         gTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAgcgc                                                                                         >  2:811615/1‑140 (MQ=255)
                                         gTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAgcgc                                                                                         >  2:307364/1‑140 (MQ=255)
                                                      gcgcATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGt                                                                            >  1:162796/1‑140 (MQ=255)
                                                       cgcATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGtt                                                                           >  2:1105730/1‑140 (MQ=255)
                                                               gcggcgTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGTTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGa                                                                   <  1:329305/140‑1 (MQ=255)
                                                                   cgTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTgg                                                               <  1:906074/140‑1 (MQ=255)
                                                                                gcgcgGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATggg                                                  <  1:1105730/140‑1 (MQ=255)
                                                                                 cgcgGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGc                                                 <  1:960634/140‑1 (MQ=255)
                                                                                 cgcgGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGc                                                 <  1:796311/140‑1 (MQ=255)
                                                                                       ccagccagCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGcc                                           <  2:88038/140‑1 (MQ=255)
                                                                                               cGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCttg                                   <  2:1035648/140‑1 (MQ=255)
                                                                                               cGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCttg                                   <  1:867423/140‑1 (MQ=255)
                                                                                                gATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCttgt                                  <  2:876671/140‑1 (MQ=255)
                                                                                                 aTGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCttgtt                                 <  2:1152032/140‑1 (MQ=255)
                                                                                                 aTGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCttgtt                                 <  1:183462/140‑1 (MQ=255)
                                                                                                 aTGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCttgtt                                 <  1:1193964/140‑1 (MQ=255)
                                                                                                     gCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATg                             >  1:427102/1‑140 (MQ=255)
                                                                                                             ttcCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATAtt                     >  2:697344/1‑140 (MQ=255)
                                                                                                              tcCCGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTg                    >  1:313628/1‑140 (MQ=255)
                                                                                                                ccGCATTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTgcg                  <  1:1006254/140‑1 (MQ=255)
                                                                                                                    aTTGCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGc              >  2:770384/1‑140 (MQ=255)
                                                                                                                       gCGAGGTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGCttt           <  2:334979/140‑1 (MQ=255)
                                                                                                                            gTTGCCATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGCTTTCAGca      <  1:811615/140‑1 (MQ=255)
                                                                                                                                ccATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGCTTTCAGCACACg  >  2:415840/1‑140 (MQ=255)
                                                                                                                                ccATGCGGAGCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGCTTTCAGCACACg  <  1:1184588/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                  
TACCGGCGCTGACCATTGCTGCGGCGAACATCAGAGGGTCGGTAAGTTTTTCCAGCGCATCCGGCGGCGTTTTTTCGCCCGCGCGGGCCAGCCAGCGATGAGCAAATTCTTCCCGCATTGCGAGGTTGCCATGCGGATCAATCACCTGTAGCCCGTCCATCGCCACGCCGTGACTGAGCGCAAACTGACGAAGTTCAAACGGATTGGCGACCAGAATGGGCGTTGCCAGACCTTGTTGATGTAAATATTGCGCGGCTTTCAGCACACG  >  NC_000913/2573163‑2573430

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