Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,303,969 G→A 8.8% R202H (CGC→CAC) ‡ cysS → cysteine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,303,9690GA 8.8% 158.8 / 2.7 56R202H (CGC→CAC) ‡cysScysteine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base G (23/28);  new base A (3/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGC  >  NC_002947/3303827‑3304116
                                                                                                                                              |                                                                                                                                                   
cTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCggg                                                                                                                                               <  1:124056/149‑1 (MQ=255)
   tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGc                                                                                                                                              <  2:87177/147‑1 (MQ=255)
   tCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGcc                                                                                                                                             <  6:227716/148‑2 (MQ=255)
        cAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGcaca                                                                                                                                       >  4:5587/1‑149 (MQ=255)
             ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAg                                                                                                                                  <  7:40255/149‑1 (MQ=255)
             ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAg                                                                                                                                  <  7:100503/149‑1 (MQ=255)
             ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGc                                                                                                                                              <  5:117997/137‑1 (MQ=255)
             ggATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGCGACCCGGGc                                                                                                                                              >  6:117997/1‑137 (MQ=255)
              gATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGt                                                                                                                                 <  3:69328/149‑1 (MQ=255)
               aTCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGAC‑TAGCCAAGGATGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGTGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCTCGGGCTGGCACATCGAGTGc                                                                                                                               <  2:251375/149‑1 (MQ=255)
                      aCCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTc                                                                                                                         <  3:106645/149‑1 (MQ=255)
                       ccTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTc                                                                                                                         <  5:143394/148‑1 (MQ=255)
                          gCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                                                                                                     >  3:292478/1‑149 (MQ=255)
                             tATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTg                                                                                                                                            >  4:10121/1‑123 (MQ=255)
                             tATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTg                                                                                                                                            <  3:10121/123‑1 (MQ=255)
                              aTCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                                                                                                   >  5:61467/1‑131 (MQ=255)
                              aTCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGa                                                                                                                                   <  6:61467/131‑1 (MQ=255)
                               tCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAcc                                                                                                                >  6:44843/1‑149 (MQ=255)
                               tCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGGCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCAc                                                                                                                 <  4:9294/148‑1 (MQ=255)
                                  gTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgc                                                                                                             <  4:292478/149‑1 (MQ=255)
                                    gCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACctgctg                                                                                                           >  3:72721/1‑149 (MQ=255)
                                         cATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTgg                                                                                                      <  3:285293/149‑1 (MQ=255)
                                                gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                                                                                                     <  3:37686/127‑1 (MQ=255)
                                                gTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGt                                                                                                                     >  4:37686/1‑127 (MQ=255)
                                                     cgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCAt                                                                                                                       >  3:150417/1‑120 (MQ=255)
                                                     cgaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCAt                                                                                                                       <  4:150417/120‑1 (MQ=255)
                                                      gaAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGa                                                                                         >  5:104405/1‑149 (MQ=255)
                                                                    cccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggcg                                                                           >  8:61428/1‑149 (MQ=255)
                                                                    cccGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggcg                                                                           >  3:220542/1‑149 (MQ=255)
                                                                       gCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGGGCTCGGTCATGGCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCa                                                                        <  3:308258/149‑1 (MQ=255)
                                                                             cTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGAcc                                                                  >  6:14358/1‑149 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgaga                                                                                                <  8:308720/114‑1 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCgaga                                                                                                >  7:308720/1‑114 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCga                                                                                                  <  4:22531/112‑1 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCga                                                                                                  >  3:22531/1‑112 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAg                                                             >  7:210756/1‑149 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCACGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggc                                                                            >  5:240935/1‑134 (MQ=255)
                                                                                  tCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCACGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTggc                                                                            <  6:240935/134‑1 (MQ=255)
                                                                                            gggCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGAcc                                                                  >  1:242134/1‑134 (MQ=255)
                                                                                            gggCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGAcc                                                                  <  2:242134/134‑1 (MQ=255)
                                                                                                  cAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCacgag                                              <  6:104405/148‑1 (MQ=255)
                                                                                                               ccAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCa                                <  5:56244/149‑1 (MQ=255)
                                                                                                                     tGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGc                            >  6:169554/1‑147 (MQ=255)
                                                                                                                         aaTCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAggcg                      <  4:220542/149‑1 (MQ=255)
                                                                                                                          aTCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTgg                                                               <  2:313729/107‑1 (MQ=255)
                                                                                                                          aTCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTgg                                                               >  1:313729/1‑107 (MQ=255)
                                                                                                                            cGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAggcggc                    >  6:78167/1‑148 (MQ=255)
                                                                                                                              ccATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGTTTGCCCAGAGCGAggcg                      <  5:65493/144‑1 (MQ=255)
                                                                                                                              ccATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGTTTGCCCAGAGCGAggcg                      >  6:65493/1‑144 (MQ=255)
                                                                                                                                tgccagCCCGGGCGACCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCGTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGCTGGCCCAGAGCGAGGCGGCCCCCg               >  6:200191/6‑149 (MQ=255)
                                                                                                                                aTGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTg                                                                >  7:15170/1‑100 (MQ=255)
                                                                                                                                aTGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTg                                                                <  8:15170/100‑1 (MQ=255)
                                                                                                                                   cagCCCGGGCGACCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGc             >  5:264721/3‑148 (MQ=255)
                                                                                                                                    ggCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCaa           <  5:44843/149‑1 (MQ=255)
                                                                                                                                      cccGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAg          >  7:153036/1‑148 (MQ=255)
                                                                                                                                       ccGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCa        <  5:14358/149‑1 (MQ=255)
                                                                                                                                              gCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGc  >  7:196281/1‑148 (MQ=255)
                                                                                                                                              gCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGc  >  3:13929/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CTGTCGCGCAAGAGGATCGAAGACCTGCGTATCGGTGCACGCATCGAGGTCGACGAAGCCAAGCAGGACCCGCTCGACTTCGTCCTGTGGAAGGGCGTCAAGCCGGGCGAGCCAAGCTGGGAATCGCCATGGGGCCCGGGTCGCCCGGGCTGGCACATCGAGTGCTCGGTCATGTCCACCTGCTGCCTGGGCGAGAGCTTCGACATTCACGGTGGCGGCAGCGACCTGGAGTTCCCGCACCACGAGAACGAGATTGCCCAGAGCGAGGCGGCCACCGGCAAGCAGTACGC  >  NC_002947/3303827‑3304116

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