Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,048,605 A→C 10.0% intergenic (+60/‑236) 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,6050AC10.0% 130.3 / 5.4 48intergenic (+60/‑236)PP_2661/PP_2662hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/22);  new base C (2/3);  total (22/26)
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.02e-01

GGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGC  >  NC_002947/3048461‑3048750
                                                                                                                                                |                                                                                                                                                    
ggATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCCCTCAGATAATGTAATTGGGGCTGCTATGCAGCCCTTCGCGGCCGAGc‑‑c‑c                                                                                                                                               >  4:122454/1‑149 (MQ=255)
 gATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGc‑‑c‑c                                                                                                                                               <  6:149975/148‑1 (MQ=255)
 gATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGc‑‑c‑c                                                                                                                                               <  6:154686/148‑1 (MQ=255)
     tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCGGCCCTTCGCGGGCTCGC‑‑C‑CGCGaa                                                                                                                                          <  1:46063/149‑1 (MQ=255)
     tgCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCTCGC‑‑C‑CGCGaa                                                                                                                                          <  2:114962/149‑1 (MQ=255)
          gCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGcc                                                                                                                                     <  5:41664/149‑1 (MQ=255)
                ctggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTATTTGGGGCTGCTTTGCAGCCCTTCGCGGTCTCGC‑‑C‑CGCGAAGGGCTGCaaag                                                                                                                               <  7:201577/149‑3 (MQ=39)
                ctggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAg                                                                                                                                         >  5:55747/1‑139 (MQ=255)
                ctggcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAg                                                                                                                                         <  6:55747/139‑1 (MQ=255)
                   gcctggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAggc                                                                                                                              <  6:91021/147‑1 (MQ=255)
                       ggcAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGcccc                                                                                                                        >  2:310574/1‑149 (MQ=255)
                              aCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGCAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGAGAAGGGCCGCAAGGCCGCCCCGGttca                                                                                                                  >  4:33304/1‑148 (MQ=255)
                              aCCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGttcat                                                                                                                 >  2:101283/1‑149 (MQ=255)
                                     cGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCg                                                                                                          >  3:271844/1‑149 (MQ=255)
                                      gTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATc                                                                                                           <  3:38213/147‑1 (MQ=255)
                                      gTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATc                                                                                                           >  4:38213/1‑147 (MQ=255)
                                          tGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTAtt                                                                                                      <  3:47836/148‑1 (MQ=255)
                                           gAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGGTCCTTCAATCGGATTGa                                                                                                    >  5:208953/1‑149 (MQ=255)
                                            aaaTGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGa                                                                                                    <  7:15136/148‑1 (MQ=255)
                                                             cGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTGGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcct                                                                                  <  7:159962/149‑1 (MQ=255)
                                                                  gCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTcctgcct                                                                              >  3:244014/1‑148 (MQ=255)
                                                                    aGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAg                                                                           <  1:310574/149‑1 (MQ=255)
                                                                     tctccTGGCTTCAGGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGc                                                                          >  2:207158/5‑149 (MQ=255)
                                                                       tccTGGCTTCAGGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGACGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCa                                                                        >  4:120624/3‑149 (MQ=255)
                                                                       cACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGCGACTCGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGa                                                                            <  1:101283/148‑1 (MQ=255)
                                                                            acTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAAGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGt                                                                   >  5:175176/2‑149 (MQ=255)
                                                                              ttCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCCGCCTGa                                                                            >  7:67899/1‑138 (MQ=255)
                                                                              ttCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCCGCCTGa                                                                            <  8:67899/138‑1 (MQ=255)
                                                                                   gAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGTCACTCCTGCCTGAGCCATAAGTCTCGAAc                                                            >  4:214335/1‑149 (MQ=255)
                                                                                    aagcatggAGTTTCCGCTCAGATATTGTAATTGGGGCTGCTTTGCAGCCGTGGGCGGGCGAGC‑‑G‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAAc                                                            <  3:33304/142‑1 (MQ=255)
                                                                                    aaGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAAGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAAc                                                            <  6:175176/148‑1 (MQ=255)
                                                                                    aaGCCAGGAGATTCCGCTCAGATAATGTAATTGGCGATGCATTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGTCCGCTAGGCGGACCCGGTTCAATCAATCATATTGAAACCCGCCTCCCACTCCTACATGAGCCATAAGTCTAGAAc                                                            >  1:102706/1‑148 (MQ=255)
                                                                                     aGCCAGGAGTTTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                          <  3:122454/149‑1 (MQ=255)
                                                                                     aGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCa                                                          <  1:207158/149‑1 (MQ=255)
                                                                                               ttCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGc                                                <  7:172222/149‑1 (MQ=255)
                                                                                                           aaTGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGAc                                    >  4:112034/1‑149 (MQ=255)
                                                                                                               gaaccggggCCGCCTTGCGGCCCTTCGCGGGCTCGC‑‑C‑CGCGAAGGGCTGCAAAGCAGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGCGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCa                                  >  2:211285/6‑147 (MQ=21)
                                                                                                                aaTTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGt                                            <  2:123678/136‑1 (MQ=255)
                                                                                                                aaTTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGt                                            >  1:123678/1‑136 (MQ=255)
                                                                                                                       gctgttTTGCACCCCTGCGCGGGCGTGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCggg                        <  6:208953/149‑1 (MQ=255)
                                                                                                                        ctgctTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGc                       <  4:244014/149‑1 (MQ=255)
                                                                                                                          gctTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACACGCCTGCCACTCCTCCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATAACCGGGCCg                     >  5:295213/1‑149 (MQ=255)
                                                                                                                           ctTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCAGTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCg                     >  7:6151/1‑148 (MQ=255)
                                                                                                                             ttGCGGCCCTTCGCGGGCTCGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTAAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGAAGGTAACGGGAACAAATCACCGGTCCGTAc                  >  6:32526/1‑149 (MQ=38)
                                                                                                                              tGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGTCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTAc                  <  3:214335/148‑1 (MQ=255)
                                                                                                                              tGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTa                   <  8:6151/147‑1 (MQ=255)
                                                                                                                                            ggCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAg   <  1:191211/149‑1 (MQ=255)
                                                                                                                                             gCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATAGTATTGATACACGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATAACCGGGCCGTACCCACCCGCCCCaaatc  >  6:70317/1‑147 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGATGTGCAGGCACGGCTGGCCTGGCAGCAACCGCACCGTTATGAAATGCCCGCTTGCCAGCGTTCGCAGGCACTGGCTTCCTGAAGCCAGGAGATTCCGCTCAGATAATGTAATTGGGGCTGCTTTGCAGCCCTTCGCGGGCGAGC‑‑C‑CGCGAAGGGCCGCAAGGCGGCCCCGGTTCATTCAATCGTATTGATACCCGCCTGCCACTCCTGCCTGAGCCATAAGTCTCGAACCAGGACCCAAGCAGGTACCGGGACCAAATCACCGGGCCGTACCCTCCCGCCCCAAAGC  >  NC_002947/3048461‑3048750

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