Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,208,773 C→T 7.0% A242A (GCC→GCT)  PP_3692 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,208,7730CT 7.0% 161.9 / 2.1 57A242A (GCC→GCT) PP_3692hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/2);  total (27/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.26e-01

ACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATCGATAAAGAGTGCTCG  >  NC_002947/4208628‑4208920
                                                                                                                                                   |                                                                                                                                                   
aCATCCCAAGCGGTCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGt                                                                                                                                                   >  7:38713/1‑148 (MQ=255)
   tCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAac                                                                                                                                               <  8:330393/149‑1 (MQ=255)
    cccAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAaca                                                                                                                                              >  1:176418/1‑149 (MQ=255)
        aGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACaga                                                                                                                                          <  4:233087/149‑1 (MQ=255)
          cGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTTGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGt                                                                                                                                        <  5:267306/149‑1 (MQ=255)
           ggCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGtt                                                                                                                                       <  3:152552/149‑1 (MQ=255)
                  acaAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg                                                                                                                                >  4:168176/1‑149 (MQ=255)
                   caAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAAGTGGGTGATATGAGCCAAGTCGCATGAAATTTTTGTCGGCACTCATAGCCACTGACCATGACCTCAACACAGAGTTTACAAg                                                                                                                                 >  1:282310/1‑149 (MQ=255)
                    aaaTCGTTGCGCCAGGCGAAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATTGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTAATCGCCACTGCCGATGCCGTCAaca                                                                                                                                              <  7:8868/133‑1 (MQ=255)
                    aaaGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg                                                                                                                                >  1:156798/1‑147 (MQ=255)
                    aaaGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg                                                                                                                                >  2:107325/1‑147 (MQ=255)
                    aaaGAGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGATCAGCAATACATGTGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTAGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGACGTCAaca                                                                                                                                              >  8:8868/1‑133 (MQ=255)
                      aGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGgaga                                                                                                                            <  1:107325/149‑1 (MQ=255)
                       gCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGaa                                                                                                                           <  8:38713/149‑1 (MQ=255)
                               ccAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGc                                                                                                                   <  2:179305/149‑1 (MQ=255)
                                     gcAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAg                                                                                                             >  8:50461/1‑149 (MQ=255)
                                     gcAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAg                                                                                                             >  2:301053/1‑149 (MQ=255)
                                        aGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAg                                                                                                          <  3:168176/149‑1 (MQ=255)
                                        aGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAg                                                                                                          <  7:50461/149‑1 (MQ=255)
                                                   gTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAAc                                                                                               <  2:282310/149‑1 (MQ=255)
                                                    tAAAGCCGCTCAGCAATACATGGGTCTAACTGAA‑‑GCTTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAgg                                                                                                                                <  8:326713/115‑1 (MQ=255)
                                                    tAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGaaaatg                                                                                                                                >  7:326713/1‑113 (MQ=255)
                                                        gCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGtt                                                                                                                                       <  5:358408/104‑1 (MQ=255)
                                                        gCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGtt                                                                                                                                       >  6:358408/1‑104 (MQ=255)
                                                           gCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAg                                                                                       <  2:156798/149‑1 (MQ=255)
                                                             tCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGct                                                                                     >  5:96899/1‑149 (MQ=255)
                                                                 cAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCATGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAAc                                                                                               >  2:308210/1‑135 (MQ=255)
                                                                 cAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCATGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAAc                                                                                               <  1:308210/135‑1 (MQ=255)
                                                                 cAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAggg                                                                                                                               >  4:224033/1‑103 (MQ=255)
                                                                 cAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAggg                                                                                                                               <  3:224033/103‑1 (MQ=255)
                                                                  aaTACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAg                                                                                                                          <  7:87938/107‑1 (MQ=255)
                                                                  aaTACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAg                                                                                                                          >  8:87938/1‑107 (MQ=255)
                                                                  aaTACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCtt                                                                                <  2:176418/149‑1 (MQ=255)
                                                                       aTGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTtcctc                                                                           <  6:96899/149‑1 (MQ=255)
                                                                              gAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCAccgcc                                                                    >  3:193781/1‑149 (MQ=255)
                                                                                            gCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAAt                                                          <  3:211336/147‑1 (MQ=255)
                                                                                            gCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAAt                                                          >  4:211336/1‑147 (MQ=255)
                                                                                                tGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCTGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTa                                                     >  5:89568/1‑148 (MQ=255)
                                                                                                tGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCTGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTa                                                     >  5:84468/1‑148 (MQ=255)
                                                                                                   gCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAAc                                                                                               <  8:113967/103‑1 (MQ=255)
                                                                                                    ccAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAAt                                                          <  7:285190/139‑1 (MQ=255)
                                                                                                    ccAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAAt                                                          >  8:285190/1‑139 (MQ=255)
                                                                                                           gCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACaa                                                  <  5:287300/140‑1 (MQ=255)
                                                                                                           gCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACaa                                                  >  6:287300/1‑140 (MQ=255)
                                                                                                           gCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTaa                                          >  7:44018/1‑148 (MQ=255)
                                                                                                              tGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAAt                                                          >  7:129300/1‑129 (MQ=255)
                                                                                                              tGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAAt                                                          <  8:129300/129‑1 (MQ=255)
                                                                                                                 aaTTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCAc                                    >  6:206267/1‑148 (MQ=255)
                                                                                                                       tgtCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACg                             <  8:44018/149‑1 (MQ=255)
                                                                                                                                tCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAg                     >  4:306308/1‑148 (MQ=255)
                                                                                                                                  aTCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTc                        >  7:303919/1‑143 (MQ=255)
                                                                                                                                  aTCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTc                        <  8:303919/143‑1 (MQ=255)
                                                                                                                                   tCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATc                 <  5:206267/149‑1 (MQ=255)
                                                                                                                                        aCTGCCGATGCTGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATCGATaa            <  6:84468/149‑1 (MQ=255)
                                                                                                                                        aCTGCCGATGCTGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATCGATaa            <  6:89568/149‑1 (MQ=255)
                                                                                                                                            ccGATGCCGGCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATCGATAAaga         <  6:193369/148‑1 (MQ=255)
                                                                                                                                                 gCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTAc                                                    >  8:23578/1‑100 (MQ=255)
                                                                                                                                                   cGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATCGATAAAGAGTGCTCg  >  8:20503/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACATCCCAAGCGGCCCTCACAAAGCGTTGCGCCAGGCGCAAGTGTTCGATAGTAAAGCCGCTCAGCAATACATGGGGCGAACTGAA‑‑GCGTGCTATGAGCCAAGTCGCATGAAATTTGTGTCGGCACTCATCGCCACTGCCGATGCCGTCAACACAGAGTTGACAAGGGAGAAGGAAAGGCTCAAAGGAGCCCTTCCTAACCTACCAAGCTCTCTTTCCTCCACCGCCGAAGCAAAATGGTTACAAGCACTTAAAAGCACGACTACGGCGTCGAGCATCGATAAAGAGTGCTCG  >  NC_002947/4208628‑4208920

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