Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 460,036 C→T G383G (GGC→GGT lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913460,0360CT98.2% 192.6 / NA 55G383G (GGC→GGTlonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/24);  minor base G (0/1);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACA  >  NC_000913/459900‑460181
                                                                                                                                        |                                                                                                                                                 
ttGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGtgatg                                                                                                                                      >  1:258079/1‑150 (MQ=255)
 tGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGtga                                                                                                                                        <  1:487869/147‑1 (MQ=255)
      ttCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCgg                                                                                                                                >  1:405289/1‑150 (MQ=255)
         aaaGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAAt                                                                                                                            >  1:370198/1‑151 (MQ=255)
           aGCCGTGTCAACCAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATcc                                                                                                                          >  2:450907/1‑151 (MQ=255)
              cGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTg                                                                                                                       >  1:573974/1‑151 (MQ=255)
                   caacaaAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCa                                                                                                                   <  2:419696/150‑1 (MQ=255)
                   caacaaAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCAc                                                                                                                  >  2:652820/1‑151 (MQ=255)
                    aacaaAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCAcc                                                                                                                 <  1:69095/151‑1 (MQ=255)
                         aaTCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCAccgccg                                                                                                             >  2:337657/1‑150 (MQ=255)
                           tCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTAc                                                                                                          >  1:283313/1‑151 (MQ=255)
                            cAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACt                                                                                                         >  2:296364/1‑151 (MQ=255)
                             aaGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACtt                                                                                                        <  2:511654/151‑1 (MQ=255)
                              aGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACtt                                                                                                        >  2:619646/1‑150 (MQ=255)
                                     gATCCTCTGCCTGGTAGGGCCGCCGGGGTTAGGTAAAACCTCTCTTGGTCATTCCATGGCCAAAGCCACCGGGTCTAAATATGCCCGTCTGGCGCTGGGTGGCGTGCTTGATGAGGCGGAAATCCGTGGTCACCGCCGTACTTACATc                                                                                                   <  1:70207/148‑1 (MQ=255)
                                           ctGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTa                                                                                            <  2:199904/149‑1 (MQ=255)
                                                 ggTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAGGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCgg                                                                                      <  1:369907/149‑1 (MQ=255)
                                                          gccgGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATc                                                                           >  2:655946/1‑151 (MQ=255)
                                                            cgGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATcc                                                                          >  2:204777/1‑150 (MQ=255)
                                                             gggggTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATcc                                                                          >  1:466131/1‑149 (MQ=255)
                                                             gggggTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAg                                                                        <  2:258355/151‑1 (MQ=255)
                                                              ggggTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGa                                                                       <  2:249859/151‑1 (MQ=255)
                                                              ggggTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGACACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATcc                                                                          >  2:548769/1‑148 (MQ=255)
                                                                    ggtaAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGGGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATg                                                                  <  1:473695/150‑1 (MQ=255)
                                                                    ggtaAAACCTCCCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGGGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATgg                                                                 <  1:193746/151‑1 (MQ=255)
                                                                       aaaaCCTCTCTTGGTGAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGAAAACTGATCCAGAAAATGGCg                                                               >  2:560034/1‑150 (MQ=255)
                                                                             tctctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTgg                                                        >  2:654284/1‑151 (MQ=255)
                                                                               tctTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGc                                                      >  2:536130/1‑151 (MQ=255)
                                                                                  tGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGc                                                      <  2:168001/148‑1 (MQ=255)
                                                                                  tGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGGGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGt                                                    <  1:252840/150‑1 (MQ=255)
                                                                                   ggTCAGTCCATTGCCAAAGCCACCGGGCTTAAATATGTCGGTAGGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTgg                                                        <  1:583977/145‑1 (MQ=255)
                                                                                         tCCATTGCCAACGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGCGATGAAGCGGAAATCCGGGGTCACCCCCGCACTTACATCGGTTCTCTGCCGGGTAAAGTGATACAGAAAATGGCGAAAGGGGGCGTGAAAAAcc                                            >  1:537497/1‑151 (MQ=255)
                                                                                          ccATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAAccc                                           >  2:287551/1‑151 (MQ=255)
                                                                                          ccATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAAccc                                           >  1:216564/1‑151 (MQ=255)
                                                                                                ccAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGtt                                     <  2:469603/151‑1 (MQ=255)
                                                                                                    aGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTg                                 >  2:539087/1‑151 (MQ=255)
                                                                                                     gCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGc                                >  1:621312/1‑151 (MQ=255)
                                                                                                      cccccGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCt                               <  2:136495/148‑1 (MQ=255)
                                                                                                      ccACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCt                               <  2:61830/151‑1 (MQ=255)
                                                                                                      ccACCGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCt                               <  1:428816/151‑1 (MQ=255)
                                                                                                          cGGGCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGAt                           >  1:120803/1‑151 (MQ=255)
                                                                                                           gggCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATg                          <  2:483624/151‑1 (MQ=255)
                                                                                                             gCGTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATg                          >  2:358615/1‑149 (MQ=255)
                                                                                                               gTAAATATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATg                          >  2:664190/1‑147 (MQ=255)
                                                                                                                 aaaTATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCa                                                                         >  1:291414/1‑98 (MQ=255)
                                                                                                                 aaaTATGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCa                                                                         <  2:291418/98‑1 (MQ=255)
                                                                                                                  aaTATGTCCGTATGGCGCGGGGTGGCGTGCTTGATGAACCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGa                   <  1:456393/151‑1 (MQ=255)
                                                                                                                   atatGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGAc                  >  1:220411/1‑151 (MQ=255)
                                                                                                                     atGTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACaa                <  2:506095/151‑1 (MQ=255)
                                                                                                                       gTCCGTATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACaaa               <  2:50799/150‑1 (MQ=255)
                                                                                                                            tATGGCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtct         >  1:474748/1‑151 (MQ=255)
                                                                                                                               ggCGCTTGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtcttct      >  1:19631/1‑151 (MQ=255)
                                                                                                                               ggCGCTGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGtcttct      <  1:450892/151‑1 (MQ=255)
                                                                                                                                 cgatGGGTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTg     >  1:513561/4‑150 (MQ=255)
                                                                                                                                     gggTGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACa  <  2:442237/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                 
TTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTAAAACCTCTCTTGGTCAGTCCATTGCCAAAGCCACCGGGCGTAAATATGTCCGTATGGCGCTGGGCGGCGTGCGTGATGAAGCGGAAATCCGTGGTCACCGCCGTACTTACATCGGTTCTATGCCGGGTAAACTGATCCAGAAAATGGCGAAAGTGGGCGTGAAAAACCCGCTGTTCCTGCTCGATGAGATCGACAAAATGTCTTCTGACA  >  NC_000913/459900‑460181

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: