Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,311,358 T→A 11.3% F12I (TTT→ATT) ‡ PP_3784 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,311,3580TA11.3% 95.0 / 2.9 34F12I (TTT→ATT) ‡PP_3784hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (17/13);  new base A (2/2);  total (19/15)
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.52e-01

GGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGC  >  NC_002947/4311209‑4311500
                                                                                                                                                     |                                                                                                                                              
ggACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAt                                                                                                                                                <  2:351787/150‑1 (MQ=255)
      aaCCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGt                                                                                                                                          >  1:63395/1‑150 (MQ=255)
         cgcgTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCAAATGCGTTGc                                                                                                                                       <  1:206062/150‑1 (MQ=255)
                      aaaCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGctc                                                                                                                          >  1:773720/1‑150 (MQ=255)
                       aaCCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGctct                                                                                                                         >  1:688384/1‑150 (MQ=255)
                                      ggCGTTTAACCGAATTCCGAAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatc                                                                                                          >  2:32607/1‑150 (MQ=255)
                                       gCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATtcatca                                                                                                         <  2:688384/150‑1 (MQ=255)
                                          tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGcc                                                                                                                                      >  1:142933/1‑118 (MQ=255)
                                          tttAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGcc                                                                                                                                      <  2:142933/118‑1 (MQ=255)
                                             aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGTAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAAtt                                                                                                              >  2:340360/1‑139 (MQ=255)
                                             aaCCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGTAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAAtt                                                                                                              <  1:340360/139‑1 (MQ=255)
                                                                       ccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcat                                                                                                                   <  1:562061/108‑1 (MQ=255)
                                                                       ccGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGtcat                                                                                                                   >  2:562061/1‑108 (MQ=255)
                                                                        cGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCa                                                                                                   >  1:508125/1‑123 (MQ=255)
                                                                        cGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCa                                                                                                   <  2:508125/123‑1 (MQ=255)
                                                                               aTTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCg                                                                                                                                           >  1:102928/1‑76 (MQ=255)
                                                                               aTTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCg                                                                                                                                           <  2:102928/76‑1 (MQ=255)
                                                                               aTTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTc                                                                 >  2:23668/1‑150 (MQ=255)
                                                                                    gcgcCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAAt                                                            <  2:76117/150‑1 (MQ=255)
                                                                                        ccTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCAtt                                                        >  2:103385/1‑150 (MQ=255)
                                                                                             cGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTCGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTgccgc                                                   >  2:666654/1‑150 (MQ=255)
                                                                                                 aGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATc                                               <  1:23668/150‑1 (MQ=255)
                                                                                                  gTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGc                                                                                          <  2:306658/106‑1 (MQ=255)
                                                                                                  gTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGc                                                                                          >  1:306658/1‑106 (MQ=255)
                                                                                                            gATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAAc                                    <  2:773720/150‑1 (MQ=255)
                                                                                                             atcatcAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACa                                   <  2:397166/150‑1 (MQ=255)
                                                                                                                       gCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGaca                                                                                    >  2:622108/1‑91 (MQ=255)
                                                                                                                       gCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGaca                                                                                    <  1:622108/91‑1 (MQ=255)
                                                                                                                       gCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTg                                                                   <  2:436286/108‑1 (MQ=255)
                                                                                                                       gCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTg                                                                   >  1:436286/1‑108 (MQ=255)
                                                                                                                           aCAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCg                     >  1:415729/1‑150 (MQ=255)
                                                                                                                                           gggCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgt                                                                                                <  2:522094/59‑1 (MQ=37)
                                                                                                                                           gggCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACgt                                                                                                >  1:522094/1‑59 (MQ=37)
                                                                                                                                            ggCAACGCAAATGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTATGCGACTGCCTGTCGCAATCATTGCCACCTTACCGCGGCGAACGGTCACGACGTGGCGTGCTGCGTCTTGTCGTg    >  1:178054/1‑150 (MQ=255)
                                                                                                                                              cAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGc  >  1:157455/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GGACACAACCGCGTGCGCTCGCAAACCACTCGTGGCAAGGCGTTTAACCGAATTCCGTAACGCTCTGGCCCCCGAATCGATTGCGCGCCCTCCCGTCAGTGAAAGGAAGATCATCAATGGCTAACAGCATTCCCTCCGGGGGCAACGCATTTGCGTTGCCCTTCGCCGGCTCTCGTCATCAATTCATCAGTGCCACGTGTACGCCTGACACATTCAGCGACTGCCTGTCGCAATCATTGCCGCCATCCAGCAACGAACAGGCACCACTTGGCGTGCTGCGTCTTGTCGTGGC  >  NC_002947/4311209‑4311500

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: