Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,503,465 G→T 10.7% 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,4650GT10.7% 166.7 / 9.5 56E90* (GAG→TAG) ‡rplM50S ribosomal protein L13
Reads supporting (aligned to +/- strand):  ref base G (26/24);  new base T (4/2);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

ACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAG  >  NC_002947/1503320‑1503602
                                                                                                                                                 |                                                                                                                                         
aCCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCgag                                                                                                                                         >  1:222643/1‑148 (MQ=255)
 ccAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGaa                                                                                                                                       <  3:8785/149‑1 (MQ=255)
  cAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAg                                                                                                                                      <  2:263721/149‑1 (MQ=255)
   aGAATACACCCCTCACGTTGACACCGGCGACTAATTCGTCGTTATCATCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGt                                                                                                                                     <  5:368319/149‑1 (MQ=255)
         cacCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATc                                                                                                                                <  8:39196/148‑1 (MQ=255)
         cacCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCg                                                                                                                               <  8:160105/149‑1 (MQ=255)
         cacCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCg                                                                                                                               >  6:152159/1‑149 (MQ=255)
            cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtttcct                                                                                                                            <  1:38710/149‑6 (MQ=255)
            cccTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGcc                                                                                                                             >  4:289577/1‑148 (MQ=255)
                 aCGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGGTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGc                                                                                                                        <  7:158401/148‑1 (MQ=255)
                                 cTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATc                                                                                                        >  6:205722/1‑148 (MQ=255)
                                    cATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTg                                                                                                                                   >  6:99240/1‑118 (MQ=255)
                                    cATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTg                                                                                                                                   <  5:99240/118‑1 (MQ=255)
                                    cATCGTCGTTAGCAACGCCGAGCAGGTTCGTGTGACAGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGaaa                                                                                                    <  8:119700/149‑1 (MQ=255)
                                        gtcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGc                                                                                                <  6:244937/149‑1 (MQ=255)
                                         tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCg                                                                                               >  8:246735/1‑149 (MQ=255)
                                         tcgtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCg                                                                                               <  3:289577/149‑1 (MQ=255)
                                           gtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCg                                                                                               <  6:336591/147‑1 (MQ=255)
                                           gtTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGt                                                                                             <  5:152159/149‑1 (MQ=255)
                                                cAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAg                                                                                        >  2:197423/1‑149 (MQ=255)
                                                cAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAg                                                                                        >  3:170276/1‑149 (MQ=255)
                                                  aCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGCGCGTGTTATCGAAACTGCGGTCAAAgg                                                                                       >  7:103956/1‑148 (MQ=255)
                                                  aCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAgg                                                                                       >  6:6774/1‑148 (MQ=255)
                                                  aCGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAgg                                                                                       >  7:403361/1‑148 (MQ=255)
                                                       gAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCAtgct                                                                                 >  3:198661/1‑149 (MQ=255)
                                                             gTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtt                                                                                                                                >  5:1699/1‑95 (MQ=255)
                                                             gTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCTCGAAGTTGAtt                                                                                                                                <  6:1699/96‑2 (MQ=255)
                                                              ttCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCt                                                                             >  1:388310/1‑146 (MQ=255)
                                                                      aCTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTTTGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGc                                                                   <  8:103956/148‑1 (MQ=255)
                                                                      aCTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGc                                                                   <  4:170276/148‑1 (MQ=255)
                                                                          gTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTgg                                                                <  5:167472/147‑1 (MQ=255)
                                                                                aGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGa                                                         >  7:144531/1‑148 (MQ=255)
                                                                                aGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGAc                                                        <  5:6774/149‑1 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGaa                                                                                                                                       >  5:219046/1‑66 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGaa                                                                                                                                       <  6:219046/66‑1 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGaaa                                                                                                    <  5:40091/101‑1 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGaaa                                                                                                    >  6:40091/1‑101 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                              >  4:307335/1‑107 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCgg                                                                                              <  3:307335/107‑1 (MQ=255)
                                                                                    ttCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATg                                                     >  7:130376/1‑148 (MQ=255)
                                                                                        gACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGGTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCg                                                >  7:245051/1‑149 (MQ=255)
                                                                                            aaaTGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAg                                            <  1:197423/149‑1 (MQ=255)
                                                                                                 tactacTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGa                                        >  6:108048/1‑148 (MQ=255)
                                                                                                   ctactCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGaaa                                      <  2:222643/148‑1 (MQ=255)
                                                                                                      ctCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAgtgt                                  <  7:246735/149‑1 (MQ=255)
                                                                                                       tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTg                                                                 >  1:280219/1‑117 (MQ=255)
                                                                                                       tCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTg                                                                 <  2:280219/117‑1 (MQ=255)
                                                                                                         ccACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTAcg                               >  6:163869/1‑149 (MQ=255)
                                                                                                                      ccAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCa                  >  4:163445/1‑149 (MQ=255)
                                                                                                                       cAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAc                 >  4:323255/1‑149 (MQ=255)
                                                                                                                       cAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCAc                 >  7:374206/1‑149 (MQ=255)
                                                                                                                                  aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTg       >  4:235691/1‑148 (MQ=255)
                                                                                                                                  aaGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGc      >  5:157490/1‑149 (MQ=255)
                                                                                                                                    ggAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTc    <  8:374206/149‑1 (MQ=255)
                                                                                                                                      cgaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcag  >  1:7832/3‑149 (MQ=255)
                                                                                                                                       gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcag  >  6:255317/2‑148 (MQ=255)
                                                                                                                                       gaTCAACTTCTCGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTcag  >  2:94997/2‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                         
ACCAGAATACACCCCTCACGTTGACACCGGCGACTACATCGTCGTTATCAACGCCGAGCAGGTTCGTGTGACTGGTGCCAAGTCTTCCGACAAAATGTACTACTCCCACTCCGGTTTCCCAGGCGGTATCAAGGAAATCAACTTCGAGAAGTTGATCGCCAAGGCCCCTGAGCGTGTTATCGAAACTGCGGTCAAAGGCATGCTGCCTAAGAACCCGCTGGGTCGCGACATGTACCGCAAGCTGAAAGTGTACGCGGGTGCTGCTCACCCACACACTGCTCAG  >  NC_002947/1503320‑1503602

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: ‑