Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,868,475 A→G 7.3% I30T (ATC→ACC)  lgt ← prolipoprotein diacylglyceryl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,868,4750AG 7.3% 157.1 / 3.2 55I30T (ATC→ACC) lgtprolipoprotein diacylglyceryl transferase
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base G (2/2);  total (26/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 = 5.72e-01

TGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGAACAA  >  NC_002947/5868328‑5868612
                                                                                                                                                   |                                                                                                                                         
tGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCgccgat                                                                                                                                          <  5:89393/149‑1 (MQ=255)
  cAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCgccgatgc                                                                                                                                        <  5:255456/149‑1 (MQ=255)
       cGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATc                                                                                                                                   <  7:185932/149‑1 (MQ=255)
        gTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTCCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCgccgatgccgat                                                                                                                                    >  3:100153/1‑147 (MQ=255)
         tAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAg                                                                                                                                 >  5:240001/1‑149 (MQ=255)
           gAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACCCCCGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGt                                                                                                                               <  5:126613/149‑1 (MQ=255)
           gAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAgg                                                                                                                                >  2:240861/1‑148 (MQ=255)
                 cacaTAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACAt                                                                                                                           >  1:163800/1‑147 (MQ=255)
                 cacaTAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACAt                                                                                                                           >  1:169488/1‑147 (MQ=255)
                    aTAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGc                                                                                                                      <  2:13816/149‑1 (MQ=255)
                                 ccgccgACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGAt                                                                                                         >  7:97537/1‑149 (MQ=255)
                                    ccgACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGAtttt                                                                                                      <  5:68027/149‑1 (MQ=255)
                                         attCACCTCGATCGCCAGCCAGACCCGGAGGTCGGACAGCTTCGCTCGGTTGCAGGTGGGGTCGAAGGGTTTGAGCCGGCGTGAGGCCAGCTGCCAGGGACCGCCGATGCGGATCAGGTACATCATGCCGTACCAATGGATTTTCAggg                                                                                                 <  4:219668/147‑3 (MQ=255)
                                           tCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGc                                                                                               >  7:47023/1‑149 (MQ=255)
                                            cACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTCCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGcc                                                                                              <  4:100153/149‑1 (MQ=255)
                                              cccccATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCg                                                                                            <  1:122351/149‑1 (MQ=255)
                                                 ccATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATg                                                                                         <  1:240861/149‑1 (MQ=255)
                                                       aCAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   >  4:15608/1‑149 (MQ=255)
                                                         aGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   >  5:251103/1‑147 (MQ=255)
                                                         aGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   <  6:251103/147‑1 (MQ=255)
                                                                 cacCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGt                                                                         <  8:97537/149‑1 (MQ=255)
                                                                 cacCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGt                                                                         <  3:15608/149‑1 (MQ=255)
                                                                  acCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTc                                                                        <  2:36151/149‑1 (MQ=255)
                                                                   ccAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCTTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCt                                                                       <  5:153428/149‑1 (MQ=255)
                                                                    cAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACg                                                                             >  1:17296/1‑142 (MQ=255)
                                                                    cAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACg                                                                             <  2:17296/142‑1 (MQ=255)
                                                                    cAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTa                                                                      >  5:332297/1‑149 (MQ=255)
                                                                     aGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTAt                                                                     <  6:240001/149‑1 (MQ=255)
                                                                          ggACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGc                                                                 <  8:47023/148‑1 (MQ=255)
                                                                          ggACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCg                                                                >  4:293140/1‑149 (MQ=255)
                                                                             cAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCggg                                                              >  1:114178/1‑148 (MQ=255)
                                                                                 ttCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAgg                                                         >  2:262505/1‑149 (MQ=255)
                                                                                   ctcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAgg                                                         <  3:293140/147‑1 (MQ=255)
                                                                                   ctcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGgc                                                        <  3:339466/148‑1 (MQ=255)
                                                                                    tcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATg                                                                                                            >  8:320872/1‑95 (MQ=255)
                                                                                    tcGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATg                                                                                                            <  7:320872/95‑1 (MQ=255)
                                                                                            ccAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAAc                                              >  3:258044/1‑149 (MQ=255)
                                                                                                        gAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTc                                           >  2:250087/1‑140 (MQ=255)
                                                                                                        gAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTc                                           <  1:250087/140‑1 (MQ=255)
                                                                                                         aaGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   <  3:151886/99‑1 (MQ=255)
                                                                                                         aaGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAg                                                                                   >  4:151886/1‑99 (MQ=255)
                                                                                                                  aGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATtcaaatcaaa                         >  5:127425/1‑148 (MQ=255)
                                                                                                                   gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATg                                                                                         >  6:110069/1‑83 (MQ=255)
                                                                                                                   gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATg                                                                                         <  5:110069/83‑1 (MQ=255)
                                                                                                                   gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCaaaaa                       <  2:163800/149‑1 (MQ=255)
                                                                                                                   gCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCaaaaa                       <  2:169488/149‑1 (MQ=255)
                                                                                                                    ccGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAg                      >  4:282901/1‑149 (MQ=255)
                                                                                                                           gAGGCCAGCAGCCAGGCACCGCCGATGCCGTTCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGc               >  1:33265/1‑149 (MQ=255)
                                                                                                                             ggCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCa             <  4:258044/149‑1 (MQ=255)
                                                                                                                              gCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGa                  <  3:295098/143‑1 (MQ=255)
                                                                                                                              gCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGa                  >  4:295098/1‑143 (MQ=255)
                                                                                                                                   cagcCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAg       <  2:114178/149‑1 (MQ=255)
                                                                                                                                     gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                >  6:67610/5‑106 (MQ=255)
                                                                                                                                     gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                <  7:116694/102‑1 (MQ=255)
                                                                                                                                     gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                <  5:67610/102‑1 (MQ=255)
                                                                                                                                     gatcggCATCGGCGGTGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGa                                                >  8:116694/5‑106 (MQ=255)
                                                                                                                                      ccAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGaa     >  3:220025/1‑148 (MQ=255)
                                                                                                                                        aGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGACCTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGaacaa  >  8:328047/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
TGCAGGTCGTAGAACAGCACATAGCCCAGGCGCCCGCCGACAATCACCCCCATCGACAGCCAGAACACCAGGTCGGACAGCTTCTCGCGGTTCCAGGTGGGGTCGAAGCGGTTCAGCCGGCGTGAGGCCAGCAGCCAGGCACCGCCGATGCCGATCAGGTACATCAGGCCGTACCAATGGATTTTCAGCGGCCCGATGGCCAAGGCCACGGGGTCTATCTGCGGGTAAGGCAGCATGGAACTTCCTCTCGATTCAAATCAAAAAGCTGAGGCCAACACAGAACAA  >  NC_002947/5868328‑5868612

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