Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,991,236 G→A 8.1% R131C (CGT→TGT)  PP_3517 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,991,2360GA 8.1% 135.1 / 2.4 49R131C (CGT→TGT) PP_3517hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (23/22);  new base A (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

GCCCGTCATGCTTGCCCGGCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGATGCTGT  >  NC_002947/3991089‑3991376
                                                                                                                                                   |                                                                                                                                            
gCCCGTCATGCTTGCCCGGCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATacga                                                                                                                                             >  6:84785/1‑149 (MQ=255)
           ttGCCCGGCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGAt                                                                                                                                  >  4:235880/1‑149 (MQ=255)
           ttGCCCGGCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGAGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGAt                                                                                                                                  >  6:91166/1‑149 (MQ=255)
                  gCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCa                                                                                                                           <  1:293550/149‑1 (MQ=255)
                  gCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCa                                                                                                                           <  7:207371/149‑1 (MQ=255)
                    tCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACt                                                                                                                         >  8:133171/1‑149 (MQ=255)
                      gCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTg                                                                                                                        <  7:133171/148‑1 (MQ=255)
                        tGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGccc                                                                                                                     >  1:231520/1‑149 (MQ=255)
                        tGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGAGTTCGGCACTGcac                                                                                                                     >  4:66455/1‑147 (MQ=255)
                                  ctGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTg                                                                                                               <  8:93139/145‑1 (MQ=255)
                                       cGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGCCGAGGATTTCGGCACTGCCCGCCATGAGGTCGaac                                                                                                      >  3:83715/1‑148 (MQ=255)
                                              cAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGc                                                                                               >  2:45625/1‑149 (MQ=255)
                                               aGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGcc                                                                                              >  5:251221/1‑149 (MQ=255)
                                                   cACCATTGCCGTTATGGTCCTGAGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCctgg                                                                                          <  5:91166/149‑1 (MQ=255)
                                                    aCCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCctggc                                                                                         >  6:293714/1‑149 (MQ=255)
                                                    aCCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCctggc                                                                                         >  1:307791/1‑149 (MQ=255)
                                                      cATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGccagg                                                                                           >  8:226775/1‑142 (MQ=255)
                                                      cATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGccagg                                                                                           <  7:226775/145‑4 (MQ=255)
                                                            cGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTg                                                                                 >  7:127977/1‑149 (MQ=255)
                                                             gTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTgg                                                                                >  2:12507/1‑149 (MQ=255)
                                                              ttATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGt                                                                               <  2:292927/149‑1 (MQ=255)
                                                                    tCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCCATACGACGGGCGCGGATTTCGGCACAGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTCGTCAGCgc                                                                         >  3:218366/1‑148 (MQ=255)
                                                                       tGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCgg                                                                                                                             <  3:196040/94‑1 (MQ=255)
                                                                       tGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCgg                                                                                                                             >  4:196040/1‑94 (MQ=255)
                                                                              cgTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAg                                                               <  1:12507/149‑1 (MQ=255)
                                                                              cgTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAg                                                               <  4:114779/149‑1 (MQ=255)
                                                                                tAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGt                                                             <  1:45625/149‑1 (MQ=255)
                                                                                    ccATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACAACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATg                                                                     >  2:259218/1‑137 (MQ=255)
                                                                                    ccATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACAACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATg                                                                     <  1:259218/137‑1 (MQ=255)
                                                                                    ccATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACAACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATg                                                                     <  1:255615/137‑1 (MQ=255)
                                                                                    ccATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACAACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATg                                                                     >  2:255615/1‑137 (MQ=255)
                                                                                         tGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGAc                                                     <  3:235880/148‑1 (MQ=255)
                                                                                                   gATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCAcc                                          >  1:114786/1‑149 (MQ=255)
                                                                                                    aTAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCg                                         <  6:251221/149‑1 (MQ=255)
                                                                                                    aTAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCg                                         >  7:185948/1‑149 (MQ=255)
                                                                                                        gCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAg                                     <  2:307791/149‑1 (MQ=255)
                                                                                                         ccAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAg                                     >  8:6620/1‑148 (MQ=255)
                                                                                                              aGCGTCTGCAGAACACACAGCTCGTTGCGCCCGATACGACGGGCGAGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCATTGACCGGCGCCCACCGCCATGATGGc                               <  4:218366/149‑1 (MQ=255)
                                                                                                                   cTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGc                          <  5:293714/149‑1 (MQ=255)
                                                                                                                             cccAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGa                <  7:6620/149‑1 (MQ=255)
                                                                                                                              ccAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGa                >  5:323539/1‑148 (MQ=255)
                                                                                                                              ccAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGAc               <  8:127977/149‑1 (MQ=255)
                                                                                                                              ccAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGACAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGccc                        <  2:148313/140‑1 (MQ=255)
                                                                                                                              ccAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGACAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGccc                        >  1:148313/1‑140 (MQ=255)
                                                                                                                               cAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCTGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACa              >  8:96910/1‑149 (MQ=255)
                                                                                                                                     gTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGa        >  6:282432/1‑149 (MQ=255)
                                                                                                                                         cgcCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCTGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGAtgct    <  7:96910/149‑1 (MQ=255)
                                                                                                                                         cgcCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGAtgct    <  2:114786/149‑1 (MQ=255)
                                                                                                                                           cccGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGAtggtgt  >  6:66443/1‑149 (MQ=255)
                                                                                                                                           cccGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGAtgctgt  <  5:282432/149‑1 (MQ=255)
                                                                                                                                                 acgacgGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCa                                                          <  7:185544/87‑1 (MQ=255)
                                                                                                                                                 acgacgGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCa                                                          >  8:185544/1‑87 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
GCCCGTCATGCTTGCCCGGCTCGCTGAAAATCTTCTGTGCGTACTCCAGCGCACCATTGCCGTTATGGTCCTGGGTCGCGTAGTCCATTTGCGCATCGTGATAGGCCAACAGCGACTGCAGAACACCCAGCTCGTTGCGCCCGATACGACGGGCGCGGATTTCGGCACTGCCCGCCTTGAGGTCGAAGCGCCAGCCCTGGCTGGTTCTGGTCAGCGGTATGGGCAAGGTCCAGTGACCGGCGCCCACCGCCAGGATGGCCTTGCCCTCGCCTGACATGTCGATGCTGT  >  NC_002947/3991089‑3991376

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