Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,261,795 A→G 8.0% Y111H (TAC→CAC)  PP_2859 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,261,7950AG 8.0% 139.4 / 3.2 50Y111H (TAC→CAC) PP_2859hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (24/22);  new base G (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

ATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGTCCTCCCATGCCAGTTGT  >  NC_002947/3261653‑3261940
                                                                                                                                               |                                                                                                                                                  
aTTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCg                                                                                                                                              >  3:272522/1‑149 (MQ=255)
               ttGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTGACGCCGTGCCAGCAGTGCctc                                                                                                                               >  6:6683/1‑149 (MQ=255)
               ttGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTGACGCCGTGCCAGCAGTGCctc                                                                                                                               >  6:1970/1‑149 (MQ=255)
                  tgCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTc                                                                                                                            >  1:317127/1‑149 (MQ=255)
                    ccgccTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTcca                                                                                                                          <  4:272522/149‑1 (MQ=255)
                            gCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCACCAGGTGGAGGTCGTAACGCCGTCCCAGCAGTGCCTCTTCCAGCCAggat                                                                                                                  >  4:287864/1‑147 (MQ=255)
                             ccACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTg                                                                                                                >  3:143476/1‑149 (MQ=255)
                              cACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTgg                                                                                                               >  3:137464/1‑149 (MQ=255)
                                cGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTGACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGc                                                                                                             <  5:1970/149‑1 (MQ=255)
                                cGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTGACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGc                                                                                                             <  5:6683/149‑1 (MQ=255)
                                   cTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGt                                                                                                          >  4:271512/1‑149 (MQ=255)
                                        tGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGcc                                                                                                      >  4:6819/1‑148 (MQ=255)
                                          tcttcgTAAAAGGGCTGACGATGACTGTGGTCGGGTTGGCTACGCTGCCCATCGGCAATCCTGTCCACCCCTTGAGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGaa                                                                                                   <  7:143161/146‑1 (MQ=255)
                                                 aaaaGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATGCCAGACCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGc                                                                                             >  1:251435/1‑149 (MQ=255)
                                                    aGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCTAACAGTGCgggg                                                                                         <  4:137464/149‑1 (MQ=255)
                                                      cgcTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCGATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAAGAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCt                                                                                       <  3:287864/149‑1 (MQ=255)
                                                      cgcTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCcagc                                                                                                                                       >  2:54624/1‑102 (MQ=255)
                                                      cgcTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCcagc                                                                                                                                       <  1:54624/102‑1 (MQ=255)
                                                           aCGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCaga                                                                                  <  3:271512/149‑1 (MQ=255)
                                                             gATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCagag                                                                                 >  7:141017/1‑148 (MQ=255)
                                                              aTCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGACCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGc                                                                                <  1:211487/149‑1 (MQ=255)
                                                                cACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATg                                                                             >  2:219549/1‑149 (MQ=255)
                                                                 aCTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGt                                                                            >  2:139961/1‑149 (MQ=255)
                                                                  cTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGtt                                                                           <  4:143476/149‑1 (MQ=255)
                                                                     aGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTg                                                                          <  2:251435/147‑1 (MQ=255)
                                                                              tGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCTGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAggtg                                                               <  2:317127/149‑1 (MQ=255)
                                                                                 cAACGCTGCCCAACGGCAAACCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCCTGCCAGCAGAGCCTCTTCCAGCCAGGTAGGGCAGTCGCCGAACAGTGACAGGATCCAGAGCATGTTGCTTAGCAggtgggt                                                            >  8:165709/1‑149 (MQ=255)
                                                                                  aaCGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTa                                                                      >  8:181979/1‑138 (MQ=255)
                                                                                  aaCGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTa                                                                      <  7:181979/138‑1 (MQ=255)
                                                                                     gCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTcc                                                        <  1:139961/149‑1 (MQ=255)
                                                                                           cATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAAtt                                                  <  3:203859/149‑1 (MQ=255)
                                                                                            aTCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAAt                                                   >  4:25183/1‑147 (MQ=255)
                                                                                             tCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGt                                                                            <  4:102262/121‑1 (MQ=255)
                                                                                             tCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGt                                                                            >  3:102262/1‑121 (MQ=255)
                                                                                                 cAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTagcagc                                            >  4:129095/1‑149 (MQ=255)
                                                                                                    tCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCggg                                         <  6:207676/149‑1 (MQ=255)
                                                                                                               cccTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCt                               <  8:272429/149‑1 (MQ=255)
                                                                                                                  ttGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTccc                            >  5:205976/1‑149 (MQ=255)
                                                                                                                       ggCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCa                       <  3:6819/149‑1 (MQ=255)
                                                                                                                        gCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCAc                      >  7:141347/1‑149 (MQ=255)
                                                                                                                               aaCAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGtcctc               >  8:275806/1‑149 (MQ=255)
                                                                                                                                  aGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGTCCTccc             >  6:192645/1‑148 (MQ=255)
                                                                                                                                   ggTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGtcctcc              >  3:51684/1‑146 (MQ=255)
                                                                                                                                   ggTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGtcctcc              <  4:51684/146‑1 (MQ=255)
                                                                                                                                    gTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCa                       <  2:41490/136‑1 (MQ=255)
                                                                                                                                    gTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCa                       >  1:41490/1‑136 (MQ=255)
                                                                                                                                       gAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGTCCTCCCATGCCa       <  8:141017/149‑1 (MQ=255)
                                                                                                                                          gTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGtcctc               >  4:55026/1‑138 (MQ=255)
                                                                                                                                          gTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGtcctc               <  3:55026/138‑1 (MQ=255)
                                                                                                                                             gTAACGCCGTGCCAGGAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGTCCTCCCATGCCAGTTGt  <  6:205976/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ATTTGGTGCGGTTGATTGTGCCGCCTCAGCCACGCCTGGCTGTCGTCGAAAAAGCGCTGACGATCACTGAGGTCGGGTTGGCAACGCTGCCCATCGGCAATCCAGTCCACCCCTTGCGGGCTGAGCAACAGGTGGAGGTCGTAACGCCGTGCCAGCAGTGCCTCTTCCAGCCAGGTTGGGCAGTCGCCGAACAGTGCGTGGCTCCAGAGCATGTTGCTTAGCAGGTGGGTGTCCAGAATTAGCAGCGGGGGCGCCTGCTCCCGCGCACGGTCCTCCCATGCCAGTTGT  >  NC_002947/3261653‑3261940

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: