Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,007,421 T→C S60S (TCA→TCG ycbX ← 6‑N‑hydroxylaminopurine resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,007,4210TC80.0% 85.0 / 17.7 40S60S (TCA→TCGycbX6‑N‑hydroxylaminopurine resistance protein
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base C (18/14);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 4.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGGCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCGTGATCATAAAGAT  >  NC_000913/1007345‑1007498
                                                                            |                                                                             
gtgtgGCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAg                                                                           >  1:190221/1‑81 (MQ=255)
 tgtgGCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAg                                                                           <  1:190189/80‑1 (MQ=255)
 tgtgGCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAgg                                                                          <  2:190190/81‑1 (MQ=255)
    ggCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAg                                                                        >  1:190225/1‑80 (MQ=255)
     gCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTa                                                                      <  2:190191/81‑1 (MQ=255)
     gCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGt                                                                       <  1:190192/80‑1 (MQ=255)
       gAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTaaa                                                                    >  2:190231/1‑81 (MQ=255)
       gAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTaaa                                                                    >  1:190226/1‑81 (MQ=255)
         aaTCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAAcc                                                                  <  2:190193/81‑1 (MQ=255)
          aTCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCg                                                                 >  1:190236/1‑81 (MQ=255)
           tCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGt                                                                <  2:190194/81‑1 (MQ=255)
                aaaaCGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCTCAGGCGAAGGAGTAAACCGTACCAt                                                           >  1:190250/1‑81 (MQ=255)
                aaaaCGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCAt                                                           <  2:190195/81‑1 (MQ=255)
                      aaCATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTggg                                                      >  1:190232/1‑80 (MQ=255)
                       aCATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTgggg                                                     >  1:190243/1‑80 (MQ=255)
                       aCATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTgggg                                                     >  1:190240/1‑80 (MQ=255)
                         atatGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGaaa                                                  >  2:190247/1‑81 (MQ=255)
                          tatGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAAc                                                 <  1:190196/81‑1 (MQ=255)
                           atGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACt                                                >  1:190235/1‑81 (MQ=255)
                            tGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACt                                                <  2:190197/80‑1 (MQ=255)
                            tGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACt                                                <  1:190198/80‑1 (MQ=255)
                             gcgcTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTgg                                              >  2:190237/1‑81 (MQ=255)
                               gcTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCg                                            <  1:190199/81‑1 (MQ=255)
                                   cTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGc                                         <  1:190201/80‑1 (MQ=255)
                                   cTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCg                                        <  1:190200/81‑1 (MQ=255)
                                          cTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTAAACCGCACCATCTGGGGAAACTGGCGAGCGGtaata                                  >  2:190272/1‑80 (MQ=255)
                                          cTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGtaata                                  >  1:190239/1‑80 (MQ=255)
                                          cTGGTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGtaata                                  >  2:190241/1‑80 (MQ=255)
                                             gTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAAc                               <  1:190202/80‑1 (MQ=255)
                                             gTGCGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACg                              >  1:190242/1‑81 (MQ=255)
                                                cGGTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTAc                           <  2:190203/81‑1 (MQ=255)
                                                 ggTGAGATGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTAcc                          <  2:190204/81‑1 (MQ=255)
                                                      gaTGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAg                     <  2:190205/81‑1 (MQ=255)
                                                       aTGCAAGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAgg                    <  2:190206/81‑1 (MQ=255)
                                                           aaGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTcc                <  1:190207/81‑1 (MQ=255)
                                                            aGCCATCATGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCg               >  1:190248/1‑81 (MQ=255)
                                                                atcatGCACAGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCGTGAt           <  1:190208/81‑1 (MQ=255)
                                                                atcatGCACAGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCGTGAt           >  1:190259/1‑81 (MQ=255)
                                                                       acaGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCGTGATCATAAAg    <  2:190209/81‑1 (MQ=255)
                                                                         aGGCGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCGTGATCATAAAGAt  >  2:190245/1‑81 (MQ=255)
                                                                            |                                                                             
GTGTGGCGAAATCAGCAAAACGAACATATGCGCTACTGCCATCTGGTGCGGTGAGATGCAAGCCATCATGCACAGGTGAAGGAGTAAACCGTACCATCTGGGGAAACTGGCGAGCGGTAATAAACGTACCGTCAGGTTCCGTGATCATAAAGAT  >  NC_000913/1007345‑1007498

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