Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,383,621 A→C 8.0% A470A (GCT→GCG)  proS ← proline‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,383,6210AC 8.0% 154.0 / 4.0 50A470A (GCT→GCG) proSproline‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (22/24);  new base C (2/2);  total (24/26)
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.86e-01

GGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCC  >  NC_002947/1383492‑1383768
                                                                                                                                 |                                                                                                                                                   
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                                                                                               >  4:68529/1‑130 (MQ=255)
                                                      ggCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCa                                                                                               <  3:68529/130‑1 (MQ=255)
                                                      ggCAGACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGac                                                                             <  3:286951/148‑1 (MQ=255)
                                                        cATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACg                                                                          >  6:172702/1‑149 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                            <  1:54733/114‑1 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGct                                                                                                            >  2:54733/1‑114 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                         >  1:249271/1‑149 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                         <  5:172702/149‑1 (MQ=255)
                                                         aTACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGc                                                                         >  2:16262/1‑149 (MQ=255)
                                                                ttGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCGGCTGCAACCACACGGGACACGCCAATGcc                                                                  >  8:319854/1‑149 (MQ=255)
                                                                   tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                               <  1:373023/149‑1 (MQ=255)
                                                                   tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                               <  5:255326/149‑1 (MQ=255)
                                                                   tCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTa                                                               <  7:95448/149‑1 (MQ=255)
                                                                                    cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                   >  4:109764/1‑144 (MQ=255)
                                                                                    cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                   <  3:109764/144‑1 (MQ=255)
                                                                                    cACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTa                                              <  1:256694/149‑1 (MQ=255)
                                                                                           gTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                   >  7:151655/1‑137 (MQ=255)
                                                                                           gTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGa                                                   <  8:151655/137‑1 (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)
                                                                                                                                 |                                                                                                                                                   
GGTTTTCTTGTCGCGGTCGTCCAGCAGCACTTCATAGCCGGCGGCAGTCAGTTCGGCATACAGCTTGTCGGTTGCTTCGCGGACCACGTCGGTTTCGTAGCGCAGTGGGACCAGGGCAATCTGGAACGGAGCCAGGGCGTCGTTCCAGATGATGCCCTTGTCGTCATAGCTCTGCTCGATGGCAGCTGCAACCACACGGGACACGCCAATGCCGTAGCAACCCATCGACAGTACGACCGGCTTGCCGTTCTCGCCCAGCACCTGGCACTTCAGCGCC  >  NC_002947/1383492‑1383768

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