Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,785,977 C→A 7.7% Y327* (TAC→TAA) ‡ PP_5068 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,785,9770CA 7.7% 178.7 / 2.1 57Y327* (TAC→TAA) ‡PP_5068hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base A (3/2);  total (28/29)
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

GCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGA  >  NC_002947/5785829‑5786120
                                                                                                                                                    |                                                                                                                                                  
gCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTAc                                                                                                                                                    >  2:4891/1‑149 (MQ=255)
     cACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGc                                                                                                                                               <  6:120726/149‑1 (MQ=255)
          cgacgaCCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaac                                                                                                                                          <  4:195591/149‑3 (MQ=255)
                            cAGCTACGCCAACCCTGTGCCCATCGCCCACTGTCACCTCAGCGCGCTGGCCCAGGAACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgcgc                                                                                                                         <  5:293472/148‑1 (MQ=255)
                             aGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCg                                                                                                                       >  3:134368/1‑149 (MQ=255)
                                tACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                                    >  8:108099/1‑149 (MQ=255)
                                tACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcg                                                                                                                    >  1:159101/1‑149 (MQ=255)
                                  cGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGgcgc                                                                                                                   >  4:9720/1‑148 (MQ=255)
                                   gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTTGGCCTGTTCCTGCCGCTGTAACAGGCCCATTACCTGGACCTGATGCGACGgcgca                                                                                                                  >  6:244348/1‑147 (MQ=255)
                                   gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGa                                                                                                                                   <  1:126027/131‑1 (MQ=255)
                                   gCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGa                                                                                                                                   >  2:126027/1‑131 (MQ=255)
                                     cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaa                                                                                                                                           <  6:263299/121‑2 (MQ=38)
                                     cAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGGTACAGCGGCAGGAACAGGCCCaa                                                                                                                                           >  5:263299/1‑120 (MQ=38)
                                         cAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCt                                                                                                           <  4:134368/149‑1 (MQ=255)
                                          aGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                                      >  2:165251/1‑121 (MQ=255)
                                          aGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCt                                                                                                                                      <  1:165251/121‑1 (MQ=255)
                                           ggTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGa                                                                                                         <  1:4891/149‑1 (MQ=255)
                                                 cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                          >  2:155112/1‑142 (MQ=255)
                                                 cATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                          <  1:155112/142‑1 (MQ=255)
                                                   tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTa                                                                                                                                         >  1:148964/1‑109 (MQ=255)
                                                   tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTa                                                                                                                                         <  2:148964/109‑1 (MQ=255)
                                                   tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                                  >  7:115976/1‑148 (MQ=255)
                                                   tCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCAccgc                                                                                                  <  8:115976/148‑1 (MQ=255)
                                                    cGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCAGCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                             <  4:67698/136‑1 (MQ=255)
                                                    cGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAggg                                                                                                             >  3:67698/1‑136 (MQ=255)
                                                      cccACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGaa                                                                                              >  1:344833/1‑149 (MQ=255)
                                                           tGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgc                                                                                                                           <  3:98220/115‑1 (MQ=255)
                                                           tGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATgc                                                                                                                           >  4:98220/1‑115 (MQ=255)
                                                                ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTAcc                                                                                                                                       >  1:139433/1‑98 (MQ=255)
                                                                ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTAcc                                                                                                                                       <  2:139433/98‑1 (MQ=255)
                                                                ccTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACa                                                                                 >  3:344653/1‑149 (MQ=255)
                                                                 cTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑tgatga                                                                                   >  2:217648/1‑146 (MQ=255)
                                                                 cTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑tgatga                                                                                   <  1:217648/146‑1 (MQ=255)
                                                                   cAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATgg                                                                              >  5:280011/1‑149 (MQ=255)
                                                                      cgcgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGc                                                                           <  7:108099/149‑1 (MQ=255)
                                                                        cgcTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTg                                                                         <  6:280011/149‑1 (MQ=255)
                                                                             gCCCCGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGa                                                                    <  4:346855/149‑1 (MQ=255)
                                                                               ccAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAgc                                                                  <  3:9720/149‑1 (MQ=255)
                                                                                 aGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                                  <  2:222101/86‑1 (MQ=255)
                                                                                 aGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGAc                                                                                                                                  >  1:222101/1‑86 (MQ=255)
                                                                                    cACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCctg                                                             <  4:344653/149‑1 (MQ=255)
                                                                                           caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                          >  7:201756/1‑100 (MQ=255)
                                                                                           caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTg                                                                                                          <  8:201756/100‑1 (MQ=255)
                                                                                           caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGt                                                      <  2:159101/149‑1 (MQ=255)
                                                                                           caccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGCGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCAGACAACCGCCGAA‑‑GA‑TGATGACATGGTGCTTGTCGTGCGGCTGCTGCAGt                                                      >  6:203901/1‑149 (MQ=255)
                                                                                             ccGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCTAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgc                                                    <  5:244348/149‑1 (MQ=255)
                                                                                                tGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGtgcatg                                                 >  3:92324/1‑149 (MQ=255)
                                                                                                     gAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAAcg                                            >  6:89036/1‑149 (MQ=255)
                                                                                                           gAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAAATAAGTGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGAgg                                         >  1:294156/1‑149 (MQ=255)
                                                                                                                 cccatcgaagCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTg                                  <  1:305358/141‑1 (MQ=255)
                                                                                                                 cTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGa                                >  7:219081/1‑149 (MQ=255)
                                                                                                                        aggtaatGGGCCTGGTACAGCGGCAGGAACAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGc                         >  7:145827/7‑149 (MQ=38)
                                                                                                                           cGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTg                      <  5:89036/149‑1 (MQ=255)
                                                                                                                                        ccTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAAc         >  5:107847/1‑149 (MQ=255)
                                                                                                                                             cGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGt    <  8:219081/149‑1 (MQ=255)
                                                                                                                                              gCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACt        <  2:131705/144‑1 (MQ=255)
                                                                                                                                              gCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACt        >  1:131705/1‑144 (MQ=255)
                                                                                                                                               cTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGa  <  6:107847/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCAACCACTCCGACGACCGTGGCCGCTACAGCTACGCCAACCAGGTGCCCATCGCCCACTGGAACCTCAGCGCGCTGGCCCAGGCACTGACCACCGTGATCGAAGTTGAGCCACTGAAGGAAGCGTTGGGCCTGTTCCTGCCGCTGTACCAGGCCCATTACCTGGACCTGATGCGCCGGCGCCTAGGGCTGACCACCGCCGAA‑‑GA‑TGATGACATGGTGCTGGTCGAGCGCCTGCTGCAGTGCATGCAACGCGGTGGCGTGGACTACAGCCTGTTCTTCCGCAAACTGGGTGA  >  NC_002947/5785829‑5786120

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