Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,995 C→G 7.0% G54A (GGC→GCC)  PP_4045 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,9950CG 7.0% 171.8 / 2.2 55G54A (GGC→GCC) PP_4045hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (26/25);  new base G (2/2);  total (28/27)
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.86e-01

CCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCA  >  NC_002947/4557859‑4558141
                                                                                                                                        |                                                                                                                                                  
ccTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATcg                                                                                                                                        >  4:278081/1‑149 (MQ=255)
 cTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGAt                                                                                                                                          >  1:284381/1‑146 (MQ=255)
 cTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGAt                                                                                                                                          <  2:284381/146‑1 (MQ=255)
                     ccACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTgg                                                                                                                                              <  5:151979/122‑1 (MQ=255)
                     ccACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTgg                                                                                                                                              >  6:151979/1‑122 (MQ=255)
                      cACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGtg                                                                                                                  <  2:93597/149‑2 (MQ=255)
                      cACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTc                                                                                                                  <  8:51278/149‑1 (MQ=255)
                         ttCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGc                                                                                                                >  7:166454/1‑148 (MQ=255)
                              ttAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAAccc                                                                                                          >  5:55840/1‑149 (MQ=255)
                                 aTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAGCGGCTCGGCGCTgatgatgga                                                                                                                                       >  8:15156/1‑114 (MQ=38)
                                 aTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAGCGGCTCGGCGCTgatgatgga                                                                                                                                       <  7:15156/117‑4 (MQ=38)
                                          ccaacTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATcgccgc                                                                                                                                    <  2:41148/111‑1 (MQ=255)
                                          ccaacTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATcgccgc                                                                                                                                    >  1:41148/1‑111 (MQ=255)
                                            aacTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTg                                                                                            >  6:234688/1‑149 (MQ=255)
                                             acTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTg                                                                                                                                               <  3:136611/97‑1 (MQ=255)
                                             acTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTg                                                                                                                                               >  4:136611/1‑97 (MQ=255)
                                              cTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGcc                                                                                          <  8:243381/149‑1 (MQ=255)
                                               tGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCg                                                                                         >  6:288871/1‑149 (MQ=255)
                                                 ttATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                                       <  8:166454/149‑1 (MQ=255)
                                                  tATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGtt                                                                                       <  4:77830/148‑1 (MQ=255)
                                                       ccATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACg                                                                                   >  7:210195/1‑147 (MQ=255)
                                                                 ggACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTg                                                                       <  6:55840/149‑1 (MQ=255)
                                                                  gACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTgg                                                                      <  2:303317/149‑1 (MQ=255)
                                                                     gccgccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCt                                                                                                                     <  7:9311/99‑1 (MQ=255)
                                                                     gccgccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCt                                                                                                                     >  8:9311/1‑99 (MQ=255)
                                                                        gccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGc                                                                                                                                                    <  2:141402/65‑1 (MQ=255)
                                                                        gccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGc                                                                                                                                                    >  1:141402/1‑65 (MQ=255)
                                                                          cATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCgctg                                                              <  5:288871/149‑1 (MQ=255)
                                                                                 ggTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCt                                                                        <  3:138386/132‑1 (MQ=255)
                                                                                 ggTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCt                                                                        >  4:138386/1‑132 (MQ=255)
                                                                                       aTCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggtgg                                                 <  8:210195/149‑1 (MQ=255)
                                                                                          gCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGc                                              <  4:67898/149‑1 (MQ=255)
                                                                                           ccTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGTGcc                                             <  5:234688/149‑1 (MQ=255)
                                                                                            cTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCg                                            >  5:248410/1‑149 (MQ=255)
                                                                                              aCCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGca                                          >  3:10138/1‑149 (MQ=255)
                                                                                                ccTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCg                                                                                                                 <  5:210703/76‑1 (MQ=255)
                                                                                                ccTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCg                                                                                                                 >  6:210703/1‑76 (MQ=255)
                                                                                                  tAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGGGGCCGCCCTGCGGTTGACGAAGAACAGGGTGGAGCCGCTGCTGACGCTGGTGGGGCCGca                                          <  4:203136/145‑1 (MQ=255)
                                                                                                  tAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGca                                          <  4:199194/145‑1 (MQ=255)
                                                                                                  tAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGca                                          >  3:203136/1‑145 (MQ=255)
                                                                                                  tAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGca                                          >  3:199194/1‑145 (MQ=255)
                                                                                                  tAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAAt                                      >  3:15934/1‑149 (MQ=255)
                                                                                                       aGCGGCTGTCCTTCATCCATCATCAGCGGCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCa                                 <  6:248410/149‑1 (MQ=255)
                                                                                                          ggCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATg                              >  7:157247/1‑149 (MQ=255)
                                                                                                            cTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGt                            >  2:1178/1‑149 (MQ=255)
                                                                                                            cTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGt                            <  1:303112/149‑1 (MQ=255)
                                                                                                             tGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGACAACACTGGCGGCCGCCATGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTCGTGGGGCCGCACAATGGGCAATgct                            >  5:96140/1‑146 (MQ=255)
                                                                                                              gTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGt                            >  2:162424/1‑147 (MQ=255)
                                                                                                               tgcTTCATGCATCATCAGCGCCGAGCCGGTGGTGATCGCGGCCCCACAGCCGATCCTGTCGCCAACGGTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGTTGGTGCCGCTGGTGACGCTGGTTGGGCCGCACAATGGGGAATGGTGa                          <  6:96140/146‑1 (MQ=255)
                                                                                                               tCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCACCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAg                         <  3:278081/149‑1 (MQ=255)
                                                                                                                        catcatCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGAcc                  >  4:63564/1‑147 (MQ=255)
                                                                                                                         atcaccAGCGGCTCGGCGCTGATGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggtg                                                  <  6:281799/114‑1 (MQ=255)
                                                                                                                         atcaccAGCGGCTCGGCGCTGATGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCtggtg                                                  >  5:281799/1‑114 (MQ=255)
                                                                                                                             tcaGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAg            >  1:89115/1‑148 (MQ=255)
                                                                                                                             tcaGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAg            >  6:19959/1‑148 (MQ=255)
                                                                                                                                gcgcCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAAc         >  1:316085/1‑148 (MQ=255)
                                                                                                                                      aGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATggca  <  1:1178/149‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                  
CCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAACGATGGCA  >  NC_002947/4557859‑4558141

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