Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,311,359 T→A 11.1% F12Y (TTT→TAT) ‡ PP_3784 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,311,3590TA11.1% 150.9 / 10.8 54F12Y (TTT→TAT) ‡PP_3784hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/23);  new base A (3/3);  total (28/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 = 3.53e-01

ACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGC  >  NC_002947/4311211‑4311507
                                                                                                                                                    |                                                                                                                                                    
acacAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAttt                                                                                                                                                     >  1:26486/1‑150 (MQ=255)
    aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGt                                                                                                                                                 >  2:550231/1‑150 (MQ=255)
    aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTACCTCCGGGGGCAACGCATTTGCGt                                                                                                                                                 >  1:95086/1‑150 (MQ=255)
      ccGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTg                                                                                                                                               >  1:862144/1‑150 (MQ=255)
          gTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGTATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCt                                                                                                                                           <  1:67350/150‑1 (MQ=255)
          gTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCt                                                                                                                                           >  2:67350/1‑150 (MQ=255)
                tCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGc                                                                                                                                                   >  1:1021755/1‑136 (MQ=255)
                tCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGc                                                                                                                                                   <  2:1021755/136‑1 (MQ=255)
                tCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCg                                                                                                                                     >  2:570116/1‑150 (MQ=255)
                           tCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcat                                                                                                                          >  2:1011033/1‑150 (MQ=255)
                                gCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAAtt                                                                                                                     >  1:1124725/1‑150 (MQ=255)
                                     gCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAttt                                                                                                                                                     >  2:602872/1‑113 (MQ=255)
                                     gCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAttt                                                                                                                                                     <  1:602872/113‑1 (MQ=255)
                                        tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTg                                                                                                             <  2:1213127/150‑1 (MQ=255)
                                        tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTg                                                                                                             >  2:818869/1‑150 (MQ=255)
                                           aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCa                                                                                                          <  2:1044779/150‑1 (MQ=255)
                                           aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCa                                                                                                          <  2:862144/150‑1 (MQ=255)
                                           aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGcc                                                                                                                                             >  1:717524/1‑115 (MQ=255)
                                           aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGcc                                                                                                                                             <  2:717524/115‑1 (MQ=255)
                                            aCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCAc                                                                                                         >  1:147312/1‑150 (MQ=255)
                                             ccGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca                                                                                                                >  1:515941/1‑142 (MQ=255)
                                             ccGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca                                                                                                                <  2:515941/142‑1 (MQ=255)
                                              cGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgt                                                                                                       <  1:668653/150‑1 (MQ=255)
                                                     cGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGt                                                                                                                             <  2:374623/121‑1 (MQ=255)
                                                     cGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGt                                                                                                                             >  1:374623/1‑121 (MQ=255)
                                                      gTAACGCTCTGGCACCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTACCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCt                                                                                               <  2:95086/150‑1 (MQ=255)
                                                        aaCGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGa                                                                                             <  1:512315/150‑1 (MQ=255)
                                                                 gCCCCCTAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGTCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAg                                                                                    <  2:644201/150‑1 (MQ=255)
                                                                    cccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGa                                                                                 <  1:550231/150‑1 (MQ=255)
                                                                    cccGAAGCGATTTCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGCCTCTCGTCAGCAATTCATCAGTGCCCCGTTGACGCCTGACACCTTCAGCgc                                                                                 >  2:119915/1‑149 (MQ=255)
                                                                             aTTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTc                                                                        <  2:1124725/150‑1 (MQ=255)
                                                                                    gcCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACACCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCa                                                                 >  1:1221566/1‑150 (MQ=255)
                                                                                         cccGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTgcc                                                            >  1:1260072/1‑150 (MQ=255)
                                                                                                   aaaGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGc                                                  >  1:3812/1‑150 (MQ=255)
                                                                                                          gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACAtt                                                                                       >  2:418755/1‑106 (MQ=255)
                                                                                                          gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACAtt                                                                                       <  1:418755/106‑1 (MQ=255)
                                                                                                          gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCt                                                                           <  2:336972/118‑1 (MQ=255)
                                                                                                          gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAAc                                           >  2:399213/1‑150 (MQ=255)
                                                                                                          gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAAc                                           <  2:1270894/150‑1 (MQ=255)
                                                                                                          gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGACACTGCCt                                                                           >  1:336972/1‑118 (MQ=255)
                                                                                                            tcatcaATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAg                                         >  2:1138736/1‑150 (MQ=255)
                                                                                                                     gcgaaCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGCCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACtt                                <  1:1240882/147‑1 (MQ=255)
                                                                                                                        aaCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACCCCTGACAAATTCAGCCACTGCTTGTCGCAATCATTGGCGTGATGCTGCAACGAACAGGCACCACTTGGc                             <  2:147312/150‑1 (MQ=255)
                                                                                                                           aGCATTCCCTCCTGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtg                          >  2:776972/1‑150 (MQ=255)
                                                                                                                           aGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtg                          >  2:723830/1‑150 (MQ=255)
                                                                                                                            gCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgc                         <  1:818869/150‑1 (MQ=255)
                                                                                                                             cATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTc                                                                        <  2:1057101/102‑1 (MQ=255)
                                                                                                                             cATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTc                                                                        >  1:1057101/1‑102 (MQ=255)
                                                                                                                               ttCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgctgc                      >  1:123051/1‑150 (MQ=255)
                                                                                                                                  ccTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTc                   <  1:1011033/150‑1 (MQ=255)
                                                                                                                                   cTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCt                  >  1:469947/1‑150 (MQ=255)
                                                                                                                                      gaagggCAACGCAAATGCGTTGCCCCTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATcc                                                     <  2:968775/109‑1 (MQ=255)
                                                                                                                                      gaagggCAACGCAAATGCGTTGCCCCTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATcc                                                     >  1:968775/4‑112 (MQ=255)
                                                                                                                                           gCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGtcatc                                                                                                                         <  2:224197/39‑1 (MQ=21)
                                                                                                                                           gCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGtcatc                                                                                                                         >  1:224197/1‑39 (MQ=21)
                                                                                                                                               cGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGt      <  1:399213/150‑1 (MQ=255)
                                                                                                                                                  aTTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAg   <  1:723830/150‑1 (MQ=255)
                                                                                                                                                   tttGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGc  <  1:1138736/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAGC  >  NC_002947/4311211‑4311507

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