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.2% 178.6 / NA 54M845I (ATG→ATT) lhrputative ATP‑dependent helicase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (28/25);  minor base C (0/1);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTAC  >  NC_000913_3_pae_pgi/1731475‑1731766
                                                                                                                                                  |                                                                                                                                                 
cagcGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCg                                                                                                                                                 <  2:874918/146‑1 (MQ=255)
 aTCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAACTCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGt                                                                                                                                                <  2:576572/149‑1 (MQ=255)
              tgggCAACAACGGCTGCTCGCCTGGCATTTGCTGGAAATCGCTCAGGGATATCCGGCTTTGCCGGTCATTGTAGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTTATTCGTCGCCCGAACGGTg                                                                                                                                   <  2:891275/147‑1 (MQ=255)
                                   cagtcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATg                                                                                                              <  1:885621/149‑1 (MQ=255)
                                   cagtcagtTGCTGGAAAGCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGCCGTCCGAACGGTGGCGAAATTCAAATATCCGATg                                                                                                              <  1:257405/149‑1 (MQ=255)
                                      tcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGCACGGTGGCGAAATTCAAATATCCGGTGTaa                                                                                                           >  1:921250/1‑149 (MQ=255)
                                      tcagtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTaa                                                                                                           >  2:275666/1‑149 (MQ=255)
                                         gtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAAcca                                                                                                        >  1:224290/1‑147 (MQ=255)
                                         gtTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGa                                                                                                        >  1:187149/1‑149 (MQ=255)
                                             cTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGCCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCAc                                                                                                    >  1:799464/1‑149 (MQ=255)
                                                  aaTCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGc                                                                                               >  1:508667/1‑149 (MQ=255)
                                                      gCTCAGGGATATCCGGATTTCCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTc                                                                                           <  1:759683/149‑1 (MQ=255)
                                                        ttatggATATCCGGATTTTCCGGTCATTCTCGAACCCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGc                                                                                         <  2:1047358/145‑1 (MQ=255)
                                                          aGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCt                                                                                                                                           >  2:759805/1‑97 (MQ=255)
                                                          aGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCt                                                                                                                                           <  1:759805/97‑1 (MQ=255)
                                                          aGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCt                                                                                       >  2:685836/1‑149 (MQ=255)
                                                           gggATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCtt                                                                                      <  1:1021543/149‑1 (MQ=255)
                                                             gattgcaggATTTACTCGTCAGTCTCGTAACCCTACGCGAATGTCTGGAAGATGTTTATGATCTTCCCGCACTGGCACGTTTGAGTGGTCGCCTGAATGGTTTCGCAATTCAAATATCCGATGTAACTACCACTACGCCCTCGCCTTTc                                                                                    <  2:368420/142‑1 (MQ=255)
                                                              atatCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCACGCCTTTCg                                                                                   >  2:262652/1‑149 (MQ=255)
                                                                            cGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGc                                                                                  <  1:889934/136‑1 (MQ=255)
                                                                            cGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGc                                                                                  >  2:889934/1‑136 (MQ=255)
                                                                             ggTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGt                                                                    >  1:789671/1‑149 (MQ=255)
                                                                                cATTCTCGAAACCCTACGCGAATGTCTGCCAGCTGTTTATGCTCTTCCCGCACTGGATCGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAGTATCCGATGTAACGACCGCTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCg                                                                 <  2:625398/149‑1 (MQ=255)
                                                                                 aTTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCgg                                                                >  2:820186/1‑149 (MQ=255)
                                                                                     tcGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTAt                                                            <  2:129798/149‑1 (MQ=255)
                                                                                              tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAGTTCAAATATCCGATGGAACGACCACTACGCCCTCGCCTTTCGCCACAAGTGTATTGTTCGTCTATGTTGGGGaa                                                   >  1:237490/1‑149 (MQ=255)
                                                                                              tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCtt                                                                                      <  1:527900/114‑1 (MQ=255)
                                                                                              tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCtt                                                                                      >  2:527900/1‑114 (MQ=255)
                                                                                              tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCtt                                                                                      <  1:323751/114‑1 (MQ=255)
                                                                                              tACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCtt                                                                                      >  2:323751/1‑114 (MQ=255)
                                                                                                   gAATGTCTGTAAGATGTTTATGCTCTTCCCGCACTGGAACTTTTGATTCGTCGTCTGAACGGTGGCGACATTCAAATATCCGATGTAACGACCCCTACGCCCTCGCCTTTCGCCACCAGTTTTTTGTTCGGCTATGTCGCGGAAtttcg                                              <  2:398967/149‑3 (MQ=255)
                                                                                                   gAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGACTTTTGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTAt                                              <  1:523004/149‑1 (MQ=255)
                                                                                                     aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTAt                                              >  1:312868/1‑147 (MQ=255)
                                                                                                     aTGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTAt                                              <  2:312868/147‑1 (MQ=255)
                                                                                                      tttCTGGAAGATGTTTATGATCTGCCCGCGCTGTAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTa                                           <  1:1025044/147‑1 (MQ=255)
                                                                                                           tcAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCaga                                      <  1:429958/148‑1 (MQ=255)
                                                                                                             actaTGTTTATGATCTCGCCGCATTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGc                                    <  2:425913/146‑1 (MQ=255)
                                                                                                             aaGATGTTTATGATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTAGGTCGCGGACTTTCTGTATCcgaga                                    >  2:429958/1‑148 (MQ=255)
                                                                                                                        gATCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTAGAATGCCACCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCt                         >  1:515977/1‑149 (MQ=255)
                                                                                                                         aTCTTCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCCCTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTg                        >  2:504863/1‑149 (MQ=255)
                                                                                                                            ttCCCGCACTGGAACGTTTGCTTCGTCGCTCGCGCGGTACGTAAACTCAAATCTCCGATGTAACGACCAGTACGCCCACGCCTTTCGCCACAAGTTTATTGTCCGGCTGTGTCGCGAAAGTTATGTACCAGAGCGACGCGCCGCTggac                     >  2:181130/1‑147 (MQ=255)
                                                                                                                            ttCCCGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCa                     <  1:275666/149‑1 (MQ=255)
                                                                                                                              cccGCACTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTCTTCGGCTATGTCGCGGAATTTATCTACCAGAGCGACGCCCCGCTGGCaga                   <  1:866677/149‑1 (MQ=255)
                                                                                                                               ccGCATTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc                                                                  >  2:549656/1‑101 (MQ=255)
                                                                                                                               ccGCATTGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTc                                                                  <  1:549656/101‑1 (MQ=255)
                                                                                                                                     tGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCgg            >  2:994379/1‑149 (MQ=255)
                                                                                                                                     tGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCgg            <  1:942572/149‑1 (MQ=255)
                                                                                                                                     tGGAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCgg            >  1:719634/1‑149 (MQ=255)
                                                                                                                                     tGGAACGATGAATCCGCCGCCCAAACGATCACGAAACCCAAATATCCGATGTAACGACCACTACGCCCTCGCCAATCGCCACAAGTCAATTCTCCGGCTATGTCGCCGAATTTATGTACCAGACCGACGCCCCGCTGGCAGAGCCCCgc            <  2:203408/149‑2 (MQ=255)
                                                                                                                                       gAACGTTTGATTCTTCGCCAGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGc          >  2:888184/1‑149 (MQ=255)
                                                                                                                                       gAACGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATATACCAGAGCGCCGACCCGCTGGCAGAGGGCgggtc          >  1:287323/1‑147 (MQ=255)
                                                                                                                                          cGTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGGTCGGCTATGTCGCGGAACTTATGTACCAGAGCAACGCCCCGCTGGCAGAGCGCCGGGCATc       >  1:29543/1‑149 (MQ=255)
                                                                                                                                           gTTTGATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATcc      >  2:280561/1‑149 (MQ=255)
                                                                                                                                               gATTCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGCATTTATGTACCAGAGCGAAGCCCCGCTGGCAGAGCGCCGGGCATCCGTAc  >  2:606395/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CATCGCACCCCGCTTTGGCAACAACGGCTGCGCGCCAGTCAGTTGCTGGAAATCGCTCAGGGATATCCGGATTTTCCGGTCATTCTCGAAACCCTACGCGAATGTCTGCAAGATGTTTATGATCTTCCCGCACTGGAACGTTTGATGCGTCGCCTGAACGGTGGCGAAATTCAAATATCCGATGTAACGACCACTACGCCCTCGCCTTTCGCCACAAGTTTATTGTTCGGCTATGTCGCGGAATTTATGTACCAGAGCGACGCCCCGCTGGCAGAGCGCCGGGCATCCGTAC  >  NC_000913_3_pae_pgi/1731475‑1731766

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: ‑