Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,466,841 T→G 14.6% intergenic (‑16/+108) PP_3081 ← / ← PP_5560 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,466,8410TG14.6% 93.7 / 6.3 52intergenic (‑16/+108)PP_3081/PP_5560hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/26);  new base G (5/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

CTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGC  >  NC_002947/3466699‑3466987
                                                                                                                                              |                                                                                                                                                  
cTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTgcggcg                                                                                                                                              >  7:104655/1‑149 (MQ=255)
        aGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTgcggcg                                                                                                                                              <  6:214080/141‑1 (MQ=255)
        aGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTgcggcg                                                                                                                                              >  5:214080/1‑141 (MQ=255)
                  aGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAggg                                                                                                                             >  6:90835/1‑148 (MQ=255)
                   gCTCTGTTGCGGGGTCCGCTCTTCACCGGCACCGTATGCGGGGTCCTGGCCTACCTGTTCCAGTGCAATCACCCCGCAATCATTTAAAGACGCGTAACTGACCGACATCGCCCCTTGGCCCCCTGCGGCGGCCAACGCTGTGCAGGgca                                                                                                                           <  4:219280/149‑1 (MQ=255)
                       aGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGgcagc                                                                                                                         <  2:212889/147‑1 (MQ=255)
                       aGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGCCCAAAGCGGTGCAGGgcagca                                                                                                                        >  6:288603/1‑148 (MQ=255)
                        gTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCaaa                                                                                                                      <  5:90835/149‑1 (MQ=255)
                        gTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCaaa                                                                                                                      <  3:112083/149‑1 (MQ=255)
                         ttGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGGCGCCCGCGGCGGCCCACCCAGTGCAGGGCCGCCAag                                                                                                                     >  6:33986/1‑149 (MQ=255)
                           gCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCCCCTGCGGCGGCCCACGCTGTGCAGGGCAGCCAagcc                                                                                                                   >  4:224339/1‑148 (MQ=255)
                                gTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGaa                                                                                                              <  5:288603/149‑1 (MQ=255)
                                            ccGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGGCAACCCGGTGGAGGGGCGCAAAGCCGCGCACAGCCCgtgtgg                                                                                                  >  5:218224/1‑147 (MQ=255)
                                             cGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCAGCCGCGGCCCAACCAGTGCAGGGCAGCAAAACAGCGAACAGcccatgg                                                                                                    >  1:331638/1‑141 (MQ=255)
                                              ggCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCt                                                                                                  >  4:83005/1‑147 (MQ=255)
                                               gCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAg                                                                                               <  3:83005/149‑1 (MQ=255)
                                                   cGCATGCCGGGTCCTGGCCGACGTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGc                                                                                           <  4:281461/149‑1 (MQ=255)
                                                         ccGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCAGCGGCGGCCCACGCGCTTCCCGGCCGCAAAGCCGCGAACACCCAGTGTCTGAGGTGCATGGTc                                                                                     >  2:110503/1‑149 (MQ=255)
                                                                 tGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTg                                                                             >  6:94358/1‑149 (MQ=255)
                                                                        ccTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGc                                                                      <  3:50177/149‑1 (MQ=255)
                                                                            ttGCCGGGCAATGACGCCGCCATCCGTCAAAGACGAGCAACTGACCCACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCa                                                                  <  1:366724/149‑1 (MQ=255)
                                                                            tggCAGGGCAATGCCGCCGCCATCAGTCAAATACCAGCAACTGACCGACAGCCCCCCTTGGCCGCCGGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCa                                                                  <  6:44723/147‑1 (MQ=255)
                                                                                   gCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCa                                                           <  5:94358/149‑1 (MQ=255)
                                                                                   gCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCa                                                           <  4:290260/149‑1 (MQ=255)
                                                                                     aaTGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCAGCGGCGGCCCACCCTGGGCAGGGCCGCCAAGCAGCGAACAGCCAGTGGCGGAGGTGCAGGGTCAATCCTTGGCAGGGCGTCACGGCACAgg                                                         >  8:91382/1‑148 (MQ=255)
                                                                                           gccgccATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGc                                                                      <  6:186752/130‑1 (MQ=255)
                                                                                             cgccATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCa                                                 <  8:104655/149‑1 (MQ=255)
                                                                                                     gTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCGGCGGCGGCCCACCCAGTGCAGGGCCCCAAAACAGCGAACACCCCCTGTCGGAGGTGGAAGGTGAATACTTTGCACGGCGTAACGGCGCACTCAAGCCCAGCTCAGca                                         >  6:356086/1‑148 (MQ=255)
                                                                                                              gAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGGAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGGCACGGCTCAGTCTAGCCCACGTCTGCGTCAGCgg                                  >  8:294592/1‑147 (MQ=255)
                                                                                                              gAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCt                                <  3:199983/149‑1 (MQ=255)
                                                                                                               aGGAACTGACCTACATCGCGCCGTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTg                               <  6:218224/149‑1 (MQ=255)
                                                                                                                 gAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAg                                                                         >  5:76925/1‑105 (MQ=255)
                                                                                                                 gAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAg                                                                         >  5:82797/1‑105 (MQ=255)
                                                                                                                 gAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAg                                                                         <  6:82797/105‑1 (MQ=255)
                                                                                                                 gAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAg                                                                         <  6:76925/105‑1 (MQ=255)
                                                                                                                 gAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTc                             <  5:8157/149‑1 (MQ=255)
                                                                                                                  aaCTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTcc                            <  1:301223/149‑1 (MQ=255)
                                                                                                                  aaCTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTcc                            <  5:356086/149‑1 (MQ=255)
                                                                                                                       aCCGACATCGCGCCTTGGCCCCCAGCGGCGGCCCACGCGGTGCAGGGCAGCCAAGCAGCGCACCGCCAGGGGCGGAGGTGGATGGGCAATCCGTGGCCGGGCGTCACGGCTCACTCTAGCCCAGGGCTTCGTCAGCGGCTGGCCTGGTc                       >  6:310036/1‑149 (MQ=255)
                                                                                                                           aCATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCg                   >  5:308512/1‑149 (MQ=255)
                                                                                                                            cATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCgg                  <  6:308512/149‑1 (MQ=255)
                                                                                                                              tCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGACCAGCCAGTGGCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCa                <  7:91382/149‑1 (MQ=255)
                                                                                                                                  gcCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTg                          >  2:66917/1‑135 (MQ=255)
                                                                                                                                  gcCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTg                          <  1:66917/135‑1 (MQ=255)
                                                                                                                                  gcCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGc            <  5:33986/149‑1 (MQ=255)
                                                                                                                                   cgttGGCCGCCGCAGGCGGCCAAGGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGc            >  6:116089/3‑148 (MQ=255)
                                                                                                                                    gttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCt          >  3:347397/2‑149 (MQ=255)
                                                                                                                                    gttGGCCGCCGCAGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCt          >  3:347390/2‑149 (MQ=255)
                                                                                                                                       ggCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCg       >  4:357565/1‑149 (MQ=255)
                                                                                                                                       ggCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCg       >  7:358598/1‑149 (MQ=255)
                                                                                                                                         ccgccTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCt     <  5:39075/149‑1 (MQ=255)
                                                                                                                                            ccTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGc  <  8:358598/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
CTGCGGGCAGGTTGTCCCAGCTCAGTTGCGGGGTCCGCTCTTCACCGGCACCGCATGCCGGGTCCTGGCCGACCTGTTGCAGGGCAATGACGCCGCCATCAGTGAAAGACGAGGAACTGACCGACATCGCGCCTTGGCCGCCTGCGGCGGCCAACGCTGTGCAGGGCAGCAAAGCAGCGAACAGCCAGTGTCTGAGGTGCATGGTCAATCCTTGGCAGGGCGTCACGGCTCAGTCTAGCCCAGGTCTGCTTCAGCGGCTGTCCTGGTCCTCGGCAGTGCCTCCGCTGGC  >  NC_002947/3466699‑3466987

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