Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,503,474 G→T 14.3% E90* (GAG→TAG) ‡ rplM → 50S ribosomal protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,503,4740GT14.3% 108.3 / 9.5 42E90* (GAG→TAG) ‡rplM50S ribosomal protein L13
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base T (3/3);  total (22/20)
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.51e-01

ACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAG  >  NC_002947/1503329‑1503614
                                                                                                                                                 |                                                                                                                                            
aCCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGaa                                                                                                                                          >  2:256431/1‑150 (MQ=255)
 ccAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAACTGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAg                                                                                                                                         <  2:583960/150‑1 (MQ=255)
 ccAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAg                                                                                                                                         >  1:386294/1‑150 (MQ=255)
 ccAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAg                                                                                                                                         <  1:81487/150‑1 (MQ=255)
  cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGt                                                                                                                                        >  1:1114139/1‑150 (MQ=255)
  cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGt                                                                                                                                        <  1:251168/150‑1 (MQ=255)
    gAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTg                                                                                                                                      <  2:416889/150‑1 (MQ=255)
     aaTACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGNTGa                                                                                                                                     <  2:704216/150‑1 (MQ=255)
         cacCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGAGAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttc                                                                                                                                 <  2:638693/150‑4 (MQ=255)
               tCACGTTGACACCGGCGACTACATCGTCGCTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAAGCTGTATCAAGGAAATCAACTTCGACAAGTCAATCGCCAAGGc                                                                                                                           >  2:690111/1‑150 (MQ=255)
                   gTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCt                                                                                                                       <  2:386294/150‑1 (MQ=255)
                         acCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCTGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCgt                                                                                                                 <  2:616578/150‑1 (MQ=255)
                                     aTCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCga                                                                                                                                             >  2:792184/1‑110 (MQ=255)
                                     aTCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCga                                                                                                                                             <  1:792184/110‑1 (MQ=255)
                                         tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                                 >  1:348703/1‑150 (MQ=255)
                                         tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                                 <  2:330467/150‑1 (MQ=255)
                                            ttATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCa                                                                                              >  1:832554/1‑150 (MQ=255)
                                               tCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTTATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAg                                                                                           <  1:256431/150‑1 (MQ=255)
                                                       gAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAtgctg                                                                                   >  2:1015195/1‑150 (MQ=255)
                                                        aGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAtgctgc                                                                                  <  2:348703/150‑1 (MQ=255)
                                                                           tGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAt                                                                                       <  1:786846/126‑1 (MQ=255)
                                                                           tGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAt                                                                                       >  2:786846/1‑126 (MQ=255)
                                                                             ccAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTcgcg                                                             >  1:368552/1‑150 (MQ=255)
                                                                              cAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAAATGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGa                                                            <  2:1114139/150‑1 (MQ=255)
                                                                               aaGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGAc                                                           >  2:736304/1‑150 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTa                                                      >  1:775637/1‑150 (MQ=255)
                                                                                        gACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                  >  1:165900/1‑150 (MQ=255)
                                                                                          cAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCaa                                                <  2:775637/150‑1 (MQ=255)
                                                                                           aaaaTGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAg                                               <  1:1015195/150‑1 (MQ=255)
                                                                                               tGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                   <  1:246505/62‑1 (MQ=255)
                                                                                               tGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                   >  2:246505/1‑62 (MQ=255)
                                                                                                       tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                   >  2:387815/1‑54 (MQ=255)
                                                                                                       tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                   <  1:387815/54‑1 (MQ=255)
                                                                                                       tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAc                                   >  1:994856/1‑150 (MQ=255)
                                                                                                         ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAATCAGAAAGTGTAccc                                 >  2:171530/1‑148 (MQ=255)
                                                                                                             tcccttttCCGAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGt                             >  2:276429/6‑150 (MQ=255)
                                                                                                                 gTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCg                                                                                                  <  1:178758/77‑1 (MQ=255)
                                                                                                                 gTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCg                                                                                                  >  2:178758/1‑77 (MQ=255)
                                                                                                                                       gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGt                                                                 <  2:529984/87‑1 (MQ=255)
                                                                                                                                       gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGt                                                                 >  1:529984/2‑88 (MQ=255)
                                                                                                                                       gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagca   >  2:410975/2‑150 (MQ=255)
                                                                                                                                        aTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagcag  >  2:767796/1‑150 (MQ=255)
                                                                                                                                              ttCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGtg                            <  2:567973/118‑1 (MQ=255)
                                                                                                                                              ttCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGtg                            >  1:567973/1‑118 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
ACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAG  >  NC_002947/1503329‑1503614

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