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,6210GT98.0% 154.8 / ‑6.3 49M845I (ATG→ATT) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (29/19);  total (29/20)
Fisher's exact test for biased strand distribution p-value = 4.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGT  >  NC_000913_3_pae_pgi/1731475‑1731764
                                                                                                                                                  |                                                                                                                                                
cATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCg                                                                                                                                                >  2:129544/1‑149 (MQ=255)
    gCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGc                                                                                                                                            >  2:674944/1‑149 (MQ=255)
            cTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACgg                                                                                                                                    >  2:389020/1‑149 (MQ=255)
                  caacaaCGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGa                                                                                                                              >  1:758341/1‑149 (MQ=255)
                   aacaacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGaa                                                                                                                             >  1:204031/1‑149 (MQ=255)
                   aacaacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCCCGAAACCCTACGCGAATGTCTGCAAAATTTTTATGATCTTCCCGCACTGGAACCTCTCATTCGTCGCCTGAACGCTGGCgca                                                                                                                             >  2:269260/1‑147 (MQ=255)
                     caacGGCTGCGGGCCAGTGAGTTGGTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCGGCAAGATGTTCCTGATCCTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAAt                                                                                                                           <  2:131125/149‑1 (MQ=255)
                     caacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGGGGCGAAAt                                                                                                                           >  1:676556/1‑149 (MQ=255)
                     caacGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTCCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAAt                                                                                                                           >  1:675983/1‑149 (MQ=255)
                         ggCTGCGCGCCAGTCAGGTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCaa                                                                                                                       <  1:469962/149‑1 (MQ=255)
                          gCTGCGCGCCAGGCAGTTGCTGGAAAGCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCaaa                                                                                                                      <  1:242418/149‑1 (MQ=255)
                            tGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTTGTCGCCTGAACGGTGGCGAGATTCAAata                                                                                                                    >  1:276192/1‑149 (MQ=255)
                                 gcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGa                                                                                                               >  1:632598/1‑149 (MQ=255)
                                 gcCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGa                                                                                                               <  1:389020/149‑1 (MQ=255)
                                         gtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGa                                                                                                       >  2:638174/1‑149 (MQ=255)
                                          ttGCTGGAAATCGCTCAGGGATCTCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGAc                                                                                                      <  1:608196/149‑1 (MQ=255)
                                          ttGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGAc                                                                                                      <  1:129544/149‑1 (MQ=255)
                                                              atatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCg                                                                                  >  2:645349/1‑149 (MQ=255)
                                                               tatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGc                                                                                 <  2:273815/149‑1 (MQ=255)
                                                               tatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGc                                                                                 <  1:567556/149‑1 (MQ=255)
                                                                 tCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGGCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCTTCGCCTTTGGCca                                                                               <  2:491166/149‑1 (MQ=255)
                                                                    ggATTTTCCGGTTATTCTCGAAACCCTACGCAACTGCATTCCAGGTTGTTATGCTCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTCGCGAAATTCAAATATCCGATGTAACGACCACTACGCCATCGCCTTTCGCCACaa                                                                            <  2:165212/149‑1 (MQ=255)
                                                                    ggATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCCAGCTGTTTATGATCTTCACGCACTGAAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACaa                                                                            >  1:227317/1‑149 (MQ=255)
                                                                       ttttCCGGTCATTCTCGAAACACTACGCGAATGTCTGCAAAATGGTTAGGATCTTCCCGAACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACACCTACGCCAACACCTTGCGCcaacattt                                                                         >  1:157134/1‑143 (MQ=255)
                                                                        tttCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGttt                                                                        <  2:204031/149‑1 (MQ=255)
                                                                              gTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTaa                                                                                                          >  1:271946/1‑109 (MQ=255)
                                                                              gTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTaa                                                                                                          <  2:271946/109‑1 (MQ=255)
                                                                              gTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGtt                                                                  <  2:632598/149‑1 (MQ=255)
                                                                                             cTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTAGCGTGGCCTCCGGAACGTGTGCGAAATGCAAATATCCGATGTACCGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGa                                                  >  2:289414/1‑149 (MQ=255)
                                                                                                   gAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTAt                                            <  1:191720/149‑1 (MQ=255)
                                                                                                   gAATGTCGGCAAGATGTTTATGTTCTTCCCGCACTGGAACGTTGGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACCACCACTACTCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTAt                                            <  2:377772/149‑1 (MQ=255)
                                                                                                          tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCag                                     >  2:256430/1‑149 (MQ=255)
                                                                                                          tGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTAAGTCCCag                                     >  1:572406/1‑149 (MQ=255)
                                                                                                             aaGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATGCAAATATCCGATGTAACGACCACTACGACCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACAAGAGc                                  >  2:686637/1‑149 (MQ=255)
                                                                                                                    ttATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TCTGTCGCGGAATTTATGTACCCGAGCGACGccc                           >  2:147525/1‑149 (MQ=255)
                                                                                                                      aTGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAGTATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCg                         >  2:467321/1‑149 (MQ=255)
                                                                                                                         aTCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTCTCTCCCAGAGCGACGCCCCGCtc                      >  2:198378/1‑148 (MQ=255)
                                                                                                                         aTCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TACGGGGGGGAATTTATGTACCAGAGCGACGCCCCGCTg                      <  2:758341/149‑1 (MQ=255)
                                                                                                                          tCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTgg                     >  2:66730/1‑149 (MQ=255)
                                                                                                                          tCTTCCCGCACTGGAACGTTGGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGTTGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCATAAG‑GGCGGAATTTATGTACCAGAGCGACGCCCCGCTgg                     <  2:723801/149‑1 (MQ=255)
                                                                                                                                cGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAgc               <  1:445016/149‑1 (MQ=255)
                                                                                                                                cGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAgc               >  2:642066/1‑149 (MQ=255)
                                                                                                                                cGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAgc               >  1:228454/1‑149 (MQ=255)
                                                                                                                                     tGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑CATGTCGCGGAATTTATGTACCAGTGCGACGCCCCGCCGGCAGAGCGCCgg          <  1:645349/149‑1 (MQ=255)
                                                                                                                                       gAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGa                                >  2:57506/1‑125 (MQ=255)
                                                                                                                                       gAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGa                                <  1:57506/125‑1 (MQ=255)
                                                                                                                                           gTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATcc    >  1:591393/1‑149 (MQ=255)
                                                                                                                                            tttGATTCGTCGCCTAATCGGTGGCAAAATTCACATATCCGATATAACGACCACTACGCCCTCGCCTATGGCCACAAGTTTATTGTTCGGC‑TACGTCGCGGAATTTATGTACCATAGCGACCCCCCGCAGGCATAGCGCCCGACATCCg   >  2:646647/1‑149 (MQ=255)
                                                                                                                                             ttGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCACGCTGGCAGAGCGCCGGGCATCCGt  >  1:772410/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGC‑TATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGT  >  NC_000913_3_pae_pgi/1731475‑1731764

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