Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,658,812 C→T 8.0% V128M (GTG→ATG)  PP_1451 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,658,8120CT 8.0% 139.0 / 2.3 50V128M (GTG→ATG) PP_1451hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/24);  new base T (2/2);  total (24/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 = 9.46e-01

CATCCCGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGCCTTCGCGTTCCTTGGCGCCCTT  >  NC_002947/1658665‑1658959
                                                                                                                                                   |                                                                                                                                                   
cATCCCGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTAATCATCGTTGATCGGCGTGCCGTGAATCAAGATCATGCTGCCGGCATTCAc                                                                                                                                                     <  5:145684/148‑1 (MQ=255)
 aTCCCGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGGATTCAccc                                                                                                                                                   >  1:88828/1‑149 (MQ=255)
 aTCCCGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCAccc                                                                                                                                                   >  3:77930/1‑149 (MQ=255)
  tCCCGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCAcccc                                                                                                                                                  >  6:322497/1‑149 (MQ=255)
   cccGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCt                                                                                                                                                 >  2:247530/1‑149 (MQ=255)
     cGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGc                                                                                                                                               >  4:261447/1‑149 (MQ=255)
     cGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGc                                                                                                                                               >  4:184248/1‑149 (MQ=255)
      gAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGcc                                                                                                                                              >  4:74714/1‑149 (MQ=255)
       aaCCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTAATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGc                                                                                                                                               >  8:281710/1‑147 (MQ=255)
        aCCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCg                                                                                                                                            <  1:247530/149‑1 (MQ=255)
         ccAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCgg                                                                                                                                           <  8:46985/149‑1 (MQ=255)
                ccACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTCCCCGGGTcgcgc                                                                                                                                    >  4:54277/1‑149 (MQ=255)
                  acacTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCgcgg                                                                                                                                  >  3:236642/1‑149 (MQ=255)
                        ccTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGggcggc                                                                                                                             <  1:29996/148‑1 (MQ=255)
                             aTGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCt                                                                                                                       >  2:292232/1‑149 (MQ=255)
                                tCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCAcc                                                                                                                                                    >  8:131712/1‑117 (MQ=255)
                                      ttGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCgg                                                                                                              >  6:170419/1‑149 (MQ=255)
                                      ttGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCgg                                                                                                              <  4:77930/149‑1 (MQ=255)
                                          gcATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGGaa                                                                                                          >  6:288347/1‑149 (MQ=255)
                                                cAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCTGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGAtg                                                                                                    <  5:322497/149‑1 (MQ=255)
                                                cAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGAtg                                                                                                    <  3:184248/149‑1 (MQ=255)
                                                 tgcgCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGAtgt                                                                                                   <  1:171292/146‑1 (MQ=255)
                                                      aaCCATCCGTCCATTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCAt                                                                                               <  7:281710/148‑1 (MQ=255)
                                                      aaCCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCAt                                                                                               <  3:261447/148‑1 (MQ=255)
                                                      aaCCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATg                                                                                              >  8:133182/1‑149 (MQ=255)
                                                             cGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGt                                                                                       <  3:74714/149‑1 (MQ=255)
                                                              gtccagtccaGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCAAGGtt                                                                                      >  6:262216/1‑149 (MQ=255)
                                                                   gtccagGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCATCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGc                                                                                 >  3:190301/1‑149 (MQ=255)
                                                                   gtcaagGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCATCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGc                                                                                 >  3:190302/1‑149 (MQ=255)
                                                                     ccagGTTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTTCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGTGTAATTGATGTTCATGGCCAGGTTGAAGCg                                                                                <  5:262216/148‑1 (MQ=255)
                                                                      cagGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGACGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGt                                                                              >  1:76599/1‑149 (MQ=255)
                                                                             tGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGctct                                                                        <  6:122792/148‑1 (MQ=255)
                                                                                    aCCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATTATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGc                                                                           <  2:246612/138‑1 (MQ=255)
                                                                                    aCCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATAATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGAGGTGCATGGCCAGGTTGAAGCGGTCGc                                                                           >  1:246612/1‑138 (MQ=255)
                                                                                     ccATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCGCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTgcgc                                                                <  3:54277/148‑1 (MQ=255)
                                                                                        ttCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGa                                                                                                                     <  1:97201/92‑1 (MQ=255)
                                                                                        ttCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGa                                                                                                                     >  2:97201/1‑92 (MQ=255)
                                                                                          ccGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTc                                                           <  1:236262/148‑1 (MQ=255)
                                                                                            gggTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTccag                                                        <  5:30151/149‑1 (MQ=255)
                                                                                              gTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTccagcc                                                      <  2:76599/149‑1 (MQ=255)
                                                                                                                  tGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGtt                                  >  4:167341/1‑149 (MQ=255)
                                                                                                                          atcatgatcatgCTGCCGGGATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTc                          <  2:88828/149‑1 (MQ=255)
                                                                                                                           tcatgatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCa                                                              >  5:215401/1‑112 (MQ=255)
                                                                                                                           tcatgatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCa                                                              <  6:215401/112‑1 (MQ=255)
                                                                                                                           tcatgatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCgggg                                     >  2:332312/1‑137 (MQ=255)
                                                                                                                           tcatgatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCgggg                                     <  1:332312/137‑1 (MQ=255)
                                                                                                                            catgatcatgCTGCCGGCATTCATCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGc                        <  4:190302/149‑1 (MQ=255)
                                                                                                                             atgatcatgCTGCCGGCATTCATCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGc                        <  4:190301/148‑1 (MQ=255)
                                                                                                                               gatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCgggg                                     >  6:286327/1‑133 (MQ=255)
                                                                                                                               gatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCgggg                                     <  5:286327/133‑1 (MQ=255)
                                                                                                                               gatcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGCCtt                     >  3:74547/1‑149 (MQ=255)
                                                                                                                                 tcatgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGCCTTcg                   <  7:133182/149‑1 (MQ=255)
                                                                                                                                  catgCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTTAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGTGTTTTCTTGTCGCCTTcgc                  >  6:171991/1‑149 (MQ=255)
                                                                                                                                           ggCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGCCTTCGCGTTCCTTgg         <  3:167341/149‑1 (MQ=255)
                                                                                                                                                  aCCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGCCTTCGCGTTCCTTGGCGCCCtt  >  8:286899/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CATCCCGAACCAGGTCCCACACTTCCTGCATGTCGCGGTTGCGCATGGCAATGCAACCATCCGTCCAGTCCAGGGTATGGAAGTACCATTCCGGGTATTCATCGTTGATCGGCGTGCCGTGAATCATGATCATGCTGCCGGCATTCACCCCTGCCCGGGTCGCGCGGGCGGCATCGCTGATGTTCGGGTAATTGATGTGCATGGCCAGGTTGAAGCGGTCGCTCTGCTTGCGCCAGTCCAGCCAGTACAAGCCCTCGGGGGTTTTCTTGTCGCCTTCGCGTTCCTTGGCGCCCTT  >  NC_002947/1658665‑1658959

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