Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_pgi 1,731,621 G→T M845I (ATG→ATT)  lhr → putative ATP‑dependent helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_pgi1,731,6210GT97.4% 119.0 / NA 38M845I (ATG→ATT) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (19/18);  minor base C (1/0);  total (20/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 = 1.00e+00

CCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCG  >  NC_000913_3_pae_pgi/1731484‑1731763
                                                                                                                                         |                                                                                                                                              
ccGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGaa                                                                                                                                     >  2:436885/1‑149 (MQ=255)
 cGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAAc                                                                                                                                    <  2:679875/149‑1 (MQ=255)
     ttGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTg                                                                                                                                <  2:729923/149‑1 (MQ=255)
     ttGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTg                                                                                                                                <  2:579152/149‑1 (MQ=255)
        gCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCg                                                                                                                             >  1:410311/1‑149 (MQ=255)
           acaacGGCTGCGGGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaaa                                                                                                                          <  1:125902/149‑1 (MQ=255)
            caacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAAt                                                                                                                         <  2:600607/149‑1 (MQ=255)
                ggcggcGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCaa                                                                                                                     <  1:903773/149‑1 (MQ=255)
                ggCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCCGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATATTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCaa                                                                                                                     <  2:737595/149‑1 (MQ=255)
                ggCTGCGCGCCAGCCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCaa                                                                                                                     <  1:3106/149‑1 (MQ=255)
                          cagtcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATCCGTCGCCTGGACGGTGGCGAACTTCAAATATCCGacg                                                                                                           >  2:478318/1‑147 (MQ=255)
                             tcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTat                                                                                                        >  2:480777/1‑148 (MQ=255)
                                      ggAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTa                                                                                               <  1:462614/149‑1 (MQ=255)
                                               tCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGTTCTTGCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGc                                                                                      <  1:548494/149‑1 (MQ=255)
                                                            gATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGGTGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTACAGACCCCTACGCCCCCGCAGTACGCCACAAg                                                                         >  2:601933/1‑149 (MQ=255)
                                                                ttCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTa                                                                     <  1:544374/149‑1 (MQ=255)
                                                                 tCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGGACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTAt                                                                    <  2:624915/149‑1 (MQ=255)
                                                                          tctcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCt                                                           >  1:441715/1‑149 (MQ=255)
                                                                             ccaaaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTTTGTCCTTACCTTTCGCCACAAGTTTATTGTTCGGCTATg                                                        <  1:480777/147‑1 (MQ=255)
                                                                               aaaCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTc                                                      >  1:886931/1‑149 (MQ=255)
                                                                                 aCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCACACGCCTTTAGCCACGAGGGTATTGTTCGGCTATGTcgc                                                    >  2:49414/1‑149 (MQ=255)
                                                                                   ccTACGCGTATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGACACAAGTTTATTGTTCGGCTATGTCGCgg                                                  >  2:414917/1‑149 (MQ=255)
                                                                                     tatgccaaTGCCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTGGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGaa                                                <  1:421411/143‑1 (MQ=255)
                                                                                          gAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAGTTCAAATATCCGATGTAACGACAACACGGGACTTGGCTTTCGCCACAAGTCTATTGTTAGGATATGTCGCGGAATTAAt                                           >  1:376041/1‑149 (MQ=255)
                                                                                          gAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACACCCTCGCCTTTCGCCACAAGTTTATAGTTCGGCTATGTCGCGGAATTTAt                                           >  1:793221/1‑149 (MQ=255)
                                                                                           aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATg                                          <  2:726336/149‑1 (MQ=255)
                                                                                           aaTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATg                                          >  1:165324/1‑149 (MQ=255)
                                                                                              gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTAc                                       <  2:268645/149‑1 (MQ=255)
                                                                                              gTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTAc                                       <  1:4044/149‑1 (MQ=255)
                                                                                                  gCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCCATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCaga                                   <  2:410311/149‑1 (MQ=255)
                                                                                                   cAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCagag                                  >  1:404062/1‑149 (MQ=255)
                                                                                                         gTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTCAACGGTGGCGAAATTCAAATATCAGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGc                            >  2:919184/1‑149 (MQ=255)
                                                                                                            tATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTGACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGCATTTATGTCACAGAGCGACGcccc                         >  1:810209/1‑149 (MQ=255)
                                                                                                               gATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGAGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGCCTCTGTCGCGGACTTGATGTACCAGAGCCACGccccccc                      >  2:927413/1‑146 (MQ=255)
                                                                                                                aTCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTg                     >  2:481143/1‑149 (MQ=255)
                                                                                                                                  gTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATcc   >  1:685611/1‑149 (MQ=255)
                                                                                                                                   tttGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCg  >  1:956889/1‑149 (MQ=255)
                                                                                                                                   tttGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCg  >  1:703506/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCG  >  NC_000913_3_pae_pgi/1731484‑1731763

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