Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,096,710 T→A 16.7% E86V (GAA→GTA)  lptB ← ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,096,7100TA16.7% 55.2 / 5.2 24E86V (GAA→GTA) lptBABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base A (2/2);  total (12/12)
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.83e-01

GAGGTTGTCCCGGATGTGGCTGATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTACCCGGCCCTGTTCGGCCTGGACCAGCCCGACAATCATGTAGAAGCA  >  NC_002947/1096565‑1096834
                                                                                                                                                 |                                                                                                                            
gAGGTTGTCCCGGATGTGGCTGATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTg                                                                                                                          >  1:567831/1‑150 (MQ=255)
    ttGTCCCGGATGTGGCTGATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGc                                                                                                                      >  1:289281/1‑150 (MQ=255)
         ccGGATGTGGCTGATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTa                                                                                                                 <  1:117356/150‑1 (MQ=255)
                   cTGATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCg                                                                                                       <  1:214866/150‑1 (MQ=255)
                     gATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGc                                                                                                     >  2:435326/1‑150 (MQ=255)
                               tCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCCGCGGCAGGTAGCCGATTCCGGCACGGGCGCgg                                                                                           >  2:355688/1‑150 (MQ=255)
                                 cTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCa                                                                                                    <  1:85251/138‑1 (MQ=255)
                                 cTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCa                                                                                                    >  2:85251/1‑138 (MQ=255)
                                  tGCAGCAGGCTTTCCAGCGCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGACTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCg                                                                                        >  2:164794/1‑150 (MQ=255)
                                      gcagGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTACCTGCCGCAGGAAGcct                                                                                                             <  1:576973/125‑2 (MQ=255)
                                      gcagGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTACCTGCCGCAGGAAGcct                                                                                                             >  2:576973/1‑124 (MQ=255)
                                            tttCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCt                                                                              >  2:248007/1‑150 (MQ=255)
                                                  gCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTTCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGc                                                                        <  1:454670/150‑1 (MQ=255)
                                                   ctctTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCt                                                                       <  1:123100/150‑1 (MQ=255)
                                                                ccGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTg                                                          <  2:567831/150‑1 (MQ=255)
                                                                            aGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCCGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTAcc                                              <  1:355688/150‑1 (MQ=255)
                                                                             ggTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTAccc                                             <  2:310829/150‑1 (MQ=255)
                                                                                 cTTGCGGGTCTCGCGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTACCCGGcc                                         >  2:473484/1‑150 (MQ=255)
                                                                                      gggTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTACCTGCCGCAGGAAGc                                                                                                               >  1:484008/1‑75 (MQ=255)
                                                                                      gggTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTACCTGCCGCAGGAAGc                                                                                                               <  2:484008/75‑1 (MQ=255)
                                                                                           tcGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTACCCGGCCCTGTTCGGcc                               <  1:572684/150‑1 (MQ=255)
                                                                                           tcGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTACCCGGCCCTGTTCGGcc                               >  2:572684/1‑150 (MQ=255)
                                                                                            cgcggATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTACCCGGCCCTGTTCGGCCt                              <  1:473484/147‑1 (MQ=255)
                                                                                                                        cgacAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATGCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGATGATCTCGTGGGTGAGCATTAGTACCCGGCCCTGTTCGGCCTGGACCAGCCCGACAATCAGGTAGAAgga  >  1:534555/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                            
GAGGTTGTCCCGGATGTGGCTGATGTGGAACTCCTGCAGCAGGCTTTCCAGCTCTTTGCGGCGGCCGTCCCGATCGAGGTCCTTGCGGGTCTCGAGGATGGCCATGATGTTGTCGGCGACCGACAGCTTGCGGAAGATCGAGGCTTCCTGCGGCAGGTAGCCGATTCCGGCACGGGCGCGGCCGTGCATGGGCTGGTGGCTGACATCGAGGTTGTCGATCAGTACCCGGCCCTGTTCGGCCTGGACCAGCCCGACAATCATGTAGAAGCA  >  NC_002947/1096565‑1096834

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