Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,750,486 A→T 12.5% 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,4860AT12.5% 106.9 / 8.1 40L378H (CTT→CAT) ‡thiItRNA sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (19/16);  new base T (3/2);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.47e-01

CTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCGGT  >  NC_002947/5750355‑5750628
                                                                                                                                   |                                                                                                                                              
cTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACgg                                                                                                                              >  1:762940/1‑150 (MQ=255)
     tGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATg                                                                                                                         >  2:695734/1‑150 (MQ=255)
      gggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGt                                                                                                                        >  1:95248/1‑150 (MQ=255)
      gggCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGt                                                                                                                        <  2:782964/150‑1 (MQ=255)
              gggTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGc                                                                                                                <  1:59723/150‑1 (MQ=255)
                gTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCtt                                                                                                              <  2:762940/150‑1 (MQ=255)
                gTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTAATTGAGCACATGATTCTTCCAATTTAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCtt                                                                                                              <  1:696048/150‑1 (MQ=255)
                   gCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTttc                                                                                                           >  2:247560/1‑150 (MQ=255)
                      aaTGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCAAGGCCATGTCGAACTGCTTttcttc                                                                                                        <  1:422293/150‑1 (MQ=255)
                        tGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGt                                                                                                      <  1:152807/150‑1 (MQ=255)
                             aTGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATTTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTcc                                                                                                 >  2:752110/1‑150 (MQ=255)
                                   aTCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCgg                                                                                           >  2:664133/1‑150 (MQ=255)
                                    tCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGt                                                                                          >  2:724460/1‑150 (MQ=255)
                                            ggcAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTg                                                                                  <  1:615793/150‑1 (MQ=255)
                                                         ccACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCAGGGGTGGTCg                                                                     >  2:73150/1‑150 (MQ=255)
                                                           aCTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCggg                                                                   >  2:718702/1‑150 (MQ=255)
                                                               ccTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                                              <  1:356985/71‑1 (MQ=255)
                                                               ccTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAaga                                                                                                                                              >  2:356985/1‑71 (MQ=255)
                                                                 tCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTAcc                                                             >  1:272856/1‑150 (MQ=255)
                                                                   gATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCga                                                           >  2:138405/1‑150 (MQ=255)
                                                                    aTTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCAAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGc                                                                                             <  2:128961/115‑1 (MQ=255)
                                                                    aTTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCAAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGc                                                                                             >  1:128961/1‑115 (MQ=255)
                                                                     tttCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCgaga                                                         >  2:383249/1‑150 (MQ=255)
                                                                               cTTGCCCAGCTCATTGATCACATAATCGATAGAAATGAACTTAGCGCGCTCAAGAGCGCGCTACAGCACGGCCATGTCGAACTGCTTTTCGTCGTGCTCCATGCGCTGACGCTTGGCATTGGTGTTCGGGTTTACCGAGATCGCGCCGCa                                               >  2:141634/1‑150 (MQ=255)
                                                                                ttGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                                                     <  1:656934/63‑4 (MQ=37)
                                                                                ttGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCTTGAgcgcgccaa                                                                                                                                     >  2:656934/1‑60 (MQ=37)
                                                                                ttGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAg                                              >  1:19612/1‑150 (MQ=255)
                                                                                     cAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTc                                         <  2:95248/150‑1 (MQ=255)
                                                                                      aGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTcc                                        <  1:247560/150‑1 (MQ=255)
                                                                                                       aTCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTg                                                                                         >  2:193714/1‑84 (MQ=255)
                                                                                                       aTCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTg                                                                                         <  1:193714/84‑1 (MQ=255)
                                                                                                         cGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCTTGCTCCATGCGGt                                                                                          <  2:33184/81‑1 (MQ=255)
                                                                                                         cGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGt                                                                                          >  1:33184/1‑81 (MQ=255)
                                                                                                                aaTGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGt              <  1:724460/150‑1 (MQ=255)
                                                                                                                      cTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGAt        >  1:584168/1‑150 (MQ=255)
                                                                                                                        tGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTAcccctat      <  1:752110/150‑6 (MQ=255)
                                                                                                                        gagcgcgcTCTTGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCcgcg                >  2:477174/3‑140 (MQ=255)
                                                                                                                        gagcgcgcTCTTGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCcgcg                <  1:477174/138‑1 (MQ=255)
                                                                                                                          gcgcgcTCTTGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCg    >  2:775247/1‑150 (MQ=255)
                                                                                                                            gcgcTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCGGt  <  1:718702/150‑1 (MQ=255)
                                                                                                                                tCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCgg               <  2:412210/133‑1 (MQ=255)
                                                                                                                                tCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCgg               >  1:412210/1‑133 (MQ=255)
                                                                                                                                   |                                                                                                                                              
CTTCCTGGGCATCGGGGTGGCGAATGTCGATGACGATCTGGCCTGGCAGCGCCTCGGCCACTTCCTCGATTTCGATGTCCTTGCCCAGCTCATCGATCACATGATCGATGGAAATGAACTTGGCGCGCTCAAGAGCGCGCTCCAGCACGGCCATGTCGAACTGCTTTTCTTCGTGCTCCATGCGGTGACGCTTGGCATGGGTGGTCGGGTTTACCGAGATCACGCCGCAGTATTCCGGCATGTGCTTGGCAAAGTCCGCGGTACCGATTTCGGT  >  NC_002947/5750355‑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: