Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,750,486 A→T 18.8% L378H (CTT→CAT) ‡ thiI ← tRNA sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,750,4860AT18.8% 110.3 / 21.3 48L378H (CTT→CAT) ‡thiItRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base T (4/5);  total (26/22)
Fisher's exact test for biased strand distribution p-value = 7.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

AGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCGGT  >  NC_002947/5750339‑5750628
                                                                                                                                                   |                                                                                                                                              
aGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                                              <  1:1134771/150‑1 (MQ=255)
    ccAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgc                                                                                                                                          >  1:183091/1‑150 (MQ=255)
     cAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAgcgcg                                                                                                                                         >  1:154008/1‑150 (MQ=255)
       gagGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgcc                                                                                                                                       <  1:300079/150‑2 (MQ=255)
        agGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgcca                                                                                                                                      <  2:1066683/150‑3 (MQ=255)
                          aTCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                                                    <  1:703227/150‑1 (MQ=255)
                                      aaTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                                                        <  1:361710/150‑1 (MQ=255)
                                      aaTGTCGATGACGATCTGGCATGGCAGAGCCTCGGCTACTTCCTCGTTTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttcttc                                                                                                        <  2:63044/150‑1 (MQ=255)
                                        tGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGt                                                                                                      <  1:1143011/150‑1 (MQ=255)
                                            gATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                                                     >  1:1116348/1‑112 (MQ=255)
                                            gATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                                                     <  2:1116348/115‑4 (MQ=255)
                                            gATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTc                                                                                                  >  1:301663/1‑150 (MQ=255)
                                               gACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGcc                                                                                                                            >  1:732537/1‑121 (MQ=255)
                                               gACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGcc                                                                                                                            <  2:732537/121‑1 (MQ=255)
                                                  gATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCtttt                                                                                                            <  2:454215/134‑1 (MQ=255)
                                                  gATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCtttt                                                                                                            >  1:454215/1‑134 (MQ=255)
                                                    tCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCtt                                                                                                              >  1:340682/1‑130 (MQ=255)
                                                    tCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCtt                                                                                                              <  2:340682/130‑1 (MQ=255)
                                                      tggcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGa                                                                                        >  2:298759/1‑150 (MQ=255)
                                                       ggcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGAc                                                                                       <  2:921159/150‑1 (MQ=255)
                                                        gcctggcAGCGCCTCGGCCGCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGTAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACg                                                                                      >  2:366179/1‑150 (MQ=255)
                                                        gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                                              >  2:50941/1‑94 (MQ=255)
                                                        gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                                              <  1:50941/94‑1 (MQ=255)
                                                        gcctggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGAcc                                                                                      >  2:789526/1‑149 (MQ=255)
                                                             gcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCtt                                                                                   >  1:1112165/1‑148 (MQ=255)
                                                             gcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCtt                                                                                   <  2:1112165/148‑1 (MQ=255)
                                                                   ccTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATggg                                                                           >  1:1185284/1‑150 (MQ=255)
                                                                               ccTCGAGTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTa                                                               <  2:557089/150‑1 (MQ=255)
                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                                                     >  1:1178002/1‑67 (MQ=38)
                                                                                         gATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                                                     <  2:1178002/70‑4 (MQ=38)
                                                                                              ccTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAcg                                                    >  2:741045/1‑146 (MQ=255)
                                                                                              ccTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCAcg                                                    <  1:741045/146‑1 (MQ=255)
                                                                                               cTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                                                    >  2:842741/1‑81 (MQ=255)
                                                                                               cTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGaa                                                                                                                    <  1:842741/81‑1 (MQ=255)
                                                                                                    ccAGCTCATCGATCTCATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTAtt                                          >  2:943110/1‑150 (MQ=255)
                                                                                                    ccAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTAtt                                          >  2:1105123/1‑150 (MQ=255)
                                                                                                     cAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTc                                         >  1:732412/1‑150 (MQ=255)
                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt                                   <  2:301663/150‑1 (MQ=255)
                                                                                                           aTCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAt                                   >  2:983881/1‑150 (MQ=255)
                                                                                                            tCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCAtg                                  >  1:522030/1‑150 (MQ=255)
                                                                                                                 cacaTGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCGTGGGGGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCtt                             >  1:378992/1‑150 (MQ=255)
                                                                                                                       aTCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCaaa                       <  2:747562/150‑1 (MQ=255)
                                                                                                                         cGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGt                     <  2:1185284/150‑1 (MQ=255)
                                                                                                                          gATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTc                    >  1:710989/1‑150 (MQ=255)
                                                                                                                                      tggagcgcgcTCTTGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGCGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGAt        >  1:511735/5‑150 (MQ=255)
                                                                                                                                          gcgcgcTCTTGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                                                  <  1:1153050/72‑1 (MQ=255)
                                                                                                                                          gcgcgcTCTTGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                                                  >  2:1153050/1‑72 (MQ=255)
                                                                                                                                          gcgcgcTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                                                    >  2:571619/1‑86 (MQ=255)
                                                                                                                                          gcgcgcTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCgg                                                                    <  1:571619/86‑1 (MQ=255)
                                                                                                                                            gcgcTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTCCCGAATTTACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCGGt  <  2:1012015/150‑1 (MQ=255)
                                                                                                                                            gcgcTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAATTCCGCGGTACCGATTTCggg  >  1:608890/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AGCACCAGAGGTTCGTCTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCGGT  >  NC_002947/5750339‑5750628

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