Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,812,772 G→A 11.5% intergenic (+23/‑87) rplT → / → pheS 50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,812,7720GA11.5% 131.9 / 7.6 52intergenic (+23/‑87)rplT/pheS50S ribosomal protein L20/phenylalanine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (23/23);  new base A (3/3);  total (26/26)
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.56e-01

GCTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAG  >  NC_002947/2812633‑2812916
                                                                                                                                           |                                                                                                                                                
gCTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAgcgc                                                                                                                                        <  2:1124977/150‑1 (MQ=255)
  tGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTa                                                                                                                                            <  2:1005843/144‑1 (MQ=255)
  tGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTa                                                                                                                                            >  1:1005843/1‑144 (MQ=255)
       tGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCGCCCGGCAGCGCTAGCACCGCCggg                                                                                                                                 <  1:32490/150‑1 (MQ=255)
        gAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGc                                                                                                                                >  2:672438/1‑150 (MQ=255)
         aaaaaGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCa                                                                                                                               >  2:1090302/1‑150 (MQ=255)
          aaaaGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCAt                                                                                                                              >  1:84867/1‑150 (MQ=255)
                cTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGTGGTGCTAGCGCTGCCGGGTATTGAACg                                                                                                                        <  1:162411/150‑1 (MQ=255)
                    gatcgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCAt                                                                                                                    >  2:543075/1‑150 (MQ=255)
                       cgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAACTAGCGAAGTGTGCGGCGATTGTCGAGAAAGCTCAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGa                                                                                                                 <  1:607390/150‑1 (MQ=255)
                       cgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCGCCCGGCAGCGCTAGCACCGCCgggtga                                                                                                                              <  2:480206/137‑4 (MQ=255)
                       cgaAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCGCCCGGCAGCGCTAGCACCGCCgggtga                                                                                                                              >  1:480206/1‑134 (MQ=255)
                          aaTCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAg                                                                                                              <  1:312156/150‑1 (MQ=255)
                           aTCGACCGTAAGGTTCTGGCTGAACTGGTAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCGGTTCTGGCTTAAGTACCCGCGACCTGCACCGCGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAgg                                                                                                             <  1:543075/150‑1 (MQ=255)
                                    aaGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGcc                                                                                                    <  1:679127/150‑1 (MQ=255)
                                                      gCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTa                                                                                  <  1:544723/150‑1 (MQ=255)
                                                        aGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTAtt                                                                                >  1:60011/1‑150 (MQ=255)
                                                           gAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTc                                                                             <  1:113438/150‑1 (MQ=255)
                                                                      gcgcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCTACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACCTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGt                                                                  <  2:460762/148‑1 (MQ=255)
                                                                        gcgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCt                                                                <  2:839012/150‑1 (MQ=255)
                                                                         cgTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGACCCTTATATGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTg                                                               <  1:91593/150‑1 (MQ=255)
                                                                                          gagaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGa                                              >  1:1022674/1‑150 (MQ=255)
                                                                                           agaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgggtg                                                                                                                               <  2:372700/68‑3 (MQ=21)
                                                                                           agaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCAGCGCTAGCACCGCCgggtg                                                                                                                               >  1:372700/1‑66 (MQ=25)
                                                                                            gaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCcggcgg                                                                                                                                                 <  1:696263/49‑1 (MQ=255)
                                                                                            gaAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCcggcgg                                                                                                                                                 >  2:696263/1‑49 (MQ=255)
                                                                                             aaaGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTa                                                             >  2:947273/1‑132 (MQ=255)
                                                                                             aaaGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTa                                                             <  1:947273/132‑1 (MQ=255)
                                                                                             aaaGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATgc                                           >  2:850751/1‑150 (MQ=255)
                                                                                                gCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCt                                        >  2:459851/1‑150 (MQ=255)
                                                                                                 cTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTg                                       <  1:1090302/150‑1 (MQ=255)
                                                                                                  tAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTgg                                      >  1:698213/1‑150 (MQ=255)
                                                                                                   aaaGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTcccc                                                                                    >  2:1082828/1‑103 (MQ=38)
                                                                                                   aaaGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCTTTGAAGGCTCTTcccc                                                                                    <  1:1082828/103‑1 (MQ=255)
                                                                                                       ccGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTcc                                 <  1:459851/150‑1 (MQ=255)
                                                                                                          ttCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCaa                              >  2:647525/1‑150 (MQ=255)
                                                                                                               gCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCt                         >  1:761099/1‑150 (MQ=255)
                                                                                                                cTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCCTCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGCTCTCCCAAGACCtc                        >  2:229691/1‑149 (MQ=255)
                                                                                                                 ttAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTag                       >  1:531838/1‑150 (MQ=255)
                                                                                                                  tAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                      <  2:84867/150‑1 (MQ=255)
                                                                                                                  tAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTaga                      <  1:850751/150‑1 (MQ=255)
                                                                                                                         ccGCGACAATCACCCGGTGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg               <  1:915707/150‑1 (MQ=255)
                                                                                                                         ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTg               >  2:416193/1‑150 (MQ=255)
                                                                                                                         ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCAATAGGGGAAGAGCTTTGAAGGCTCTTCCCCTAtcga                                                                              >  1:95822/1‑84 (MQ=25)
                                                                                                                         ccGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCAATAGGGGAAGAGCTTTGAAGGCTCTTCCCCTAtcga                                                                              <  2:95822/87‑4 (MQ=21)
                                                                                                                           gcgACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGa             >  2:142085/1‑150 (MQ=255)
                                                                                                                             gACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAc           <  1:416193/150‑1 (MQ=255)
                                                                                                                               cAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAAcgc         >  2:731293/1‑150 (MQ=255)
                                                                                                                                  atgcccGGCAGCGCTAGCACCGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCGCCTATTTTCGTATCTGGAGTCTGTACATGGAACAGCGGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCt      >  1:1239137/4‑150 (MQ=255)
                                                                                                                                    aCCCGGCGGTGCTAGCGCTTCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCCGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGa    >  1:1208103/1‑150 (MQ=255)
                                                                                                                                     cccGGCGGTGCTAGCGCTTCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGaa   <  2:1208103/150‑1 (MQ=255)
                                                                                                                                      ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAggg                                                                                                            >  2:160027/1‑44 (MQ=255)
                                                                                                                                      ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAggg                                                                                                            <  1:160027/44‑1 (MQ=255)
                                                                                                                                      ccGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAg  >  1:182839/1‑150 (MQ=255)
                                                                                                                                           |                                                                                                                                                
GCTGGCCTGAAAAAGGCTTCGATCGAAATCGACCGTAAGGTTCTGGCTGATCTGGCAGTGAACGAAAAAGCGGCGTTTGCTGCGATTGTCGAGAAAGCTAAAGCCGTTCTGGCTTAAGTACCCGCGACAATCACCCGGCGGTGCTAGCGCTGCCGGGCATTGAACGTCATCGATAGGGGAAGAGCCTTCAAAGCTCTTCCCCTATTTTCGTATCTGGAGTCTGTACATGGAAAACCTGGATGCGCTGGTCTCCCAAGCCCTAGAGGCCGTGGAACGCGCTGAAG  >  NC_002947/2812633‑2812916

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