Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,008,134 C→T 14.3% intergenic (‑112/‑132) ycbX ← / → ycbY predicted 2Fe‑2S cluster‑containing protein/predicted methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,008,1340CT14.3% 130.9 / 12.2 49intergenic (‑112/‑132)ycbX/ycbYpredicted 2Fe‑2S cluster‑containing protein/predicted methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (24/18);  new base T (4/3);  total (28/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

ATGACGCGCTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCG  >  AP009048/1008064‑1008196
                                                                      |                                                              
atgacGCGCTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTc                                                                <  1:3958621/71‑1 (MQ=255)
 tgacGCGCTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCa                                                               <  1:1686882/71‑1 (MQ=255)
   acGCGCTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTTAGt                                                             <  1:744039/71‑1 (MQ=255)
       gcTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGttt                                                           >  1:2382601/1‑69 (MQ=255)
       gcTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTgg                                                         >  1:2834007/1‑71 (MQ=255)
           attaGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGccc                                                      <  1:2824807/70‑1 (MQ=255)
            ttaGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCt                                                     <  1:1971466/70‑1 (MQ=255)
              aGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTggg                                                  >  1:2052843/1‑71 (MQ=255)
               gCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAg                                                              <  1:2564096/58‑1 (MQ=255)
                   taATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGt                                                             <  1:730793/55‑1 (MQ=255)
                        cgcAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTa                                               >  1:1034021/1‑64 (MQ=255)
                          caacaaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGt                                                 <  1:3158675/60‑1 (MQ=255)
                           aacaaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGaaa                                      <  1:2545960/70‑1 (MQ=255)
                            acaaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGtt                                                            >  1:1503744/1‑47 (MQ=255)
                            acaaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTg                                                          <  1:1862124/49‑1 (MQ=255)
                            acaaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTc                                          >  1:3142928/1‑65 (MQ=255)
                             caaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTTAGTTTGGCCCTGGGTTATGCTCGAAATcc                                   >  1:3829435/1‑71 (MQ=255)
                             caaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATcc                                   >  1:869175/1‑71 (MQ=255)
                              aaTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCt                                  <  1:1977395/71‑1 (MQ=255)
                               aTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTTAGTTTGGCCCTGGGTTATGCTCGAAATCCTc                                 <  1:2875169/71‑1 (MQ=255)
                               aTTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCt                                  >  1:2119541/1‑70 (MQ=255)
                                ttttCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCa                                >  1:3935212/1‑71 (MQ=255)
                                     ttAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGa                                        <  1:3911432/58‑1 (MQ=255)
                                       aGCGCCTGTCCCATTAGGCTATTTTATTTGTTAGTTTGGCCCTGGGTTa                                               >  1:3215538/1‑49 (MQ=255)
                                        gcgcCTGTCCCATTAGGCTATTTTATTTGTTAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACtt                         <  1:3391616/70‑1 (MQ=255)
                                        gcgcCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACtt                         <  1:480070/70‑1 (MQ=255)
                                         cgcCTGTCCCATTAGGCTATTTTATTTGTTAGTTTGGcc                                                       >  1:2437247/1‑39 (MQ=255)
                                         cgcCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCt                                           >  1:3564239/1‑51 (MQ=255)
                                           ccTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGa                                        >  1:2911969/1‑52 (MQ=255)
                                           ccTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTgtgt                     <  1:4285775/71‑1 (MQ=255)
                                            cTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCt                                  <  1:329896/57‑1 (MQ=255)
                                            cTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTa                    >  1:4281646/1‑71 (MQ=255)
                                             tGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTAc                   >  1:2459482/1‑71 (MQ=255)
                                             tGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTAc                   >  1:1984108/1‑71 (MQ=255)
                                               tCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGa                                        >  1:1492523/1‑48 (MQ=255)
                                               tCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGa                                        >  1:1979089/1‑48 (MQ=255)
                                               tCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACg                  <  1:3518310/70‑1 (MQ=255)
                                                cccATTAGGCTATTTTATTTGTCAGTTTGGCCCTg                                                    >  1:2264262/1‑35 (MQ=255)
                                                 ccATTAGGCTATTTTTTTTGTCAGTTTGGCCCTGGGTTATGCTCGaaa                                      >  1:715284/1‑48 (MQ=255)
                                                 ccATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGaaa                                      >  1:2585565/1‑48 (MQ=255)
                                                   aTTAGGCTATTTTATTTGTTAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTgtgt                     >  1:4308692/1‑63 (MQ=255)
                                                   aTTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGaaa                                      <  1:4175231/46‑1 (MQ=255)
                                                      aGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCa                                >  1:262286/1‑49 (MQ=255)
                                                         cTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTtct       <  1:4551545/71‑1 (MQ=255)
                                                           attttatttGTCAGTTTGGCCCTGGGTTATGCTCGaaa                                      <  1:3127992/38‑1 (MQ=255)
                                                            ttttatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCcg    >  1:2910172/1‑71 (MQ=255)
                                                             tttatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTAc                           <  1:277807/47‑1 (MQ=255)
                                                             tttatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCg             >  1:1062503/1‑61 (MQ=255)
                                                             tttatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCcgc   >  1:2417405/1‑71 (MQ=255)
                                                              ttatttGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCcgcg  >  1:1826602/1‑71 (MQ=255)
                                                                      |                                                              
ATGACGCGCTTATTAGCTATAATGCGCAACAATTTTCTTAGCGCCTGTCCCATTAGGCTATTTTATTTGTCAGTTTGGCCCTGGGTTATGCTCGAAATCCTCACGTACTTGTGTACGCTCCGGTTTCTCCGCG  >  AP009048/1008064‑1008196

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