Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,763,502 A→G 6.9% intergenic (+22/+1653) PP_3326 → / ← PP_3328 hypothetical protein/ring‑cleaving dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,763,5020AG 6.9% 162.9 / 2.8 58intergenic (+22/+1653)PP_3326/PP_3328hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (2/2);  total (28/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 = 7.75e-01

GCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTGCA  >  NC_002947/3763359‑3763650
                                                                                                                                               |                                                                                                                                                    
gCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACg                                                                                                                                                 <  4:162830/149‑1 (MQ=255)
  aacgAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGc                                                                                                                                               >  4:11020/4‑149 (MQ=255)
       aGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAgg                                                                                                                                          <  5:270974/149‑1 (MQ=255)
         gACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAGGTa                                                                                                                                        <  5:139866/149‑1 (MQ=255)
          aCTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTAt                                                                                                                                       <  3:148032/149‑1 (MQ=255)
            tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATgg                                                                                                                                     >  6:95253/1‑149 (MQ=255)
            tCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATgg                                                                                                                                     >  4:85357/1‑149 (MQ=255)
              aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATgg                                                                                                                                     <  5:95253/147‑1 (MQ=255)
                 atcatcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGtt                                                                                                                                >  8:83068/1‑149 (MQ=255)
                     tcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                            >  1:52378/1‑149 (MQ=255)
                      cgcgGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                            <  4:259859/148‑1 (MQ=255)
                           cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTc                                                                                                                      <  7:83068/149‑1 (MQ=255)
                               ggCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGaa                                                                                                                  <  3:147963/149‑1 (MQ=255)
                                  tACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                            <  4:29616/136‑1 (MQ=255)
                                  tACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                            >  3:29616/1‑136 (MQ=255)
                                   aCCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                                            <  4:65770/119‑1 (MQ=255)
                                   aCCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGa                                                                                                                                            >  3:65770/1‑119 (MQ=255)
                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAg                                                                                                                                  <  5:326150/127‑1 (MQ=255)
                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAg                                                                                                                                  >  6:326150/1‑127 (MQ=255)
                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGa                                                                                                               <  1:123270/146‑1 (MQ=255)
                                     ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGa                                                                                                               >  2:123270/1‑146 (MQ=255)
                                        gCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                         >  5:79031/1‑149 (MQ=255)
                                         ccAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGcc                                                                                                          <  2:70887/147‑1 (MQ=255)
                                               gCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCa                                                                                                  <  4:163767/149‑1 (MQ=255)
                                               gCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCa                                                                                                  <  2:276443/149‑1 (MQ=255)
                                                           cccGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCa                                                                                      >  5:155616/1‑149 (MQ=255)
                                                            ccGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCAt                                                                                     <  2:3632/149‑1 (MQ=255)
                                                                   ccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGGCGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCAtt                                                                              <  8:296869/149‑1 (MQ=255)
                                                                   ccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACATCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGAGCCGGCCACCAGCATCCGCCCAATGACCCtt                                                                              >  8:89925/1‑149 (MQ=255)
                                                                                   agttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCt                                                               <  7:118764/148‑1 (MQ=255)
                                                                                   agttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTg                                                              >  5:39001/1‑149 (MQ=255)
                                                                                        aCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGc                                                          <  6:106055/148‑1 (MQ=255)
                                                                                        aCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGc                                                          <  6:79031/148‑1 (MQ=255)
                                                                                        aCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGc                                                          <  6:39001/148‑1 (MQ=255)
                                                                                          tGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCa                                                                                      <  8:279808/118‑1 (MQ=255)
                                                                                          tGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCa                                                                                      >  7:279808/1‑118 (MQ=255)
                                                                                            tCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCacac                                                     >  7:147894/1‑149 (MQ=255)
                                                                                                       tAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAac                                          <  3:11020/149‑1 (MQ=255)
                                                                                                       tAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAac                                          <  6:155616/149‑1 (MQ=255)
                                                                                                           ccGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGa                                      >  3:108822/1‑149 (MQ=255)
                                                                                                              caccacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATgg                                   <  3:85357/149‑1 (MQ=255)
                                                                                                                 cacGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCa                                >  1:87437/1‑149 (MQ=255)
                                                                                                                  acGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGTGGTATTGCATTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGc                                                                           <  2:268255/105‑1 (MQ=255)
                                                                                                                  acGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACAATGACCATTGGc                                                                           >  1:268255/1‑105 (MQ=255)
                                                                                                                   cGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATg                              >  4:284877/1‑149 (MQ=255)
                                                                                                                     aGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                         <  7:301631/72‑1 (MQ=255)
                                                                                                                     aGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                         >  8:301631/1‑72 (MQ=255)
                                                                                                                      gTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTc                            >  7:36176/1‑148 (MQ=255)
                                                                                                                            tAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATg                              >  6:15162/1‑140 (MQ=255)
                                                                                                                            tAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATg                              <  5:15162/140‑1 (MQ=255)
                                                                                                                                      aCCTCGGCCGTCATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATgg           >  1:134899/1‑149 (MQ=255)
                                                                                                                                      aCCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGc                  >  1:121079/1‑142 (MQ=255)
                                                                                                                                      aCCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGc                  <  2:121079/142‑1 (MQ=255)
                                                                                                                                       ccTCGGCCGTCATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg          >  6:165445/1‑149 (MQ=255)
                                                                                                                                       ccTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg          <  3:284877/149‑1 (MQ=255)
                                                                                                                                        cTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg          >  5:239769/1‑148 (MQ=255)
                                                                                                                                        cTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAACCCCGGCTGCGGCCCCACGGCAAACGAACACGATGGGCATGTCGAATACCTGCAAGAGGGGc         >  4:138492/1‑149 (MQ=255)
                                                                                                                                              cATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGGCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTg    >  8:185807/1‑148 (MQ=255)
                                                                                                                                               aTGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTGCa  >  3:119041/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTGCA  >  NC_002947/3763359‑3763650

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑