Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,897,783 C→G 8.7% intergenic (‑183/+51) PP_5172 ← / ← PP_5173 hypothetical protein/RND family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,897,7830CG 8.7% 168.1 / 6.0 57intergenic (‑183/+51)PP_5172/PP_5173hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base C (21/30);  new base G (2/3);  total (23/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

AAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC  >  NC_002947/5897641‑5897926
                                                                                                                                              |                                                                                                                                                
aaGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                            <  8:168273/149‑1 (MQ=255)
         tGCGCCTGGCTTGAGGGCACTTTTCGTGAGTCGGGATTATTCTTAAAAATACTGTGTAAGAATTTTTATATAACCATGTATGACCACACGCCTCTGGTAGGAGCGGGTTAACCCGCGATTACGGTGGTCGATGCACCATGGCATTcgc                                                                                                                                    <  6:420650/148‑1 (MQ=255)
          gcgcATGGATTGAGGGCAATAGTAGTGAGGCGTGATTATTCTTAAAAATACTCTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgcg                                                                                                                                   <  2:239766/148‑1 (MQ=255)
           cgcCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                  <  5:314377/148‑1 (MQ=255)
           cgcCTGGATTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                                  <  6:150067/148‑1 (MQ=255)
                     gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCg                                                                                                                         <  5:163688/147‑1 (MQ=255)
                        ggCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgctg                                                                                                                                  <  3:124980/135‑3 (MQ=25)
                         gaaaTAGAAGTTAGGCGTGATTATTCTTAAAAATAGTTTTTAATAATTTTTATATAACCATTTAAGAGACAACGCGTCTGGTTGGAGCGGGTTTACCCGGGAATACGGGGGTCGATGCACCATGGCATTCGCGGGTAAACCCGCTCCTa                                                                                                                   <  6:407919/147‑1 (MQ=37)
                             tagtagTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAg                                                                                                               >  2:88892/1‑149 (MQ=255)
                                   tGAGGGGGTATTATTATTAAAAATAATGTTTAAGAATTTTTATATAACAATTTAAGAGAACACGCGTGTGGTAGGAGGGGGTTTACCGGCGAATACGGTGGTCGATGCACCATCGAATTCGCGGGTAAACCCGCTCCTACCAGGGGcc                                                                                                          <  2:328303/148‑1 (MQ=37)
                                    gAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCACGGGCCa                                                                                                         >  6:146583/1‑148 (MQ=255)
                                         gggATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCTGGTTTACCAGCGAAAACGGTGGTCGATGCAACATCTCAATCGCAGGTAAAACAGCTCCTACCAGGGGCCATATggg                                                                                                   >  8:129959/1‑149 (MQ=255)
                                               tttCTTAAAAATACTTTTTAATATTTTTTATATAACAATTAAAGACAACACGCCTGTGTTAGTAGGGGTTTTACCCGCTAATAGGGTGGTGGATGCACCATCGCATTCGGGGTAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTc                                                                                              <  2:169111/147‑1 (MQ=255)
                                                 tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAg                                                                                            <  1:88892/148‑1 (MQ=255)
                                                 tCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTAAAAGCGAATAAGGTGGTGGATGCACCATAGCATTAGCGGGTAAACCCGCTCATAACAGGGGCCATATGGGCCTTCAgg                                                                                           >  7:141642/1‑149 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                    >  6:147303/1‑106 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                                    <  5:147303/106‑1 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                                 >  1:140124/1‑109 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCggg                                                                                                                                 <  2:140124/109‑1 (MQ=255)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAAGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTAc                                                                                                                  <  6:158281/124‑1 (MQ=17)
                                                   ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAAGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTAc                                                                                                                  >  5:158281/1‑124 (MQ=17)
                                                     aaaaaTACTGTTAAAGAATTTTTAAAAAACCATAAAAGACAACACGCCTCTGGTAGGAGCGGGTTAACCCGCCAATACGGTGGGCGATGCACCATCGCATTCGCGGGAAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGc                                                                                          <  5:146239/146‑1 (MQ=38)
                                                      aaaaTACTGTTTAAGAATTTTTATATAACCATTAAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCGGGCgtgt                                                                                      <  3:153607/149‑1 (MQ=255)
                                                          tACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCTTTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCt                                                                                  <  6:183203/149‑1 (MQ=255)
                                                               tttAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGc                                                                             >  5:12818/1‑149 (MQ=255)
                                                               tttAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCAGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGTGACTTCATGAGTGTTTCTGCAGc                                                                             >  6:289957/1‑149 (MQ=38)
                                                                    gAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTcgcg                                                                        >  6:79119/1‑149 (MQ=255)
                                                                       tttttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGAt                                                                     <  7:129959/149‑1 (MQ=255)
                                                                        ttttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGAt                                                                     <  3:285837/148‑1 (MQ=255)
                                                                             tataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                               <  2:287883/149‑1 (MQ=255)
                                                                                       tAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTAACAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCAc                                                      <  5:289957/148‑1 (MQ=255)
                                                                                         acacaaCACGCCTCTTGTAGGCGCGGTTTTACCCGCTAATACGGTGGTCTATGAACCATCGCATTCGCGTTTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                                   <  7:259539/145‑1 (MQ=38)
                                                                                          gataaCTCGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCACCATCTCATTCGCGGTTAAACCCGCTCCTACCAGGGGc                                                                                                           <  2:303157/89‑1 (MQ=21)
                                                                                          gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGc                                                                                                           >  1:303157/1‑92 (MQ=255)
                                                                                           agaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                                 <  1:235253/149‑1 (MQ=255)
                                                                                                        ggTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTTAACCCGCTGCTACCACGGGGCATATGGGCCTTCAGGCGTGTTTCTGCGGCTGGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGaa                                   >  5:201398/1‑149 (MQ=38)
                                                                                                               gCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                            >  3:198469/1‑149 (MQ=255)
                                                                                                                 gggTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                                 <  3:191940/111‑1 (MQ=255)
                                                                                                                 gggTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                                 >  4:191940/1‑111 (MQ=255)
                                                                                                                  ggTTTACCGGCGAATACGGTGGTCGATGCACCATCGCCTTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTAc                                                              <  5:332320/113‑1 (MQ=255)
                                                                                                                  ggTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCATGCGTGGTTCTGCGGCTCGCGGCTCTTGtcc                                                              >  6:332320/1‑111 (MQ=255)
                                                                                                                   gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGGGTAGAACAGCAGCGt                          >  1:140231/1‑148 (MQ=39)
                                                                                                                         ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                            >  3:142419/1‑139 (MQ=25)
                                                                                                                         ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACagcagc                            <  4:142419/139‑1 (MQ=25)
                                                                                                                         ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                      <  5:379614/129‑1 (MQ=255)
                                                                                                                         ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                      >  6:379614/1‑129 (MQ=255)
                                                                                                                         ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg                  >  8:391942/1‑149 (MQ=255)
                                                                                                                          cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                 >  1:263191/1‑149 (MQ=255)
                                                                                                                          cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt                 <  5:79119/149‑1 (MQ=255)
                                                                                                                           gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                               >  5:156384/1‑103 (MQ=255)
                                                                                                                           gcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTa                                                               <  6:156384/103‑1 (MQ=255)
                                                                                                                             gAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGc              <  7:391942/149‑1 (MQ=255)
                                                                                                                               aTAGGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGTTGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCa            <  1:117571/149‑1 (MQ=255)
                                                                                                                                tACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAc           >  6:298759/1‑149 (MQ=255)
                                                                                                                                  cGGTGGTCTATCCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACga         <  6:12818/149‑1 (MQ=255)
                                                                                                                                       gtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATc    >  6:236328/1‑149 (MQ=255)
                                                                                                                                         cGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATccc  <  5:146583/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                
AAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATGATCCC  >  NC_002947/5897641‑5897926

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: