Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,578,431 T→A 7.1% intergenic (+108/‑216) PP_3158 → / → PP_5589 hypothetical protein/rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,578,4310TA 7.1% 168.5 / 3.1 56intergenic (+108/‑216)PP_3158/PP_5589hypothetical protein/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (20/32);  new base A (2/2);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATG  >  NC_002947/3578284‑3578575
                                                                                                                                                    |                                                                                                                                                
gcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCGCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                 >  2:148075/1‑149 (MQ=255)
 cAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTg                                                                                                                                                <  8:8876/149‑1 (MQ=255)
  agagCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                               <  5:125658/149‑1 (MQ=255)
             gtgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTGCAGc                                                                                                                                     <  8:277905/147‑1 (MQ=255)
              tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGGATTTTTTGCAGTTGGAGCg                                                                                                                                    >  4:577/1‑148 (MQ=255)
              tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCAAAAAAGAGTTTCGTCGGGCTGCAACAGGGATTTTTTGCAGtttgagc                                                                                                                                     >  1:147699/1‑142 (MQ=255)
               gatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAg                                                                                                                                  >  7:282174/1‑149 (MQ=255)
                atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCGCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAg                                                                                                                                  >  1:284992/1‑149 (MQ=255)
                atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                   <  4:293089/147‑1 (MQ=255)
                atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAg                                                                                                                                  <  7:183895/148‑1 (MQ=255)
                atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGACCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAATTTTGACAGTTGCAGCGa                                                                                                                                   >  6:434844/1‑147 (MQ=255)
                   aaCTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTa                                                                                                                              <  4:217303/149‑1 (MQ=255)
                    aCTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCTTCGGGCTGCAACTTGCATTTTTGCCAGTTGCAGCGAGTGTAt                                                                                                                             <  6:142839/149‑1 (MQ=255)
                            ccgtggCAGTAAGAGGCCTCTGCCGCAGGCCGTTTGCGAAATGCCCTTACTTGAATGTGTATTGTGTTGATTAGTGTTCTTTAGCGATTCAACGTTTCGTCTGGCTGAAACTGTCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCt                                                                                                                     <  8:347868/144‑1 (MQ=255)
                            ccGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGc                                                                                                                       <  3:577/147‑1 (MQ=255)
                               gTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAg                                                                                                                   >  3:178231/1‑148 (MQ=255)
                                  cAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAAtt                                                                                                               <  1:148075/149‑1 (MQ=255)
                                    gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCa                                                                                                             <  7:161052/149‑1 (MQ=255)
                                    gTAAGACGCCTGTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATTAGGATTGAGTTGATGAGTGTTCATTAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCa                                                                                                             <  6:15208/149‑1 (MQ=255)
                                      atgagGCCTCTTCCGCACGCCCTTTCCTAAATGCCCATACATGAATGAGGATTGAGTTGATTAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                            <  4:258341/146‑1 (MQ=255)
                                      aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                            <  3:38548/148‑1 (MQ=255)
                                       agagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAg                                                                                                           <  8:316699/148‑1 (MQ=255)
                                        gacgCCGCTGCCGCACGCCCGTGGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAggg                                                                                                         <  6:300706/146‑1 (MQ=255)
                                           gCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGGATTTTTTGCAGTTGCAGCGAGTGTATTAAAGCCGAGAATTAAAGGGaaa                                                                                                      >  8:55974/1‑149 (MQ=255)
                                           gCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGaaa                                                                                                      <  2:284992/149‑1 (MQ=255)
                                                        cGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAAATGGCATTTTTTGCAGTTGCAGCGAGTGTATGACAGCCTATAATTTAAGGGAAATgttgttgtat                                                                                           >  7:144374/1‑144 (MQ=255)
                                                           ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAgg                                                                                       >  5:153949/1‑149 (MQ=255)
                                                             tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGGTGCAACAGGCAATTTTTGCAGTTGCGACGGAGGTATTTCAGCCTAGAATTAAAGGGAAAATGTGTTATGTAAGGt                                                                                      >  6:40355/1‑148 (MQ=255)
                                                                gCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCAt                                                                                   <  6:153949/148‑1 (MQ=255)
                                                                  gAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAAGCATTTTTTCCAGTTGAAGCGAGTGTATGACAACCTACAATTCAAGTGAAAAGGTGTTATCTAAGGTCATTa                                                                                 >  4:334456/1‑148 (MQ=255)
                                                                    aaTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATc                                                                              <  1:14266/149‑1 (MQ=255)
                                                                      tGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                            <  3:334456/149‑1 (MQ=255)
                                                                               atgaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                   <  7:260112/149‑1 (MQ=255)
                                                                               atgaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTGCAGCGAGTGTATTAAAGCCTAGAATTCAAAGGAAATTGTTTTATATAAGGTCATTAATAGTGAAACGTTc                                                                   >  2:122030/1‑149 (MQ=255)
                                                                                tgaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCg                                                                  <  4:178231/149‑1 (MQ=255)
                                                                                  aatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGc                                                                 <  1:436489/148‑1 (MQ=255)
                                                                                   atgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTACGCaa                                                               <  8:282174/149‑1 (MQ=255)
                                                                                   atgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTAAAAGGAAAATGTTTTGTGTAATGTCATTTATAGTGAAAAGTTCGCaa                                                               >  6:387225/1‑149 (MQ=255)
                                                                                      aGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGc                                                             <  1:122030/148‑1 (MQ=255)
                                                                                       ggATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGt                                                           <  7:286731/149‑1 (MQ=255)
                                                                                             aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTTCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAATGGTGTTGTCTAAGGTAATTAATCGTGAAACGTTAGCAAGCGTaaga                                                       >  8:106420/1‑144 (MQ=255)
                                                                                             aGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGt                                                           <  7:106420/143‑1 (MQ=255)
                                                                                              gTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGcc                                                    >  5:184968/1‑149 (MQ=255)
                                                                                                  atgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGGATTTTTGCCAGTTGCAGCGAGTGTATGAAAGCCTAGAATTCAAGGGAAATGGTGTTATCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTAc                                                >  5:194383/1‑149 (MQ=255)
                                                                                                     aGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGGATTTTTTGCAGTTTCAGCGAGTGTATTACAGCCTAGAATTAAAGGGAAATGGTGTTCTCTAAAGTCAATTATAGTGAAAAGTTCGCAAGCGTACACCCCCTACGgc                                             >  8:48512/1‑149 (MQ=255)
                                                                                                        gtTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGGGTTTTTGGCAGTTGGGGGGGGTGGTTGGCAGCGTAGAATTGAAGGGGAATGGGGTTGTGTAAGGGGATTAATGGTGGAACGTTGGGAAGAGTAGACGGCCGACGgcgcg                                          >  8:72531/1‑148 (MQ=255)
                                                                                                             tGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTg                                     <  5:40355/149‑1 (MQ=255)
                                                                                                                gCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg                                                                                                          <  3:32751/77‑1 (MQ=255)
                                                                                                                     tCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                   <  6:202894/143‑1 (MQ=255)
                                                                                                                     tCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                   >  5:202894/1‑143 (MQ=255)
                                                                                                                         cGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCa                         <  5:387225/149‑1 (MQ=255)
                                                                                                                                      tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt             >  8:382159/1‑148 (MQ=255)
                                                                                                                                       gCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGtt                                                                    <  5:127104/92‑1 (MQ=38)
                                                                                                                                       gCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGtt                                                                    >  6:127104/1‑92 (MQ=38)
                                                                                                                                         aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg         <  6:184968/149‑1 (MQ=255)
                                                                                                                                            tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa       >  4:284398/1‑148 (MQ=255)
                                                                                                                                                aTTTTTGCCAGTTGCAGCGAGTGTTTGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTAATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATg  <  3:170356/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATG  >  NC_002947/3578284‑3578575

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: