Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,031,923 T→A 10.0% L5* (TTG→TAG) ‡ PP_0889 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,031,9230TA10.0% 112.0 / 4.6 40L5* (TTG→TAG) ‡PP_0889hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (15/21);  new base A (2/2);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTC  >  NC_002947/1031775‑1032062
                                                                                                                                                    |                                                                                                                                           
aaTATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTg                                                                                                                                            >  1:90400/1‑150 (MQ=255)
    tctcGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGCTGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCcg                                                                                                                                        <  2:553296/150‑1 (MQ=255)
          ccACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgctca                                                                                                                                  <  1:544220/150‑5 (MQ=255)
            aCCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATTGATGATGTTGGCCTGTTCGTGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGc                                                                                                                                >  2:178223/1‑150 (MQ=255)
            aCCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATTGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGc                                                                                                                                >  2:178235/1‑150 (MQ=255)
                       gcTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGt                                                                                                                     >  2:712725/1‑150 (MQ=255)
                            ccATGCACCGACGGACGTTGCAGCGCAACCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCaa                                                                                                                <  1:642333/149‑1 (MQ=255)
                             cATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGCTGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAAt                                                                                                               <  2:460446/150‑1 (MQ=255)
                                gCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCaa                                                                                                            <  1:470589/150‑1 (MQ=255)
                                  aCCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAgg                                                                                                          >  1:833331/1‑150 (MQ=255)
                                      acggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                             <  2:677353/127‑1 (MQ=255)
                                      acggacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGc                                                                                                                             >  1:677353/1‑127 (MQ=255)
                                         gacgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACAGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTc                                                                                                   >  1:791949/1‑150 (MQ=255)
                                           cgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgc                                                                                                                                     <  2:181839/114‑2 (MQ=255)
                                           cgTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCAAGCCcgcgc                                                                                                                                     >  1:181839/1‑113 (MQ=255)
                                                cAGCGCAAGCTGGAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCATTCATCTGTGAGCGCGGTCTTGCCCGGGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                            <  2:130296/150‑1 (MQ=255)
                                                cAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAg                                                                                            <  2:791949/150‑1 (MQ=255)
                                                                        gTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTg                                                                    >  1:302212/1‑150 (MQ=255)
                                                                             gcgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGc                                                               <  2:90400/150‑1 (MQ=255)
                                                                              cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                              <  2:262591/150‑1 (MQ=255)
                                                                              cgcTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCa                                                              <  2:709449/150‑1 (MQ=255)
                                                                                         tATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               >  1:572368/1‑74 (MQ=255)
                                                                                         tATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               <  2:572368/74‑1 (MQ=255)
                                                                                               tgatgTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGa                                             <  2:302212/150‑1 (MQ=255)
                                                                                                    ttGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACt                                        <  1:712725/150‑1 (MQ=255)
                                                                                                       gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTGGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  2:256566/150‑1 (MQ=255)
                                                                                                       gCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCCTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  2:200368/150‑1 (MQ=255)
                                                                                                       gCCTGTTCGCGGGCGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGc                                     <  1:168038/150‑1 (MQ=255)
                                                                                                                 gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               <  2:389841/50‑1 (MQ=255)
                                                                                                                 gggggAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCt                                                                                                                               >  1:389841/1‑50 (MQ=255)
                                                                                                                           gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                 >  1:579230/1‑150 (MQ=255)
                                                                                                                           gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                 <  1:178223/150‑1 (MQ=255)
                                                                                                                           gCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAAc                 <  1:178235/150‑1 (MQ=255)
                                                                                                                               cTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGAAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGAcc             >  1:532362/1‑150 (MQ=255)
                                                                                                                                tGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTCCAACGCACCACCAATAAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGTCTGACCCGCAGTGGCAGAACGAccc            >  1:149986/1‑150 (MQ=255)
                                                                                                                                     agcccttcgcgGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCa       >  2:304499/8‑150 (MQ=255)
                                                                                                                                        tgAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             <  2:747122/93‑1 (MQ=255)
                                                                                                                                        tgAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAg                                                             >  1:747122/1‑93 (MQ=255)
                                                                                                                                         cttcgcgGGCAAGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCt   >  1:639575/4‑150 (MQ=255)
                                                                                                                                          ttcgcgGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTc  <  1:304499/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
AATATCTCGGCCACCGCACGGGCGCTGGGCATGCACCGACGGACGTTGCAGCGCAAGCTGCAGAAGCGGCCGGTTCGGCGCTGATTCCTTATCGATGATGTTGGCCTGTTCGCGGGGGAATCCGCTCCTGCAACTGTGAGCGCGGGCTTGCCCGCGAAGGGCTGCAAAGCAGTCCCAATCAAGGCGAGCTCGATGAAGGTACTGTTGGTCATCCTTGCGCTGGTTGCAGGTAGCGCCAAAGCTGAAGACTGGCCTGACCCGCAGTGGCAGAACGACCCTGCCACCCTC  >  NC_002947/1031775‑1032062

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