Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 459,751 G→T 12.7% A288A (GCG→GCT)  lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,7510GT12.7% 161.8 / 11.4 55A288A (GCG→GCT) lonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (11/37);  new base T (3/4);  total (14/41)
Fisher's exact test for biased strand distribution p-value = 3.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAA  >  NC_000913/459612‑459879
                                                                                                                                           |                                                                                                                                
gTGAAATGGACGACGCGCCGGACGAAAACGACGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCg                                                                                                                                  <  2:980451/140‑1 (MQ=255)
    aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCTACCg                                                                                                                              >  2:1006084/1‑140 (MQ=255)
    aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCTACCg                                                                                                                              >  2:1006044/1‑140 (MQ=255)
    aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCTACCg                                                                                                                              >  2:1005925/1‑140 (MQ=255)
    aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCg                                                                                                                              >  1:21523/1‑140 (MQ=255)
    aaTGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCg                                                                                                                              >  1:21567/1‑140 (MQ=255)
         acgacgCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTg                                                                                                                         <  2:612746/140‑1 (MQ=255)
         acgacgCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTg                                                                                                                         >  1:947248/1‑140 (MQ=255)
           gacgCGCCGGACGAAAACGACGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCg                                                                                                                       <  2:862866/140‑1 (MQ=255)
                          aaCGAACCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAACCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTg                                                                                                        <  2:1194045/140‑1 (MQ=255)
                           acgacgCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg                                                                                                       <  1:877415/140‑1 (MQ=255)
                           aCGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTgg                                                                                                       <  1:877219/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCTACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  1:1006044/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCTACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  1:1005925/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCTGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCTACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  1:1006084/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCACAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGAAGTCGGCAGAAGCGACCGTAGTGCGAGGTTATATCGACTGGa                                                                                                      >  1:851350/1‑140 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  2:503904/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  2:457765/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  2:457137/140‑1 (MQ=255)
                            cGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGa                                                                                                      <  2:504375/140‑1 (MQ=255)
                               aGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATgg                                                                                                   >  1:210263/1‑140 (MQ=255)
                                           cAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGt                                                                                       >  1:984660/1‑140 (MQ=255)
                                           cAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGt                                                                                       >  1:147615/1‑140 (MQ=255)
                                           cAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGt                                                                                       >  1:1228802/1‑140 (MQ=255)
                                                      cggcgAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGt                                                                            >  2:1204476/1‑140 (MQ=255)
                                                             gATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGATGAAAATTATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAgg                                                                     <  2:1130054/140‑1 (MQ=255)
                                                              aTGCCGAACGAGGCAAAAGTGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAATTGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:815541/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:643253/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:947248/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:102592/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:456867/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:451037/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:433958/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:433667/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:530590/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:813778/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:270910/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:201978/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:186816/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:1177663/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:688044/140‑1 (MQ=255)
                                                              aTGCCGAAAGAGGCAAAAGAGAAACCGGAAGCAGAGTTGCAGAAGCTGAAAACGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTTGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGt                                                                    <  2:885532/140‑1 (MQ=255)
                                                                          gCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCt                                                        >  2:1092183/1‑140 (MQ=255)
                                                                                    aaGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGgcg                                              <  2:296614/140‑1 (MQ=255)
                                                                                          aaGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGaa                                        <  1:29697/140‑1 (MQ=255)
                                                                                               tagTTCCAGAAGCTCAAAATGATGTCTCCGATTTCGGCAGAAGCGACCGTAGTTCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCt                                   <  2:395236/139‑1 (MQ=255)
                                                                                                      aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATAcc                            <  2:299829/140‑1 (MQ=255)
                                                                                                      aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATAcc                            <  1:1053728/140‑1 (MQ=255)
                                                                                                      aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATAcc                            <  2:1148514/140‑1 (MQ=255)
                                                                                                      aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATAcc                            <  2:1148420/140‑1 (MQ=255)
                                                                                                      aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATAcc                            <  1:422375/140‑1 (MQ=255)
                                                                                                             gAAAATGATGTCTCCGATGTCGGCAGAAGCTACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAATGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCAtt                     <  1:338282/140‑1 (MQ=255)
                                                                                                             gAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATCACGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGTCGCAGGAAATCCTTGATACCGACCAtt                     <  2:95829/140‑1 (MQ=255)
                                                                                                             gAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCAtt                     <  2:1166173/140‑1 (MQ=255)
                                                                                                                        ctcCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAg          >  1:579788/1‑140 (MQ=255)
                                                                                                                                tCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGaa  <  1:329937/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                
GTGAAATGGACGACGCGCCGGACGAAAACGAAGCCCTGAAGCGCAAAATCGACGCGGCGAAGATGCCGAAAGAGGCAAAAGAGAAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAA  >  NC_000913/459612‑459879

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑