Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,624 G→A 8.7% intergenic (+79/‑217) 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,6240GA 8.7% 127.6 / 2.5 46intergenic (+79/‑217)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base A (2/2);  total (23/23)
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.82e-01

CTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCACCC  >  NC_002947/3048482‑3048761
                                                                                                                                              |                                                                                                                                           
ctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGcc                                                                                                                                       >  3:25144/1‑149 (MQ=255)
      gcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGtt                                                                                                                                 >  4:143927/1‑149 (MQ=255)
      gcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGtt                                                                                                                                 >  7:11853/1‑149 (MQ=255)
      gcaACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGt                                                                                                                                  >  7:7306/1‑148 (MQ=255)
           cGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGttcatt                                                                                                                             >  5:308394/1‑148 (MQ=255)
                cGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCg                                                                                                                       <  3:242548/149‑1 (MQ=255)
                               gCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt                                                                                                        <  8:7306/149‑1 (MQ=255)
                               gCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt                                                                                                        <  8:11853/149‑1 (MQ=255)
                                     cAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACt                                                                                                  >  5:5762/1‑149 (MQ=255)
                                       gCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGc                                                                                                      <  5:291206/143‑1 (MQ=255)
                                       gCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGc                                                                                                      >  6:291206/1‑143 (MQ=255)
                                             gCAGGCACTGGCTTCCTGAAGCCTGGATATTCCGCTCACATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcctgcct                                                                                           <  8:47568/148‑1 (MQ=255)
                                              cAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGa                                                                                         <  6:5762/149‑1 (MQ=255)
                                                    cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGCAGGGCCGCAAGGCGGCCCCGGTTCAGTCAAGCGTATTGATACCCGCCTGCCACGCCAGCCTGAGCCaa                                                                                    >  2:276972/1‑147 (MQ=255)
                                                    cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt                                                                                                        >  4:215710/1‑128 (MQ=255)
                                                    cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCt                                                                                                        <  3:215710/128‑1 (MQ=255)
                                                    cTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATa                                                                                   <  4:25144/149‑1 (MQ=255)
                                                    cTGGCTTCAGGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCAt                                                                                    >  1:158121/1‑148 (MQ=255)
                                                                 gCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAg                                                                    >  2:52800/1‑149 (MQ=255)
                                                                        gaTTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGAccc                                                               >  2:54292/1‑147 (MQ=255)
                                                                          ttCCTCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAg                                                            <  7:261841/148‑1 (MQ=255)
                                                                            ccGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCa                                                          <  1:52800/148‑1 (MQ=255)
                                                                              gCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAgg                                                        >  1:224405/1‑148 (MQ=255)
                                                                                tCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTAc                                                     >  3:48902/1‑149 (MQ=255)
                                                                                  aGATTATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTAc                                                     <  1:276972/147‑1 (MQ=255)
                                                                                         tgaaccggggCCGCCTTGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCAGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCa                                             >  3:247543/7‑148 (MQ=21)
                                                                                         gTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCa                                             >  2:245812/1‑148 (MQ=255)
                                                                                                gggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCAcc                                      <  2:224405/148‑1 (MQ=255)
                                                                                                 ggCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCgg                                    <  3:143927/149‑1 (MQ=255)
                                                                                                      ccttGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGa                                                >  5:270672/3‑132 (MQ=37)
                                                                                                      ccttGCGGCCCTTCGCGGGCTCGCCCGCGAAGGGCTGCAAAGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGa                                                <  6:270672/130‑1 (MQ=37)
                                                                                                           cAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAAt                                          >  8:66450/1‑133 (MQ=255)
                                                                                                           cAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAAt                                          <  7:66450/133‑1 (MQ=255)
                                                                                                            aGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTc                         <  4:48902/149‑1 (MQ=255)
                                                                                                                 ttCGCGGGCGCGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGcc                                 <  1:294621/136‑1 (MQ=255)
                                                                                                                 ttCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGGTCATTCAATCGTATTGATACCCGCCTGCCACTCCAGCCTGAGCCATAAGTCTC‑‑GACCCAGGCCCCACGCAGGTACCGGGGCCAAATCACCGGGcc                                 >  2:294621/1‑136 (MQ=255)
                                                                                                                      gggCGAGCCCGCGAAGGGCCGCAAGGCGGCCCGTGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa               <  6:308394/149‑1 (MQ=255)
                                                                                                                      gggCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa               <  2:158121/149‑1 (MQ=255)
                                                                                                                      gggCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCaaa               <  1:54292/149‑1 (MQ=255)
                                                                                                                          gAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCa            >  6:66272/1‑148 (MQ=255)
                                                                                                                          gAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCAc           >  7:29310/1‑149 (MQ=255)
                                                                                                                            gCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCCAGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCAc           <  1:245812/149‑1 (MQ=255)
                                                                                                                               cgcgAAGGGCTGCAAAGCAGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGcc      <  4:247543/149‑1 (MQ=255)
                                                                                                                                  gAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCAcc   >  3:277878/1‑149 (MQ=255)
                                                                                                                                  gAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCAcc   >  5:199898/1‑149 (MQ=255)
                                                                                                                                   aaGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCAccc  <  6:199898/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCCCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTC‑‑GAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGCACTCGCCACCC  >  NC_002947/3048482‑3048761

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