Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,836,736 A→T 7.1% T62S (ACG→TCG)  ttcA → tRNA 2‑thiocytidine(32) synthetase TtcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,836,7360AT 7.1% 160.1 / 2.5 56T62S (ACG→TCG) ttcAtRNA 2‑thiocytidine(32) synthetase TtcA
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base T (2/2);  total (27/29)
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.55e-01

CTGGCGAAGCCATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTAAAGGAGCTGGTACCCGAAGG  >  NC_002947/1836608‑1836884
                                                                                                                                |                                                                                                                                                    
cTGGCGAAGCCATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGt                                                                                                                                   >  8:162383/1‑148 (MQ=255)
   gCGAAGCCATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAAc                                                                                                                               <  1:126250/149‑1 (MQ=255)
        gcaaTCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGa                                                                                                                          <  8:261116/146‑1 (MQ=255)
        gCCATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATCATCGACGTTCAGTTGAACCTGCAGAAGGTGGCACCGATAACGTTCTAGATCGTTGCGGTGAACATGGa                                                                                                                          >  7:141064/1‑149 (MQ=255)
         ccATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGAc                                                                                                                         <  2:90699/149‑1 (MQ=255)
         ccATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGAc                                                                                                                         >  8:278871/1‑149 (MQ=255)
                ctaCTTCATCATGATCGAGGATGGCGACAAGGTCATGGTCTGTCTGTCCGGCGGCAAGGTCAGCTACACCATGCTCGACGTTGTGTTTCACCTGCAGAAGGTGGCACCGATCACGTTCGTGATCGTTGCGGTGAACATGGACGAGAAAc                                                                                                                  <  4:316336/147‑1 (MQ=255)
                     acaacaTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCg                                                                                                             <  6:266721/149‑1 (MQ=255)
                       aacaTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACAt                                                                                                                             >  8:291673/1‑131 (MQ=255)
                       aacaTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTCAACAt                                                                                                                             <  7:291673/131‑1 (MQ=255)
                        acaTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCgg                                                                                                            >  5:322408/1‑147 (MQ=255)
                        acaTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCgg                                                                                                            <  6:322408/147‑1 (MQ=255)
                         caTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCTCGAACGTGATCGGTg                                                                                                                                       >  7:71302/1‑119 (MQ=255)
                         caTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCTCGAACGTGATCGGTg                                                                                                                                       <  8:71302/119‑1 (MQ=255)
                          aTGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCtt                                                                                                        <  7:260604/149‑1 (MQ=255)
                            gATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGTACCAGAAACAGCCGGGCTTcc                                                                                                      >  3:121299/1‑149 (MQ=255)
                               cgcggATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCCCCTGCAGAAGGTGGCACCGCTCCCGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCt                                                                                                    >  5:27662/4‑148 (MQ=255)
                                gttgATGGCGACCAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCCTGCTCGACGTTCTGTTGCAACTGCAGAAGGGGGCACCGATCACGTTCGAGATCGGGGCGGTGAACAGTGACCCGAAACAGCCGGGCGTCCCTGa                                                                                                  <  6:27662/146‑1 (MQ=255)
                                gAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCt                                                                                                    >  6:84424/1‑147 (MQ=255)
                                             aGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCg                                                                                     <  6:177579/149‑1 (MQ=255)
                                              ggTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCgg                                                                                    <  7:162383/149‑1 (MQ=255)
                                               gTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAg                                                                                                                      >  1:271517/1‑114 (MQ=255)
                                               gTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAg                                                                                                                      <  2:271517/114‑1 (MQ=255)
                                                tCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATTCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGc                                                                                   >  6:207042/1‑148 (MQ=255)
                                                 cATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTg                                                                                        >  1:125584/1‑142 (MQ=255)
                                                       cTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTc                                                                            >  5:118812/1‑148 (MQ=255)
                                                           cTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCg                                                                                     >  1:30118/1‑135 (MQ=255)
                                                           cTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCg                                                                                     <  2:30118/135‑1 (MQ=255)
                                                                       aaGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTggggg                                                               >  2:135255/1‑145 (MQ=255)
                                                                       aaGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTggggg                                                               <  1:135255/145‑1 (MQ=255)
                                                                             aGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACca                                                     >  1:17913/1‑149 (MQ=255)
                                                                                  cacCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACAtc                                                 <  8:273798/148‑1 (MQ=255)
                                                                                         cTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTTGGGGTCGAGTACCACATCGTCGAGaa                                         <  5:207042/149‑1 (MQ=255)
                                                                                         cTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCgag                                           >  1:306713/1‑147 (MQ=255)
                                                                                         cTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGaa                                         <  6:118812/149‑1 (MQ=255)
                                                                                          tCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACAtcgt                                               <  1:172416/142‑1 (MQ=255)
                                                                                          tCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACAtcgt                                               >  2:172416/1‑142 (MQ=255)
                                                                                          tCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAg                                        >  4:289484/1‑149 (MQ=255)
                                                                                           cGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAg                                        >  6:21074/1‑148 (MQ=255)
                                                                                               gTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATAc                                   <  5:84424/149‑1 (MQ=255)
                                                                                                ttCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATAcc                                  <  2:17913/149‑1 (MQ=255)
                                                                                                  cTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTa                                >  3:50138/1‑149 (MQ=255)
                                                                                                     ttGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACt                              <  5:21074/148‑1 (MQ=255)
                                                                                                        cACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTcc                                                                                                      <  3:254468/73‑1 (MQ=255)
                                                                                                        cACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTcc                                                                                                      >  4:254468/1‑73 (MQ=255)
                                                                                                         aCCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGa                                                           <  3:239656/115‑1 (MQ=255)
                                                                                                         aCCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGa                                                           >  4:239656/1‑115 (MQ=255)
                                                                                                              cAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTaaa                    <  2:306713/149‑1 (MQ=255)
                                                                                                              aaGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGCTGAACATGGACCAGAAACAGACGGGCTTCCCTGACCACGTACTGCATGCCTACCTCAAGGAGCTGGGCTTCGAGTACCACATCGTCGAGAACCATACCTACTCGGAGGTaa                     >  4:21070/2‑148 (MQ=255)
                                                                                                                      cgCAACGATCTCGAACGTGATCGGTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGaa                                         <  2:120947/119‑1 (MQ=255)
                                                                                                                      cgCAACGATCTCGAACGTGATCGGTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGaa                                         >  1:120947/2‑120 (MQ=255)
                                                                                                                       gCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGt                                                         <  1:270492/103‑1 (MQ=255)
                                                                                                                       gCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGt                                                         >  2:270492/1‑103 (MQ=255)
                                                                                                                       gCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTAAAGGAGCTgg            >  3:283456/1‑148 (MQ=255)
                                                                                                                         aCCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTAAAGGAGCTGGTAc         <  3:289484/149‑1 (MQ=255)
                                                                                                                          ccGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTAAAGGAGCTGGTAcc        <  3:21070/149‑1 (MQ=255)
                                                                                                                                aCGGTCGAGATCGTTGCGGTGAACATGCACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCTGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTAAAGGAGCTGGTACCCGAAgg  <  8:141064/149‑1 (MQ=255)
                                                                                                                                |                                                                                                                                                    
CTGGCGAAGCCATCACCGACTACAACATGATCGAGGATGGCGACAAGGTCATGGTCTGCCTGTCCGGCGGCAAGGACAGCTACACCATGCTCGACGTTCTGTTGCACCTGCAGAAGGTGGCACCGATCACGTTCGAGATCGTTGCGGTGAACATGGACCAGAAACAGCCGGGCTTCCCTGAGCACGTACTGCCGGCCTACCTCAAGGAGCTGGGGGTCGAGTACCACATCGTCGAGAAGGATACCTACTCGGTGGTAAAGGAGCTGGTACCCGAAGG  >  NC_002947/1836608‑1836884

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