Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,604 G→T 7.3% intergenic (+59/‑237) PP_2661 → / → PP_2662 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,048,6040GT 7.3% 156.2 / 2.1 53intergenic (+59/‑237)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/28);  new base T (2/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

TGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAA  >  NC_002947/3048457‑3048748
                                                                                                                                                   |                                                                                                                                                
tggtggATGTGCAGTCACGGCTGGACTGTCAGCAACCGCACCGTTATGAAATGCCCGCTTGCAAGCTTTAGCAGGCACTGGCTTCCTGAAGCCAGGAGAATCCGCTCAGATAATTTAATTGGGGCTGCTTTGCAGCCCTTAGCGGGCGa                                                                                                                                                 <  2:61197/149‑1 (MQ=255)
                aCGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGc                                                                                                                                 <  6:181416/149‑1 (MQ=255)
                  ggCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCaa                                                                                                                               <  1:169026/149‑1 (MQ=255)
                       gcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggc                                                                                                                            <  6:181323/147‑1 (MQ=255)
                            gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                          >  4:302084/1‑144 (MQ=255)
                            gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAggcgg                                                                                                                          <  3:302084/144‑1 (MQ=255)
                            gcagcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCg                                                                                                                     >  5:115347/1‑149 (MQ=255)
                              agcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCaaa                                                                                                                              >  1:329384/1‑137 (MQ=255)
                              agcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCaaa                                                                                                                              <  2:329384/138‑2 (MQ=255)
                              agcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGt                                                                                                                   <  5:163301/149‑1 (MQ=255)
                                 aaCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCcgcg                                                                                                                                          >  1:359394/1‑123 (MQ=255)
                                 aaCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGCCcgcg                                                                                                                                          <  2:359394/123‑1 (MQ=255)
                                 aaCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttca                                                                                                                >  6:3233/1‑149 (MQ=255)
                                   ccGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttcatt                                                                                                              <  4:258330/149‑1 (MQ=255)
                                        ccGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATc                                                                                                         <  6:115347/149‑1 (MQ=255)
                                          gTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGt                                                                                                       <  5:3233/149‑1 (MQ=255)
                                               gAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGATTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGCGCCGCAAGGCGGCCCCGGTTCAATCAATCGTAATg                                                                                                   >  4:158550/1‑148 (MQ=255)
                                                     cccGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCg                                                                                            <  4:173771/149‑1 (MQ=255)
                                                      gccagtGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttcat                                                                                                               >  6:360326/6‑129 (MQ=255)
                                                      gccagtGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttcat                                                                                                               <  5:360326/124‑1 (MQ=255)
                                                            gCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCa                                                                                     >  1:77240/1‑149 (MQ=255)
                                                             ccAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCAc                                                                                    <  7:323635/149‑1 (MQ=255)
                                                             ccAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATACCGCTCAAATAATGTAATTGGGGCTGCTTTGCAGCCCTTCCCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCAGGTTCATTCAATCGTATTGAAACCCCCCTGCCa                                                                                     >  1:78584/1‑148 (MQ=255)
                                                                         ggCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATTTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGc                                                                        <  5:177724/149‑1 (MQ=255)
                                                                             cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCTGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATa                                                                    <  7:289985/149‑1 (MQ=255)
                                                                               ggCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAg                                                                  >  2:241299/1‑149 (MQ=255)
                                                                                  ttCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtct                                                               >  7:20768/1‑149 (MQ=255)
                                                                                  ttCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGtct                                                               >  4:333623/1‑149 (MQ=255)
                                                                                     cTGAAGCCAGGAGATTCCGCTCATATAATGTAATTTGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGa                                                                          <  1:172198/135‑1 (MQ=255)
                                                                                     cTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTTGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGGTCATTGAAAAGTATTTATAACCGCCTGCCACACCAGCCTGa                                                                          >  2:172198/1‑135 (MQ=255)
                                                                                             aGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                       >  1:211092/1‑66 (MQ=255)
                                                                                             aGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAg                                                                                                                                       <  2:211092/66‑1 (MQ=255)
                                                                                              ggAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa                                                                   <  4:172394/133‑1 (MQ=255)
                                                                                              ggAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATaa                                                                   >  3:172394/1‑133 (MQ=255)
                                                                                               gagaTTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTAAATAGTATTGATAACCGCCTGCCACTACTGCCTGAGCCATAAGTCTCGAACCAGTAcc                                                   >  3:220964/1‑148 (MQ=255)
                                                                                                   ttCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAg                                               <  3:333623/148‑1 (MQ=255)
                                                                                                      cGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAg                                            <  1:241299/148‑1 (MQ=255)
                                                                                                       gCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAGGGGCCGCAAGGCGGCCCCGGTTGATTCAGTCGTATTGATACCAGCCTGCCACTCCTGCCTGCGGCGTAAGTGTCGACCCAGGACCCAAGCAgg                                           >  6:109710/1‑148 (MQ=255)
                                                                                                              tAATGTAATTGGGGATGCTTTGCAGACCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGa                                   <  3:158550/149‑1 (MQ=255)
                                                                                                                     aTTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCaa                               <  4:305611/146‑1 (MQ=255)
                                                                                                                     aTTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCaa                               >  3:305611/1‑146 (MQ=255)
                                                                                                                     aTTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATc                            <  8:20768/149‑1 (MQ=255)
                                                                                                                         gggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGACGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAAtagatg                        >  3:288742/1‑144 (MQ=255)
                                                                                                                         gggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGTCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCAcc                         <  5:109710/148‑1 (MQ=255)
                                                                                                                         gggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAACCAGGTACCGGGACCaaa                              <  4:288742/143‑1 (MQ=255)
                                                                                                                           gctgctTTGCAGCCCTTCGCGGGCGAGCCCGCGACGGGCCGCAAGGCGGCCCCGGGTCATTCAATCGGATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGCCCCACGACCCAAGCAGGTACCGGGACCAAATCACCggg                      >  7:117680/1‑149 (MQ=255)
                                                                                                                              gctTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                        >  5:324832/1‑112 (MQ=255)
                                                                                                                              gctTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                        <  6:324832/112‑1 (MQ=255)
                                                                                                                                     aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTc                                                                                  <  3:34282/79‑1 (MQ=255)
                                                                                                                                     aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTc                                                                                  >  4:34282/1‑79 (MQ=255)
                                                                                                                                     aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTc            >  3:177445/1‑149 (MQ=255)
                                                                                                                                           cggcgGGCGAGCCCGCTATGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGccc      <  4:177445/147‑1 (MQ=255)
                                                                                                                                               gggCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa  <  2:77240/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGGTGGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAA  >  NC_002947/3048457‑3048748

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