Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,199,890 C→T 7.1% L42L (CTC→CTT)  PP_2807 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,199,8900CT 7.1% 162.5 / 2.8 56L42L (CTC→CTT) PP_2807hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/28);  new base T (2/2);  total (26/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 = 9.34e-01

AAACCAGACAGGTAGAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGTTCCTTGGCG  >  NC_002947/3199743‑3200035
                                                                                                                                                   |                                                                                                                                                   
aaaCCAGACAGGTAGAAAGCAGATGTATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCa                                                                                                                                                    <  8:273517/149‑1 (MQ=255)
aaaCCAGACAGGTAGAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCa                                                                                                                                                    <  7:210278/149‑1 (MQ=255)
    cAGACAGGTAGAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCg                                                                                                                                                 >  6:323516/1‑148 (MQ=255)
      gACAGGTAGAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCgaag                                                                                                                                              <  4:234719/149‑1 (MQ=255)
           gTAGAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACt                                                                                                                                          <  3:280131/148‑1 (MQ=255)
            tAGAAAGCAGATGGATACACGAGCTCTCTCCACAGTCAGCAGCGATGACTCGATCGAGCAGGTCGTGCACATCATTCAGCGTGACGGCGGCGTCAACATTGACGACATTGTTTCTCCGGCAACGCTGAACTCCCACACCGAACAACTg                                                                                                                                         >  4:280131/1‑148 (MQ=255)
              gAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGa                                                                                                                                       <  5:265163/148‑1 (MQ=255)
               aaaGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGa                                                                                                                                       >  3:42376/1‑147 (MQ=255)
                   cAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTa                                                                                                                                 <  3:27513/149‑1 (MQ=255)
                    acatcgATACCCTAGCCCTCTCCCCACTCAGCATCGATGCCTCGATCGAGCAGGTGTTGGACATCATTGGGCTGGACGCCGGAGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTAc                                                                                                                                <  8:68393/144‑1 (MQ=255)
                         gATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGaa                                                                                                                           <  6:126225/149‑1 (MQ=255)
                              acGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTg                                                                                                                                         >  7:61479/1‑130 (MQ=255)
                              acGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTg                                                                                                                                         <  8:61479/130‑1 (MQ=255)
                              acGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCa                                                                                                                      >  4:284259/1‑149 (MQ=255)
                              acGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCa                                                                                                                      >  2:139786/1‑149 (MQ=255)
                              acGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCa                                                                                                                      >  7:322270/1‑149 (MQ=255)
                                gAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAg                                                                                                                          <  2:349878/143‑1 (MQ=255)
                                gAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAg                                                                                                                          >  1:349878/1‑143 (MQ=255)
                                     ctctCCACAGTCAGCAGCGATGCCTCGATCCAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTtg                                                                                                               >  1:29691/1‑149 (MQ=255)
                                          cacaGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTgtggtg                                                                                                          <  4:48280/149‑1 (MQ=255)
                                          cacaGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTgtggt                                                                                                           <  3:284259/148‑1 (MQ=255)
                                             aGTCAGCAGCGATGGCGCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGGCTGTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGATGAACTGGACAGCTGCCTGAAGGTCACGGCTTGTGgtgt                                                                                                         >  4:107459/1‑147 (MQ=255)
                                             aGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGgtgt                                                                                                         <  3:107459/147‑1 (MQ=255)
                                                              cgatcgaGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGc                                                                                      >  6:307244/1‑149 (MQ=255)
                                                               gatcgaGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGc                                                                                      <  8:322270/148‑1 (MQ=255)
                                                                   gaGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCAccc                                                                                  <  3:125944/148‑1 (MQ=255)
                                                                     gCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTTCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTTACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCAAGa                                                                               >  8:101590/1‑149 (MQ=255)
                                                                     gCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGATGAAGTCCCGTACCGAAGAACTGGACAGCTAACTGAAGGTCACGCCTTGTGGGGTCGATCCATATTTCGCCTGCACACAGa                                                                               >  8:96319/1‑149 (MQ=255)
                                                                        ggtggtGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTTTCGATCCATATTTCGCCGGCACCCAGAccc                                                                            <  7:180617/149‑1 (MQ=255)
                                                                                  atcctTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGACGGCAACGCTGAAGTCCCTTACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGc                                                                  <  7:96319/149‑1 (MQ=255)
                                                                                  atcatTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTTACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGc                                                                  <  7:101590/149‑1 (MQ=255)
                                                                                        cGGCGTGACGGCGGCGTCATCATTGCGTACTTTGTTTCGCCGGCAACGGTGAAGTCCCTCACCGAATAACTGTACAGCTACCTTAAGGTCACGCATTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGaa                                                           <  1:22360/148‑1 (MQ=255)
                                                                                          gCGTGACTGCGGCTTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGATGGTCACGGCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGtgtg                                                                    <  2:130400/139‑1 (MQ=255)
                                                                                          gCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCAGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGtgtg                                                                    >  1:130400/1‑139 (MQ=255)
                                                                                             tGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTg                                                     <  5:323516/149‑1 (MQ=255)
                                                                                             tGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTg                                                     <  5:307244/149‑1 (MQ=255)
                                                                                             tGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTg                                                     <  2:29691/149‑1 (MQ=255)
                                                                                               aCGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTggg                                                   >  3:259495/1‑149 (MQ=255)
                                                                                                        tcatcaTTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGa                                           >  4:268445/1‑148 (MQ=255)
                                                                                                                 ccGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTg                                 >  2:106472/1‑149 (MQ=255)
                                                                                                                               ccGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTCCCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATGGGGCGCCGCGCCACCGCCGTGGAGGTGGCagtagc                   >  3:109237/1‑144 (MQ=255)
                                                                                                                               ccGGCAACGCGGAAGTCCCGCCCCGACGAACTGGACAGCTACCTGAAGGGCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCAt                       <  4:109237/145‑1 (MQ=255)
                                                                                                                                cGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCAc                  >  7:237544/1‑149 (MQ=255)
                                                                                                                                 ggCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCAcc                 <  1:139786/149‑1 (MQ=255)
                                                                                                                                  gCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCActg                >  1:315319/1‑147 (MQ=255)
                                                                                                                                    aaCGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGc              >  3:21225/1‑149 (MQ=255)
                                                                                                                                    aaCGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGc              >  5:127961/1‑149 (MQ=255)
                                                                                                                                      cGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTgg                                                    >  5:48234/1‑109 (MQ=255)
                                                                                                                                      cGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTgg                                                    <  6:48234/109‑1 (MQ=255)
                                                                                                                                        cTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGtt          <  4:259495/149‑1 (MQ=255)
                                                                                                                                              tCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGTTCCTTgg    >  6:277279/1‑149 (MQ=255)
                                                                                                                                                ccTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCT‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGTTCCTTGgcg  >  8:83661/1‑149 (MQ=255)
                                                                                                                                                ccTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGTTCCTTGgc   <  4:42376/148‑1 (MQ=255)
                                                                                                                                                 cTCACCGAACAAAAGGCCAGCTACCCGTAAGGGACGCCATGTGGGGTCGATACATATTTTGCCCGCACCCAGACACGCCCTGTGGCCCTCGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGTTCCTTGgc   <  1:106472/149‑1 (MQ=255)
                                                                                                                                                  tCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAg                              >  8:57829/1‑119 (MQ=255)
                                                                                                                                                  tCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAg                              <  7:57829/119‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAACCAGACAGGTAGAAAGCAGATGGATACACGAGCGCTCTCCACAGTCAGCAGCGATGCCTCGATCGAGCAGGTGGTGGACATCATTCGGCGTGACGGCGGCGTCATCATTGCCGACTTTGTTTCGCCGGCAACGCTGAAGTCCCTCACCGAAGAACTGGACAGCTACCTGAAGGTCACGCCTTGTGGTGTCGATCCATATTTCGCCGGCACCCAGACCCGCCGTGTGGCC‑‑CGAATCATTGGGCGCAGCGACACCGCCGTGGAGGTGGCATTGCACCCGCTGTTCCTTGGCG  >  NC_002947/3199743‑3200035

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