Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,480 2 bp→AC intergenic (+57/‑106) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,4800CA97.6% 119.8 / ‑6.8 42intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (17/24);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4810AC95.2% 117.3 / ‑5.6 42intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base C (17/23);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑‑TTT‑‑CGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG  >  NC_002947/1499339‑1499611
                                                                                                                                             ||                                                                                                                                             
cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCtt                                                                                                                                         <  1:394220/149‑1 (MQ=255)
 gACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCttt                                                                                                                                        <  6:60390/149‑1 (MQ=255)
  aCGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                                       <  4:178160/149‑1 (MQ=255)
    gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                                     >  8:254113/1‑149 (MQ=255)
     ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                   <  2:210677/149‑1 (MQ=255)
           aGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgcca                                                                                                                             <  2:304059/149‑4 (MQ=38)
             ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑ttcgc                                                                                                                        >  2:114419/1‑146 (MQ=255)
             ttCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCAAAGCCGgtgccatg                                                                                                                           <  6:363392/149‑6 (MQ=38)
                 aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑ttcgcggg                                                                                                                     >  7:111455/1‑142 (MQ=255)
                       cTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCAcgccc                                                                                                         <  2:128283/149‑2 (MQ=37)
                            gCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTcc                                                                                                    >  5:364679/1‑149 (MQ=25)
                             cccGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTccc                                                                                                   <  7:221284/149‑1 (MQ=25)
                               cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCcac                                                                                                 <  1:114419/149‑1 (MQ=25)
                                  gATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGa                                                                                              >  7:133774/1‑149 (MQ=21)
                                     gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                   >  6:274637/1‑117 (MQ=255)
                                     gCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                                   <  5:274637/117‑1 (MQ=255)
                                                   ccTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTTCCCGCGAAACACCCGCCCCCGTGCAAGGCACCCGCTATGCCCGTG‑‑TTAG‑‑CGGGCACCCCCGCTCCCCCAGACAGAATGCCATTAAGCg                                                                             >  6:221192/1‑149 (MQ=255)
                                                       aaCGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGTGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGc                                                                          <  8:133774/148‑1 (MQ=21)
                                                           aaGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATc                                                                     >  1:297915/1‑149 (MQ=21)
                                                            cgCCTTGGCAGTCGAGCAGGTCTGAGGGTGCCATACGGTGGGAGCGGGGGTGGCCGCGAAACAGGCGCCGGGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                    <  5:221192/148‑1 (MQ=18)
                                                            aGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                    <  2:297915/149‑1 (MQ=21)
                                                            aGCCTTGGCACCACAGCAGGACTGAGGCGGCCAAACGGTGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCg                                                                    <  3:294754/149‑1 (MQ=21)
                                                                ttGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTa                                                                 >  7:2431/1‑148 (MQ=21)
                                                                  ggCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCTTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACa                                                               <  6:364679/148‑1 (MQ=21)
                                                                                    aGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGtttttt                                              <  8:111455/147‑1 (MQ=21)
                                                                                      gCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCCGGTTTGGCGGAGGTTTTACCGGGCACCCCCCAACCCACACAATGATTGGCATTAAGCGGGCCTATAGCTACAGTGACGGCCCGttttttt                                             >  1:134977/1‑148 (MQ=255)
                                                                                        tGCCAAACGGTGGGAGCGGGCGTGGCCGCGCAACACGCGCCCCGGTGGATTGCACCCGCGTTTCCGGTG‑‑TTAG‑‑CGGGGACGCCCGCGCCCACAGACAGATTGGCAATTAACGGGGCGATAGGGAACATGACGGCCCGGTTTtttttggc                                        >  8:14115/1‑147 (MQ=255)
                                                                                                     gAGGGGGCGGGCCCGCTAAACAGGCCCCGGGGTGCATGGCACGGGCTTTGCGGTTG‑‑TTC‑‑GGGGCCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATc                           <  3:402084/149‑1 (MQ=255)
                                                                                                      cgCGGGCTTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  8:371006/148‑1 (MQ=25)
                                                                                                      aGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTTAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCa                          <  1:11476/149‑1 (MQ=25)
                                                                                                            cGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTc                    <  7:254113/149‑1 (MQ=25)
                                                                                                             gTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTcc                   <  8:22796/149‑1 (MQ=37)
                                                                                                                cccGCTAAACCGGCGCCGCGGTGCAGGGCACGGGCTTGCCGGGTG‑‑TGCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt                <  5:40226/149‑1 (MQ=255)
                                                                                                                 ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt                 >  8:6996/1‑147 (MQ=37)
                                                                                                                   gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAg              >  2:333454/1‑148 (MQ=37)
                                                                                                                        aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc         >  5:90910/1‑148 (MQ=37)
                                                                                                                        aCAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGc         >  5:131438/1‑148 (MQ=37)
                                                                                                                          aGGCGCCGCTGTGCTTGTCAACGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGCGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt      <  7:384999/149‑1 (MQ=37)
                                                                                                                           ggCGCCGCGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATc     <  7:6996/149‑1 (MQ=37)
                                                                                                                              gccgcGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg  >  5:389653/1‑149 (MQ=37)
                                                                                                                              gccgcGGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg  <  8:2431/149‑1 (MQ=37)
                                                                                                                                  cgGTGCATGGCACCGGCTTTGCCGGTG‑‑TTCG‑‑CGGGCACGCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGGTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAg  >  3:227702/1‑145 (MQ=37)
                                                                                                                                             ||                                                                                                                                             
CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTG‑‑TTT‑‑CGGGCACGCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTATCAAG  >  NC_002947/1499339‑1499611

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