Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,017,454 T→A 7.4% S136S (TCA→TCT)  PP_5461 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,017,4540TA 7.4% 155.0 / 2.9 54S136S (TCA→TCT) PP_5461hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base A (2/2);  total (28/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 = 5.59e-01

GTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGTCAT  >  NC_002947/1017306‑1017599
                                                                                                                                                    |                                                                                                                                                 
gTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCtgat                                                                                                                                                   <  6:88462/149‑1 (MQ=255)
      aTTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACACCTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAacac                                                                                                                                             >  6:303752/1‑149 (MQ=255)
          tGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCa                                                                                                                                           >  7:258862/1‑147 (MQ=255)
            aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCaaa                                                                                                                                         <  8:309186/147‑1 (MQ=255)
            aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAg                                                                                                                                        <  7:159378/148‑1 (MQ=255)
            aTCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCaaa                                                                                                                                         <  8:75569/147‑1 (MQ=255)
                ggTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTc                                                                                                                                   >  6:92033/1‑149 (MQ=255)
                         tGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACAttt                                                                                                                          >  2:48452/1‑149 (MQ=255)
                               gagaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTttcttc                                                                                                                     <  1:150058/148‑1 (MQ=255)
                               gagaTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCg                                                                                                                    >  8:130378/1‑149 (MQ=255)
                                  aTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGaa                                                                                                                                                <  8:124907/118‑1 (MQ=255)
                                  aTTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGaa                                                                                                                                                >  7:124907/1‑118 (MQ=255)
                                        ggggCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATaaa                                                                                                           <  6:38900/149‑1 (MQ=255)
                                             cGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTa                                                                                                      <  8:258862/149‑1 (MQ=255)
                                             cGCTTTGCGGACCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGACTCGTGACATTTTCTTAGGTATATAAATTGTa                                                                                                      >  5:168827/1‑149 (MQ=255)
                                                    cGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGAtttt                                                                                                <  1:48452/148‑1 (MQ=255)
                                                      gCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTgg                                                                                             >  7:112465/1‑149 (MQ=255)
                                                      gCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCaaa                                                                                                                                         >  2:272892/1‑105 (MQ=255)
                                                      gCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAACGCAGGGCCCTGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCATCAGAACACCaaa                                                                                                                                         <  1:272892/105‑1 (MQ=255)
                                                        ccTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTgggg                                                                                           >  5:321127/1‑149 (MQ=255)
                                                            ttGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGa                                                                                         <  5:92033/147‑1 (MQ=255)
                                                               cTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAgtgt                                                                                     <  3:53445/148‑1 (MQ=255)
                                                                               cAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAtttt                                                                     <  6:321127/148‑1 (MQ=255)
                                                                               cAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAAAAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTAttttt                                                                    >  5:14716/1‑149 (MQ=255)
                                                                                         gtgCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                          >  1:223289/1‑149 (MQ=255)
                                                                                          tgCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                          <  2:223289/148‑1 (MQ=255)
                                                                                               ggCAGGGCCCGGATGGCCGGCTTGAGATCAATAAATCCAGGTTGGAGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGCGGAGTGTGCCAATGGCGCATTTTTCGCTGCGGGCCTGTAg                                                    <  8:112465/149‑1 (MQ=255)
                                                                                                 cAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGat                                                  >  2:58739/1‑149 (MQ=255)
                                                                                                  aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATATATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                          >  3:101744/1‑140 (MQ=255)
                                                                                                  aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATATATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGc                                                          <  4:101744/140‑1 (MQ=255)
                                                                                                  aGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGATTGTGACAATGGCGTATTTTTCGCTGCGGTCCTGTAGata                                                 >  3:39831/1‑149 (MQ=255)
                                                                                                    ggCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGGCAATGGCGTATTTTTCGCTGCGGGCCTGTAGatat                                                >  8:77919/1‑148 (MQ=255)
                                                                                                            tGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCa                                       <  7:130378/149‑1 (MQ=255)
                                                                                                               tgtgCTTTAGATTAATAAATTCAGTTTGGTGTTCTGATTAACACGAAAGGCCTAGTGACATTTTGTTCGGTATATTAATTTTATATTTTTGGGGGAGTGTTCCAATGTCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTAATa                                     <  6:168827/145‑1 (MQ=255)
                                                                                                                         aTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTa                          >  2:18712/1‑149 (MQ=255)
                                                                                                                         aTCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTa                          <  5:303752/149‑1 (MQ=255)
                                                                                                                           cAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTAt                         >  5:277132/1‑148 (MQ=255)
                                                                                                                            aaTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATc                        >  3:329991/1‑148 (MQ=255)
                                                                                                                             aTAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCa                       >  8:103320/1‑148 (MQ=255)
                                                                                                                              tAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGTTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCaaa                     <  1:58739/149‑1 (MQ=255)
                                                                                                                              aaaaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCaaa                     <  6:14716/148‑1 (MQ=255)
                                                                                                                               aaaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAat                    >  2:130883/1‑149 (MQ=255)
                                                                                                                                aaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCAt                                      >  7:218472/1‑130 (MQ=255)
                                                                                                                                aaTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCAt                                      <  8:218472/130‑1 (MQ=255)
                                                                                                                                 aTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAAt                                <  7:75916/135‑1 (MQ=255)
                                                                                                                                 aTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAAt                                >  8:75916/1‑135 (MQ=255)
                                                                                                                                 aTCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatc                   >  3:54344/1‑148 (MQ=255)
                                                                                                                                     aGTTTGGTGTTCTGATGAACAACATAGGCCTCGTTAGATTTTCTTAGTTATATAAATTGTAGATTTTTGTGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatctatct              <  4:39831/149‑1 (MQ=255)
                                                                                                                                       tttGGTGTTCATCAGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTg            >  8:194697/1‑149 (MQ=255)
                                                                                                                                           gtgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatct                  >  4:152416/1‑139 (MQ=255)
                                                                                                                                           gtgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAatct                  <  3:152416/139‑1 (MQ=255)
                                                                                                                                            tgtTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCaaaa       <  4:54344/149‑1 (MQ=255)
                                                                                                                                               tCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGt     >  2:81410/1‑148 (MQ=255)
                                                                                                                                                cTGATGAACACCAAAGGCCTCGTGACA‑TTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCgg                                                            >  4:53220/1‑91 (MQ=255)
                                                                                                                                                cTGATGAACACCAAAGGCCTCGTGACA‑TTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCgg                                                            <  3:53220/91‑1 (MQ=255)
                                                                                                                                                   atgaACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGTCAt  <  7:103320/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTGATCATTCTGATCTGGTTCTTTATGCATTGAGATTTTGGGGGCCGCTTTGCGGCCCTTTTGCTACAACTCTTCGCTACAAAAGGCGTGTGCAAGGCAGGGCCCGGATGTCGGGCTTGAGATCAATAAATCCAGTTTGGTGTTCTGATGAACACCAAAGGCCTCGTGACATTTTCTTCGGTATATAAATTGTAGATTTTTGGGGGAGTGTGCCAATGGCGTATTTTTCGCTGCGGGCCTGTAGATATGAGTGCTCATAGTAATGGACTATCAAATCTATCTTGCAAAAGTCAT  >  NC_002947/1017306‑1017599

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