Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,557,990 C→G 8.4% A56P (GCG→CCG)  PP_4045 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,557,9900CG 8.4% 175.5 / 4.6 57A56P (GCG→CCG) PP_4045hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/25);  new base G (2/3);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAAC  >  NC_002947/4557844‑4558134
                                                                                                                                                  |                                                                                                                                                
tGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCg                                                                                                                                                >  2:303112/1‑149 (MQ=255)
 gATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGcc                                                                                                                                                 >  7:51278/1‑147 (MQ=255)
  aTCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGa                                                                                                                                               <  3:172344/148‑1 (MQ=255)
               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                                                                                                                                <  7:15156/117‑4 (MQ=38)
                                                aTCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCACCAGCGGCTCGGCGCTgatgatgga                                                                                                                                >  8:15156/1‑114 (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                                                                                                              >  8:9311/1‑99 (MQ=255)
                                                                                    gccgccATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCt                                                                                                              <  7:9311/99‑1 (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                                                                                                          >  6:210703/1‑76 (MQ=255)
                                                                                                               ccTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCg                                                                                                          <  5:210703/76‑1 (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                     <  1:303112/149‑1 (MQ=255)
                                                                                                                           cTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGt                     >  2:1178/1‑149 (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)
                                                                                                                                                  |                                                                                                                                                
TGATCTTGCCCATTCCCTGTGGGGGCTGAACCGCTCCCACTTCGCTTAATCGAGCAACCAACTGTTATCGCCATTTTCCAGGACGCCGCCATGACTGGTAGCATCGCCTACCCTAGCAAGCGGCTGTCCTTCATCCATCATCAGCGCCGAGCCGCTGGTGATCGCCGCCCCACAGCCGAGCCTGTCGCCAACCCTGGCGGCCGCCCTGCCGTTGACGAAGAACAGGCTGGAGCCGCTGCTGACGCTGGTGGGGCCGCACAATGGGCAATGGTGAGGGTGACCAAGCAGAAC  >  NC_002947/4557844‑4558134

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