Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,997,194 Δ1 bp 14.0% intergenic (‑163/+37) PP_5252 ← / ← PP_5253 hypothetical protein/arylesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,997,1940C.14.0% 145.1 / 23.4 57intergenic (‑163/+37)PP_5252/PP_5253hypothetical protein/arylesterase
Reads supporting (aligned to +/- strand):  ref base C (17/27);  new base . (4/4);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

GGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCC  >  NC_002947/5997053‑5997342
                                                                                                                                              |                                                                                                                                                    
ggCAGGGGGACTGCGAAGCATGCCGGGTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGa                                                                                                                                               <  6:306716/149‑1 (MQ=255)
        gACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTGTGTAGGAGCGGGttt                                                                                                                                       >  2:14858/1‑149 (MQ=255)
          cTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTAc                                                                                                                                     <  4:73206/149‑1 (MQ=255)
              gAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTTAACCCGCTACTACAACGGATGTAGGAGCGGGTTTACCcgc                                                                                                                                 >  8:54366/1‑149 (MQ=255)
                 gCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                              <  7:175066/149‑1 (MQ=255)
                    tGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGt                                                                                                                           <  6:83774/149‑1 (MQ=255)
                     gCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCggg                                                                                                                                          >  4:204683/1‑133 (MQ=255)
                     gCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGACGTTGTAGGAGCggg                                                                                                                                          <  3:204683/133‑1 (MQ=255)
                     gCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGtt                                                                                                                          >  2:69062/1‑149 (MQ=255)
                               gCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                           <  6:255081/122‑1 (MQ=255)
                               gCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                           >  5:255081/1‑122 (MQ=255)
                               gCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAAt                                                                                                                 >  5:177178/1‑148 (MQ=255)
                                cAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGt                                                                                                                                                   >  1:232180/1‑113 (MQ=255)
                                cAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGt                                                                                                                                                   <  1:326778/113‑1 (MQ=255)
                                cAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGt                                                                                                                                                   >  2:326778/1‑113 (MQ=255)
                                cAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGt                                                                                                                                                   <  2:232180/113‑1 (MQ=255)
                                 aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCggg                                                                                                                                          >  1:201158/1‑121 (MQ=255)
                                 aaGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCggg                                                                                                                                          <  2:201158/121‑1 (MQ=255)
                                      cTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCAGACAACGTAGGTAGGAGCGGGTTTTCCCGCGAAATTTTGCCCCAATTAACggg                                                                                                          >  8:145714/1‑146 (MQ=37)
                                                aTGGAGGGGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGc                                                                                                <  7:54366/148‑1 (MQ=255)
                                                   gATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAg                                                                                             <  2:203891/148‑1 (MQ=255)
                                                    aTGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTGCAACAAAGTGTTTATGAGCGGGTTTTACCGCGAATGTTTACCACAATTAACGCTGAATGAACGCCAGc                                                                                            >  7:53054/1‑148 (MQ=255)
                                                      gtgCGGTGGGTTGATGTAGGTTAGCCCTGGTCAAGCTTTTTGGTGCGTGGGCGGGCCGCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTTGTTGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCaa                                                                                          <  4:318043/148‑1 (MQ=38)
                                                        gCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                       >  6:97752/1‑149 (MQ=255)
                                                         cGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                       <  1:69062/148‑1 (MQ=255)
                                                         cGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                       <  5:97752/148‑1 (MQ=255)
                                                         cGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTACCACAATTACCGCTGAAGGAACGCCAGCAAGTc                                                                                       >  2:83249/1‑148 (MQ=255)
                                                                ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTACTACAACGTATGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTAACCCTGAAGGAACCCCAGCAATTCCTCATTCa                                                                               >  6:144501/1‑149 (MQ=38)
                                                                ttGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAAACCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCAtt                                                                                 >  5:265101/1‑147 (MQ=255)
                                                                          ttAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAACTGCTgcgc                                                                     <  5:23871/149‑1 (MQ=255)
                                                                            aGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCAt                                                                   >  4:120951/1‑149 (MQ=255)
                                                                                      aaGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGc                                                         <  3:120951/149‑1 (MQ=255)
                                                                                           ttttGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCca                                                                                                                     <  7:333501/84‑1 (MQ=255)
                                                                                                ctgcCTGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgc                                               <  7:330728/149‑1 (MQ=255)
                                                                                                 tgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGt                                                                                        <  5:312253/107‑1 (MQ=255)
                                                                                                 tgcCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGt                                                                                        >  6:312253/1‑107 (MQ=255)
                                                                                                     tGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                           <  8:17331/52‑1 (MQ=255)
                                                                                                     tGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCgg                                                                                                                                           >  7:17331/1‑52 (MQ=255)
                                                                                                          cggccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCgg                                     <  5:144501/149‑1 (MQ=255)
                                                                                                            gccTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCa                                                                               <  3:244941/105‑1 (MQ=255)
                                                                                                               tctTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                              <  3:224998/55‑1 (MQ=255)
                                                                                                                 ttCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggc                                            <  4:272132/135‑1 (MQ=255)
                                                                                                                 ttCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTgcggc                                            >  3:272132/1‑135 (MQ=255)
                                                                                                                   cgcgGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATacac                              <  6:177178/147‑1 (MQ=255)
                                                                                                                       ggTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAg                        >  5:227284/1‑149 (MQ=255)
                                                                                                                           aaCCCGCTCCTACCACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTc                     <  7:145505/148‑1 (MQ=255)
                                                                                                                           aaCCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTc                     <  1:83249/148‑1 (MQ=255)
                                                                                                                             cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc                             <  3:66234/138‑1 (MQ=255)
                                                                                                                             cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc                             >  4:60975/1‑138 (MQ=255)
                                                                                                                             cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc                             <  3:60975/138‑1 (MQ=255)
                                                                                                                             cccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACAcc                             >  4:66234/1‑138 (MQ=255)
                                                                                                                             cccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgc                  <  6:265101/149‑1 (MQ=255)
                                                                                                                              ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCAtt                                                                                 <  3:181003/85‑1 (MQ=25)
                                                                                                                              ccGCTCCTACAGACGTTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCAtt                                                                                 >  4:181003/1‑85 (MQ=25)
                                                                                                                              ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg                 >  6:222530/1‑149 (MQ=255)
                                                                                                                              ccGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGgcg                 >  2:300216/1‑149 (MQ=255)
                                                                                                                                     tACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGaa          >  5:262131/1‑149 (MQ=255)
                                                                                                                                              cTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTcc  >  7:53635/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGCAGGGGGACTGCGAAGCATGCCGGTTGATGCAAGGCCTGGGCTTGAATGGATGTGCGGTGGGTTGAAGAAGGTTAGCCCTGGTCAAGCTTTTTGCTGCCTGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAACGTCTGTAGGAGCGGGTTTACCCGCGAATGTTTCCCACAATTACCGCTGAAGGAACGCCAGCAAGTCCTCATTCAGCTGCTGCGCATGAGTCACTGCAAAACCGTGCGGCGCCCCGGCATACACCTTCAGCTCGGCGCCACGAATCAGCTCC  >  NC_002947/5997053‑5997342

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