Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,503,475 A→C 13.5% E90A (GAG→GCG) ‡ rplM → 50S ribosomal protein L13

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,503,4750AC13.5% 119.6 / 14.3 52E90A (GAG→GCG) ‡rplM50S ribosomal protein L13
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base C (3/4);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

AACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGC  >  NC_002947/1503328‑1503615
                                                                                                                                                   |                                                                                                                                            
aaCCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCgaga                                                                                                                                            >  2:334323/1‑150 (MQ=255)
 aCCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGaa                                                                                                                                           <  2:1095345/150‑1 (MQ=255)
    aGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGtt                                                                                                                                        >  1:129196/1‑150 (MQ=255)
      aaTACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGa                                                                                                                                      <  1:245987/150‑1 (MQ=255)
        tACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                    <  2:129196/150‑1 (MQ=255)
           acCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGa                                                                                                                                      >  1:382505/1‑145 (MQ=255)
           acCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGa                                                                                                                                      <  2:382505/145‑1 (MQ=255)
           acCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGcc                                                                                                                                 >  2:946180/1‑150 (MQ=255)
             cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCCCTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCaa                                                                                                                               <  2:1175786/150‑1 (MQ=255)
             cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcctt                                                                                                                               <  2:509143/150‑7 (MQ=255)
             cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcctt                                                                                                                               <  2:1112593/150‑7 (MQ=255)
             cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCaa                                                                                                                               >  2:202210/1‑150 (MQ=255)
              ccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAg                                                                                                                              >  2:276595/1‑150 (MQ=255)
                                  cTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGa                                                                                                          >  1:582171/1‑150 (MQ=255)
                                         gtcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCga                                                                                                                                              >  2:468669/1‑107 (MQ=255)
                                         gtcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCga                                                                                                                                              <  1:468669/107‑1 (MQ=255)
                                          tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                                  >  2:972857/1‑150 (MQ=255)
                                             ttATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCa                                                                                               >  1:436572/1‑150 (MQ=255)
                                               aTCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAgg                                                                                                                             >  1:672275/1‑118 (MQ=255)
                                               aTCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAgg                                                                                                                             <  2:672275/118‑1 (MQ=255)
                                                  aaCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGc                                                                                          <  1:939926/150‑1 (MQ=255)
                                                     gCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAtg                                                                                       >  2:44545/1‑150 (MQ=255)
                                                            atgTTCGTGTGACTGGTGACAAGTCTTCCGACAAAATGTACTACTACCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTGCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGACTa                                                                                <  2:149603/148‑1 (MQ=255)
                                                              gTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcc                                                                                                                                 >  1:1150521/1‑95 (MQ=255)
                                                              gTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcc                                                                                                                                 <  2:1150521/99‑5 (MQ=255)
                                                              gTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTAATACCAAACAGGTTTCCCAGGCGGTATCAAGTAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAg                                                                              >  1:328817/1‑150 (MQ=255)
                                                                 cGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAAc                                                                           <  1:276595/150‑1 (MQ=255)
                                                                       aCTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTg                                                                     <  1:334323/150‑1 (MQ=255)
                                                                        cTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTgg                                                                    >  2:667356/1‑150 (MQ=255)
                                                                         tGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTggg                                                                   <  1:44545/150‑1 (MQ=255)
                                                                             gCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTcgc                                                               <  1:946180/150‑1 (MQ=255)
                                                                              ccAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAAccc                                                                         >  1:678373/1‑139 (MQ=255)
                                                                              ccAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAAccc                                                                         <  2:678373/139‑1 (MQ=255)
                                                                              ccAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTggg                                                                   >  2:1158799/1‑145 (MQ=255)
                                                                              ccAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTggg                                                                   <  1:1158799/145‑1 (MQ=255)
                                                                                         gACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                   >  1:397187/1‑150 (MQ=255)
                                                                                         gACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGc                                                   >  1:182082/1‑150 (MQ=255)
                                                                                          aCAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGc                                                                                          >  2:966698/1‑110 (MQ=255)
                                                                                          aCAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGc                                                                                          <  1:966698/110‑1 (MQ=255)
                                                                                                 gTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTAt                                                                                                             <  1:71097/84‑1 (MQ=255)
                                                                                                 gTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTAt                                                                                                             >  2:71097/1‑84 (MQ=255)
                                                                                                    ctactCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAgt                                        >  2:1181451/1‑150 (MQ=255)
                                                                                                       ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCg                                                                        <  1:939874/115‑1 (MQ=255)
                                                                                                       ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCg                                                                        >  2:939874/1‑115 (MQ=255)
                                                                                                        tCCCACTCCGGTTTCCCAGGCGGTATCAAGGTAATAGATTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAc                                    <  2:400888/150‑1 (MQ=255)
                                                                                                        tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAc                                    <  1:939434/150‑1 (MQ=255)
                                                                                                          ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcgc                                  >  2:291134/1‑150 (MQ=255)
                                                                                                          ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTGCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcga                                  >  1:400888/1‑149 (MQ=255)
                                                                                                           cACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGCTCGCGACATGTACCGCAAGCTGAAAGTGTAcgcg                                 >  1:503957/1‑150 (MQ=255)
                                                                                                                        cAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAcc                    <  1:972857/150‑1 (MQ=255)
                                                                                                                                   aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCt                                                                      <  1:843949/89‑1 (MQ=255)
                                                                                                                                   aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCt                                                                      >  2:843949/1‑89 (MQ=255)
                                                                                                                                    aGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTc        <  1:668574/150‑1 (MQ=255)
                                                                                                                                        gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGGGGGCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcagca    >  2:834790/2‑150 (MQ=255)
                                                                                                                                          tCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCagcagc  >  2:821650/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                            
AACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAGCAGC  >  NC_002947/1503328‑1503615

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: