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,6210GT95.8% 74.3 / NA 24M845I (ATG→ATTlhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (17/6);  minor base C (1/0);  total (18/6)
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

CCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATC  >  NC_000913_3_pae_pgi/1731484‑1731761
                                                                                                                                         |                                                                                                                                            
ccGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCAGGTCATTCTCGAAACCCTACGCGAATCTCTGCAAGATGTTTATGACCTTCCCGCACTGCAACGTTCGATCCGTCGCCTgca                                                                                                                                   >  2:100654/1‑147 (MQ=255)
     ttGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTg                                                                                                                              <  2:95637/149‑1 (MQ=255)
        gCAACAACGGCTGCGTGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACggccgcc                                                                                                                           >  2:403309/1‑144 (MQ=255)
         caacaaCGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGa                                                                                                                          >  1:365695/1‑149 (MQ=255)
          aacaacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGAACTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCCACTGAACGGTGGCGaa                                                                                                                         >  2:292206/1‑149 (MQ=255)
             aacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAAtt                                                                                                                      <  1:58701/149‑1 (MQ=255)
              acGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTc                                                                                                                     >  1:366256/1‑149 (MQ=255)
                 gCTGGGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGGAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCaaa                                                                                                                  <  2:266966/149‑1 (MQ=255)
                          cagtcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATg                                                                                                         >  1:47181/1‑149 (MQ=255)
                                      ggAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCACATATCCGATGTAACGACCGCTa                                                                                             >  2:307471/1‑149 (MQ=255)
                                        aaaTCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACg                                                                                           >  2:463806/1‑149 (MQ=255)
                                            cGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTc                                                                                                                                          >  1:305797/1‑98 (MQ=255)
                                            cGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTc                                                                                                                                          <  2:305797/98‑1 (MQ=255)
                                                           ggATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                                                                        <  1:395876/149‑1 (MQ=255)
                                                                    ggTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGt                                                               >  1:278744/1‑149 (MQ=255)
                                                                        aTTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGCTTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAGGTTTGTTGTTCgg                                                           >  2:273500/1‑149 (MQ=255)
                                                                        aTTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCgg                                                           >  1:71875/1‑149 (MQ=255)
                                                                                    cTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTCTGCCGCGGa                                               >  2:166110/1‑149 (MQ=255)
                                                                                                 tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCag                                  >  1:374687/1‑149 (MQ=255)
                                                                                                                      ccGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCagag             >  1:242540/1‑149 (MQ=255)
                                                                                                                          aCTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCCGGGCGcc         >  2:279771/1‑149 (MQ=255)
                                                                                                                            tGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCgg       <  1:463806/149‑1 (MQ=255)
                                                                                                                             ggAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATATATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCggg      >  1:370408/1‑149 (MQ=255)
                                                                                                                                 cGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATc  >  1:148485/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                            
CCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATC  >  NC_000913_3_pae_pgi/1731484‑1731761

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: