Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,311,315 T→A 10.9% 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,3150TA10.9% 123.7 / 5.5 46F12Y (TTT→TAT) ‡PP_3784hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/19);  new base A (3/2);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.12e-01

CCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAG  >  NC_002947/4311173‑4311462
                                                                                                                                              |                                                                                                                                                   
ccGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGtt                                                                                                                                               <  5:169137/149‑1 (MQ=255)
   cgTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGcc                                                                                                                                            <  2:140562/149‑1 (MQ=255)
       cgcTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTc                                                                                                                                        <  1:244551/149‑1 (MQ=255)
        gcTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCg                                                                                                                                       >  6:284024/1‑149 (MQ=255)
            gCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGcccccggagg                                                                                                                                    <  1:111394/148‑8 (MQ=255)
              aaaCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGct                                                                                                                                 <  7:141525/149‑1 (MQ=255)
                  cACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCg                                                                                                                             >  3:266970/1‑149 (MQ=255)
                  cACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGAGGGCAACGCAAATGCGtt                                                                                                                                               <  2:211893/131‑1 (MQ=255)
                  cACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGAGGGCAACGCAAATGCGtt                                                                                                                                               >  1:211893/1‑131 (MQ=255)
                  cACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCg                                                                                                                             >  4:225393/1‑148 (MQ=255)
                                  tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             <  6:164273/149‑1 (MQ=255)
                                  tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             <  4:142103/149‑1 (MQ=255)
                                  tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTg                                                                                                            <  3:225393/149‑1 (MQ=255)
                                     aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGc                                                                                                           >  1:82056/1‑148 (MQ=255)
                                        cGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCCCg                                                                                                       >  4:254362/1‑149 (MQ=255)
                                          aaTTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCAGCAATTCATCAGTGCCACgtg                                                                                                     >  5:383976/1‑149 (MQ=255)
                                           aTTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgtgt                                                                                                    <  5:284024/149‑1 (MQ=255)
                                                gTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             <  3:67793/135‑1 (MQ=255)
                                                gTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGt                                                                                                             >  4:67793/1‑135 (MQ=255)
                                                       tctGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGAcaca                                                                                        >  1:362948/1‑149 (MQ=255)
                                                                  aaTCGATTGCGCGCCCTCCCGTCAGTGAAACGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTg                                                                             >  1:32303/1‑149 (MQ=255)
                                                                          gcgcgcCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCg                                                                      <  4:266970/148‑1 (MQ=255)
                                                                                ccTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCAt                                                               >  3:69187/1‑149 (MQ=255)
                                                                                  tCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTg                                                             >  8:105404/1‑149 (MQ=255)
                                                                                             aaaGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCa                                                   <  4:69187/148‑1 (MQ=255)
                                                                                               aGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGc                                                 >  3:360587/1‑148 (MQ=255)
                                                                                                ggAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCaa                                               >  5:306645/1‑149 (MQ=255)
                                                                                                 gAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAAc                                              >  6:352520/1‑149 (MQ=255)
                                                                                                         tcaATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGc                                      >  2:143830/1‑149 (MQ=255)
                                                                                                                   aCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGac                                                                                           <  4:193993/86‑1 (MQ=255)
                                                                                                                   aCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGac                                                                                           >  3:193993/1‑86 (MQ=255)
                                                                                                                     aGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGt                          >  2:341980/1‑149 (MQ=255)
                                                                                                                      gCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtg                         <  7:105404/149‑1 (MQ=255)
                                                                                                                       cATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgc                        <  6:306645/149‑1 (MQ=255)
                                                                                                                        aTTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgct                       <  4:360587/149‑1 (MQ=255)
                                                                                                                           cccTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGtgctgc                     <  3:254362/148‑1 (MQ=255)
                                                                                                                            ccTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGt                   >  5:229188/1‑149 (MQ=255)
                                                                                                                                gaagggCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTg               >  4:107973/4‑149 (MQ=255)
                                                                                                                                     gCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTg          >  1:19472/1‑149 (MQ=255)
                                                                                                                                       aaCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg         >  2:307523/1‑148 (MQ=255)
                                                                                                                                       aaCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTgg         >  2:263021/1‑148 (MQ=255)
                                                                                                                                       aaCGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCATCTTGTCGTgg         >  7:11632/1‑148 (MQ=255)
                                                                                                                                          gCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTTCGTCTTGTCGTGGCAGt     <  3:339731/149‑1 (MQ=255)
                                                                                                                                          gCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGt     <  2:362948/149‑1 (MQ=255)
                                                                                                                                          gCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGt     >  1:281767/1‑149 (MQ=255)
                                                                                                                                            aTTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCATCTTGTCGTGGCAGTGa   <  8:11632/149‑1 (MQ=255)
                                                                                                                                             tttGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAg  <  7:106819/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGCAGTGAG  >  NC_002947/4311173‑4311462

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: