Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,383,631 T→C 7.5% D467G (GAC→GGC)  proS ← proline‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,383,6310TC 7.5% 150.0 / 2.6 53D467G (GAC→GGC) proSproline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (24/25);  new base C (2/2);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCT  >  NC_002947/1383492‑1383773
                                                                                                                                           |                                                                                                                                              
ggTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGACGCCCTGGCTCCGTTCCAga                                                                                                                                       <  8:94511/149‑1 (MQ=255)
 gTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAg                                                                                                                                        <  3:254073/147‑1 (MQ=255)
         gTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGCGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGccc                                                                                                                               <  5:11540/148‑1 (MQ=255)
                 cgtcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTGGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCa                                                                                                                      <  4:205297/149‑1 (MQ=255)
                   tcCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGACGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATa                                                                                                                    <  3:201535/149‑1 (MQ=255)
                    ccAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCa                                                                                                                      >  2:357818/1‑146 (MQ=255)
                       gcagcaCTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGctc                                                                                                                >  7:188416/1‑149 (MQ=255)
                             cTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTc                                                                                                           >  2:226424/1‑148 (MQ=255)
                              ttCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCAt                                                                                                                     <  5:99728/137‑1 (MQ=255)
                              ttCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCAt                                                                                                                     >  6:99728/1‑137 (MQ=255)
                              ttCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGa                                                                                                         >  6:255326/1‑149 (MQ=255)
                                cATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGa                                                                                                         <  8:209076/147‑1 (MQ=255)
                                 aTAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                         >  3:159207/1‑114 (MQ=255)
                                 aTAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCa                                                                                                                                         <  4:159207/114‑1 (MQ=255)
                                     ccGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATgg                                                                                                      >  6:380963/1‑145 (MQ=255)
                                      cggcggCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCt                                                                                                 <  7:310234/149‑1 (MQ=255)
                                       ggcggcAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTg                                                                                                >  3:123586/1‑149 (MQ=255)
                                             agtcagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACca                                                                                          >  8:95448/1‑149 (MQ=255)
                                                cagtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACcaca                                                                                        <  1:357818/148‑1 (MQ=255)
                                                 agtTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTAGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACg                                                                                      >  4:286951/1‑149 (MQ=255)
                                                   ttCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACggg                                                                                    <  4:123586/149‑1 (MQ=255)
                                                      ggCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                                    <  3:68529/130‑1 (MQ=255)
                                                      ggCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                                    >  4:68529/1‑130 (MQ=255)
                                                      ggCAGACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGac                                                                                  <  3:286951/148‑1 (MQ=255)
                                                        cATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACg                                                                               >  6:172702/1‑149 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                 >  2:54733/1‑114 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                                 <  1:54733/114‑1 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                              >  1:249271/1‑149 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                              >  2:16262/1‑149 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                              <  5:172702/149‑1 (MQ=255)
                                                                ttGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCGGCTGCAACCACACGGGACACGCCAATGcc                                                                       >  8:319854/1‑149 (MQ=255)
                                                                   tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                    <  5:255326/149‑1 (MQ=255)
                                                                   tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                    <  1:373023/149‑1 (MQ=255)
                                                                   tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                                    <  7:95448/149‑1 (MQ=255)
                                                                                    cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                        <  3:109764/144‑1 (MQ=255)
                                                                                    cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                        >  4:109764/1‑144 (MQ=255)
                                                                                    cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTa                                                   <  1:256694/149‑1 (MQ=255)
                                                                                           gTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                        <  8:151655/137‑1 (MQ=255)
                                                                                           gTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                        >  7:151655/1‑137 (MQ=255)
                                                                                                gTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTg                                        <  6:168683/148‑1 (MQ=255)
                                                                                                   gcgcAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGt                                    >  7:163840/1‑149 (MQ=255)
                                                                                                      cAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGTTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTtct                                 >  1:186136/1‑149 (MQ=255)
                                                                                                          gggACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGTTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGcc                             <  2:186136/149‑1 (MQ=255)
                                                                                                               cAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCa                        >  4:34430/1‑149 (MQ=255)
                                                                                                                   gCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTg                    <  8:163840/149‑1 (MQ=255)
                                                                                                                    cAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTgg                   <  5:380963/149‑1 (MQ=255)
                                                                                                                       tCTGGAACGACGCCCTGGCTCCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCACTGCCGTAGCAACCCAGCGACAGTCCGACCGGCTTGCCGGTCTCGCCCAGCACCTGGccc                >  4:268970/1‑147 (MQ=255)
                                                                                                                       tCTGGAACGACGCCCTGGCTCCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTtct                                 <  7:89805/132‑1 (MQ=255)
                                                                                                                       tCTGGAACGACGCCCTGGCTCCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTtct                                 >  8:89805/1‑132 (MQ=255)
                                                                                                                                gAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACCGGACCCGCCAATGCCGTAGCAACCCATAGACAGTACGACCGGCTTGCCGTTCTCGCCACGCACCTGTCACTTCAGAGcc       >  6:319895/1‑149 (MQ=255)
                                                                                                                                  gCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTc     >  7:334407/1‑149 (MQ=255)
                                                                                                                                     aGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCt  >  3:165958/1‑149 (MQ=255)
                                                                                                                                     aGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCt  <  1:16262/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                              
GGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCCTCGCT  >  NC_002947/1383492‑1383773

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