Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 403,332 A→G T18A (ACC→GCC)  psiF → PsiF family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913403,3320AG98.2% 192.7 / ‑4.3 54T18A (ACC→GCC) psiFPsiF family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (31/22);  total (32/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGC  >  NC_000913/403197‑403463
                                                                                                                                       |                                                                                                                                   
ggtTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCCCGGGCGCATACGCTCATGGTTAAAACATGAAGGGGGATGGTGCTATGAAAATAACATTCCTGGTTCCCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGt                                                                                                                                 <  2:355823/140‑1 (MQ=255)
     ttttGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGgcg                                                                                                                            <  2:976335/140‑1 (MQ=255)
           tGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCg                                                                                                                      >  2:761183/1‑140 (MQ=255)
               aGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCgagag                                                                                                                  <  2:4299/140‑1 (MQ=255)
                     cAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACttt                                                                                                            <  2:613303/140‑1 (MQ=255)
                            aaTGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCca                                                                                                     >  1:850883/1‑140 (MQ=255)
                                  aGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaa                                                                                               >  2:445860/1‑140 (MQ=255)
                                  aGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaa                                                                                               >  1:975331/1‑140 (MQ=255)
                                  aGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCacaacaa                                                                                               <  1:358771/140‑1 (MQ=255)
                                         gggCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTa                                                                                        <  1:1157994/140‑1 (MQ=255)
                                             gcATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAc                                                                                    <  2:850883/140‑1 (MQ=255)
                                              cATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcc                                                                                   >  1:765797/1‑140 (MQ=255)
                                                tACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                                 >  2:160911/1‑140 (MQ=255)
                                                         ggTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                                 >  2:308860/1‑131 (MQ=255)
                                                         ggTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGAcctc                                                                                 <  1:308860/131‑1 (MQ=255)
                                                            tAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATcagca                                                                     <  2:561254/140‑1 (MQ=255)
                                                             aaaaCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATcagcag                                                                    >  2:1431930/1‑140 (MQ=255)
                                                              aaaCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAgg                                                                   >  2:1007383/1‑140 (MQ=255)
                                                              aaaCATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAgg                                                                   >  1:653456/1‑140 (MQ=255)
                                                                   tGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACg                                                              >  1:267936/1‑140 (MQ=255)
                                                                   tGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACg                                                              >  2:786062/1‑140 (MQ=255)
                                                                    gAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACgg                                                             >  1:1322774/1‑140 (MQ=255)
                                                                      agagGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcg                                                           <  2:1322774/140‑1 (MQ=255)
                                                                         gggATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCCCCTGTAATCAGCAGGCGACGGCGCAg                                                        >  1:851807/1‑140 (MQ=255)
                                                                          ggATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAgg                                                       >  2:574138/1‑140 (MQ=255)
                                                                            aTGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCg                                                     <  1:754407/140‑1 (MQ=255)
                                                                            aTGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCg                                                     <  2:653456/140‑1 (MQ=255)
                                                                               gTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTg                                                  >  1:430483/1‑140 (MQ=255)
                                                                               gTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTg                                                  >  2:754916/1‑140 (MQ=255)
                                                                                  cTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGaaa                                               >  2:1402831/1‑140 (MQ=255)
                                                                                  cTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGaaa                                               >  1:1265283/1‑140 (MQ=255)
                                                                                    aTGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGTTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcg                                                           >  2:1065086/1‑126 (MQ=255)
                                                                                    aTGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGTTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGgcg                                                           <  1:1065086/126‑1 (MQ=255)
                                                                                    aTGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAgg                                             >  1:311431/1‑140 (MQ=255)
                                                                                    aTGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAgg                                             >  1:826275/1‑140 (MQ=255)
                                                                                     tGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGa                                                 <  2:733256/135‑1 (MQ=255)
                                                                                     tGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGa                                                 >  1:733256/1‑135 (MQ=255)
                                                                                     tGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAggg                                            >  1:739676/1‑140 (MQ=255)
                                                                                     tGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAggg                                            >  2:304519/1‑140 (MQ=255)
                                                                                      gAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAgggg                                           <  1:304519/140‑1 (MQ=255)
                                                                                         aaTAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATg                                        <  1:1402831/140‑1 (MQ=255)
                                                                                         aaTAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATg                                        <  2:172668/140‑1 (MQ=255)
                                                                                             aCATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCg                                    >  1:1254499/1‑140 (MQ=255)
                                                                                              cATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGt                                   <  2:311431/140‑1 (MQ=255)
                                                                                               aTTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTa                                  <  1:1007383/140‑1 (MQ=255)
                                                                                                 tACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAg                                <  2:975331/140‑1 (MQ=255)
                                                                                                    tGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGAcc                             <  1:445860/140‑1 (MQ=255)
                                                                                                             tGCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGt                    >  2:311359/1‑140 (MQ=255)
                                                                                                              gCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTg                   >  1:1456582/1‑140 (MQ=255)
                                                                                                              gCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTg                   >  1:1266610/1‑140 (MQ=255)
                                                                                                              gCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTg                   >  2:1306065/1‑140 (MQ=255)
                                                                                                              gCTTTTCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTg                   >  1:1137850/1‑140 (MQ=255)
                                                                                                                  ttCGGTCTGGTTTTTTTAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTg               >  1:1449158/1‑140 (MQ=255)
                                                                                                                               ttttAACCGCCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGc  <  2:765797/140‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
GGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAACTTTAACCCCACAACAACAGCGTATGACCTCCTGTAATCAGCAGGCGACGGCGCAGGCGTTGAAAGGGGATGCTCGTAAGACCTACATGAGTGATTGCCTGAAGAACAGC  >  NC_000913/403197‑403463

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