Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,210,758 C→A 8.3% intergenic (+35/+34) PP_3692 → / ← PP_5621 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,210,7580CA 8.3% 164.8 / 2.1 56intergenic (+35/+34)PP_3692/PP_5621hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/23);  new base A (3/2);  total (31/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

GCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCTTC  >  NC_002947/4210615‑4210905
                                                                                                                                               |                                                                                                                                                   
gCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCgg                                                                                                                                                >  5:270332/1‑149 (MQ=255)
     cccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGATAAATGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGAATAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAAACGCGTCAAACTCGCTTCTg                                                                                                                                           >  4:234683/1‑149 (MQ=255)
      ccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                          >  5:232874/1‑149 (MQ=255)
      ccAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTgg                                                                                                                                          >  5:38280/1‑149 (MQ=255)
            ggaggaGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTg                                                                                                                                           >  5:254474/1‑142 (MQ=255)
            ggaggaGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGGTTGACTCGGTTCTg                                                                                                                                           <  6:254474/142‑1 (MQ=255)
             gaggagCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGc                                                                                                                                    >  7:84465/1‑148 (MQ=255)
                    ttCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCt                                                                                                                                            <  5:211641/133‑1 (MQ=255)
                    ttCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCt                                                                                                                                            >  6:211641/1‑133 (MQ=255)
                      cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCACCCTCGGTTCTGGCTTAGCCaaa                                                                                                                                <  7:47375/143‑1 (MQ=255)
                      cccAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCaaa                                                                                                                                >  8:47375/1‑143 (MQ=255)
                        cagcagTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgcagc                                                                                                                        <  3:107041/149‑1 (MQ=255)
                                  ttGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCCCAGAACCGAGGTTGACTCGGTTCTgagcccg                                                                                                                                     >  2:111835/1‑120 (MQ=38)
                                  ttGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCCCAGAACCGAGGTTGACTCGGTTCTgagcccg                                                                                                                                     <  1:111835/126‑7 (MQ=38)
                                        cTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACg                                                                                                         >  5:313530/1‑148 (MQ=255)
                                              gTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACATCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAgcggc                                                                                                  >  1:146672/1‑149 (MQ=255)
                                                 ttCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGAt                                                                                               >  2:236148/1‑149 (MQ=255)
                                                  tCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGAt                                                                                               <  6:232874/148‑1 (MQ=255)
                                                    gAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTcc                                                                                            >  2:268596/1‑149 (MQ=255)
                                                     aGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCg                                                                                           <  4:78222/149‑1 (MQ=255)
                                                        cTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGAt                                                                                         >  1:20586/1‑148 (MQ=255)
                                                         tGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGAtt                                                                                        >  4:289346/1‑148 (MQ=255)
                                                           aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTa                                                                                       <  6:313530/147‑1 (MQ=255)
                                                           aCGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTa                                                                                       <  1:268596/147‑1 (MQ=255)
                                                                    cactgcaTAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCaa                                                                              <  4:151249/147‑1 (MQ=255)
                                                                                 gATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgca                                                                                                                          <  4:107553/90‑1 (MQ=255)
                                                                                 gATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGgca                                                                                                                          >  3:107553/1‑90 (MQ=255)
                                                                                      aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCg                                                                                                                                                 >  2:25354/1‑62 (MQ=255)
                                                                                      aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGAt                                                                       >  6:199267/1‑136 (MQ=255)
                                                                                      aCGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGAt                                                                       <  5:199267/136‑1 (MQ=255)
                                                                                      aCGAATTCTCAGCACGCAAGTAGAAGCATACCGTGAGCGGGCTCAGAACCGAGTCAACCTCg                                                                                                                                                 <  1:25354/62‑1 (MQ=255)
                                                                                       cGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                            >  1:126359/1‑146 (MQ=255)
                                                                                       cGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCgg                                                            <  2:126359/146‑1 (MQ=255)
                                                                                        gAATTCTCAGAACGCAAGTAGAAGCAAACGGTGAGCGGGCTCAGAACCGAGTAAACCTCGGTTCTGCCTTAGCCAAACTGGCAGCATAGTCACCAAAACGTAGAGGCGATTCAGATTACCATCCCAATACCGATCAAGCTACCGGGGAt                                                        >  2:274354/1‑149 (MQ=255)
                                                                                         aaTTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGTGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCTATg                                                       <  8:84465/149‑1 (MQ=255)
                                                                                          aTTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGa                                                      <  6:38280/149‑1 (MQ=255)
                                                                                           ttCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAg                                                     <  3:289346/149‑1 (MQ=255)
                                                                                           ttCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAg                                                     <  2:146672/149‑1 (MQ=255)
                                                                                            tctcAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGTACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGc                                                    >  1:272596/1‑149 (MQ=255)
                                                                                             ctcAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCt                                                   >  6:260479/1‑149 (MQ=255)
                                                                                              tcAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTc                                                  >  2:98167/1‑149 (MQ=255)
                                                                                               cAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTcc                                                 >  5:61931/1‑149 (MQ=255)
                                                                                                  cACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCACACTGGCAGCATACTCAGCAAAACGTAGCGGCGATTCCGAGTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCg                                               >  8:181798/1‑148 (MQ=255)
                                                                                                  cACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCATCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCg                                               >  8:177132/1‑148 (MQ=255)
                                                                                                    cGCAAGTAGAAGCATACTGTGAGCTGGCTCAGAACCTAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGAATACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGAt                                            <  3:234683/149‑1 (MQ=255)
                                                                                                     gCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATg                                           <  2:20586/149‑1 (MQ=255)
                                                                                                          tAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGt                                      >  5:229869/1‑149 (MQ=255)
                                                                                                                cATACGGTGAGCGGGCTCAGTACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTc                                <  2:272596/149‑1 (MQ=255)
                                                                                                                 aTACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTcc                               <  1:98167/149‑1 (MQ=255)
                                                                                                                 aTACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTcc                               <  1:236148/149‑1 (MQ=255)
                                                                                                                  tACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTccc                                                <  8:9632/131‑1 (MQ=255)
                                                                                                                  tACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTccc                                                >  7:9632/1‑131 (MQ=255)
                                                                                                                                         aGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTg         >  5:89850/1‑147 (MQ=255)
                                                                                                                                           tCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGc     >  8:172361/1‑149 (MQ=255)
                                                                                                                                             aaCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACAGGCGATGAGCTACCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCtt   >  5:163317/1‑149 (MQ=255)
                                                                                                                                              aCCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCTTc  <  7:172361/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCATGCCCAGTCGGAGGAGCTTCCCAGCAGTGCCTTGGAGCTAAAGGTATTCGAGGCTGACGTCACTGCACTGCATAACTGGATTAACGAATTCTCAGCACGCAAGTAGAAGCATACGGTGAGCGGGCTCAGAACCGAGTCAACCTCGGTTCTGGCTTAGCCAAACTGGCAGCATACTCACCAAAACGTAGCGGCGATTCCGATTACCATCCCAATACCGATCAAGCTACCGGCGATGAGCTCCCGGATGGCTGTCGCCTCCATGAAGTGACGCAGCCGACTTGCAGCTTC  >  NC_002947/4210615‑4210905

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