Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,659,612 A→T 5.1% intergenic (+28/‑171) ynfD → / → ynfE hypothetical protein/oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,659,6120AT 5.1% 223.2 / 3.4 78intergenic (+28/‑171)ynfD/ynfEhypothetical protein/oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (2/2);  total (39/39)
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.15e-01

CGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATT  >  AP009048/1659544‑1659681
                                                                    |                                                                     
cgcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGAc                                                                      >  1:1571254/1‑70 (MQ=255)
 gcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                     >  1:4592458/1‑70 (MQ=255)
 gcAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                     >  1:260662/1‑70 (MQ=255)
   aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAAt                                                                  <  1:3397337/71‑1 (MQ=255)
   aGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAAt                                                                  <  1:3116916/71‑1 (MQ=255)
     tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                >  1:635609/1‑71 (MQ=255)
     tCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                >  1:4124584/1‑71 (MQ=255)
      cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                               >  1:6769075/1‑71 (MQ=255)
      cTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                               >  1:2762958/1‑71 (MQ=255)
       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                              >  1:6438097/1‑71 (MQ=255)
       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                              >  1:5752596/1‑71 (MQ=255)
       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                              >  1:6447306/1‑71 (MQ=255)
       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                              >  1:3602005/1‑71 (MQ=255)
       tAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCgg                                                              >  1:2946915/1‑71 (MQ=255)
         gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCg                                                               <  1:2303064/68‑1 (MQ=38)
         gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                                 <  1:5692329/66‑1 (MQ=255)
         gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                                 <  1:297128/66‑1 (MQ=255)
         gCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                               <  1:3496815/68‑1 (MQ=255)
          ccAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACa                                                                     <  1:7320326/61‑1 (MQ=255)
           cAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCg                                                               >  1:5551894/1‑66 (MQ=255)
            aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                         >  1:1312980/1‑71 (MQ=255)
            aGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAtt                                                         >  1:7750757/1‑71 (MQ=255)
             ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATc                                                                 >  1:6250750/1‑62 (MQ=255)
             ggATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                        >  1:6264807/1‑71 (MQ=255)
                  gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAAc                                                    <  1:5361663/70‑1 (MQ=255)
                  gcgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                   >  1:6994427/1‑71 (MQ=255)
                   cgcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACa                                                   >  1:2898138/1‑70 (MQ=255)
                    gcgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTg                                                       <  1:2857769/65‑1 (MQ=255)
                     cgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                <  1:6228825/71‑1 (MQ=255)
                     cgcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                <  1:367781/71‑1 (MQ=255)
                      gcCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATac                                                >  1:3141447/1‑70 (MQ=255)
                         tGCCGATCCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaac                                             <  1:4611776/70‑1 (MQ=255)
                         tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                <  1:3853946/51‑1 (MQ=255)
                         tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATcc                                                                <  1:2242572/51‑1 (MQ=255)
                         tGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                            >  1:2764904/1‑71 (MQ=255)
                          gCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAt                                                          >  1:5674882/1‑56 (MQ=255)
                           ccGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                            <  1:2885520/69‑1 (MQ=255)
                           ccGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATacaaca                                            <  1:4965419/69‑1 (MQ=255)
                             gAACCGCAGTAACACATTATCGACTGAACGCCGGATTTGTCATATCCg                                                               >  1:7017804/1‑48 (MQ=25)
                                   cAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGa                                                      >  1:2981791/1‑51 (MQ=255)
                                     gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtt                                <  1:5381865/71‑1 (MQ=255)
                                     gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCtt                                <  1:7088610/71‑1 (MQ=255)
                                     gTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCt                                 >  1:2740863/1‑70 (MQ=255)
                                      tAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGAttt                                                        <  1:748844/46‑1 (MQ=255)
                                        acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                             <  1:3862338/71‑1 (MQ=255)
                                        acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                             <  1:4439356/71‑1 (MQ=255)
                                        acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                              <  1:7410301/70‑1 (MQ=255)
                                        acacATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                              <  1:4177389/70‑1 (MQ=255)
                                        acacATTATCGACTGAACGCCGGATAAGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttat                              >  1:3523984/1‑70 (MQ=255)
                                           caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACAt                                           <  1:4263754/54‑1 (MQ=255)
                                           caTTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                             <  1:7610873/68‑1 (MQ=255)
                                             ttATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATc                                       >  1:4969510/1‑56 (MQ=255)
                                             ttATCGACTGAACGCCGGATATGAAAAATCCGGCGATTTGAACATACAACATAATCCCACCttatta                            >  1:2404060/1‑67 (MQ=255)
                                               aTCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                       <  1:7021548/70‑1 (MQ=255)
                                                tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATa                      <  1:5709111/70‑1 (MQ=255)
                                                tCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATAc                     >  1:3895719/1‑71 (MQ=255)
                                                    cTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCtt                 >  1:4468933/1‑71 (MQ=255)
                                                       aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATaa                                         >  1:5493212/1‑44 (MQ=255)
                                                       aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                <  1:4140405/69‑1 (MQ=255)
                                                       aaCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                <  1:5800203/69‑1 (MQ=255)
                                                        aCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa              <  1:4793747/70‑1 (MQ=255)
                                                        aCGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAt             <  1:2855376/71‑1 (MQ=255)
                                                         cGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTATTCATACCCTTCTAtt            >  1:7247639/1‑71 (MQ=255)
                                                         cGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt            <  1:4912921/71‑1 (MQ=255)
                                                           ccGGATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                          <  1:5161990/55‑1 (MQ=37)
                                                           ccGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGa          <  1:3412591/71‑1 (MQ=255)
                                                            cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                             <  1:5670731/51‑1 (MQ=255)
                                                            cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCttatt                             <  1:600856/51‑1 (MQ=255)
                                                            cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat         <  1:5232603/71‑1 (MQ=255)
                                                            cGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGat         <  1:1922881/71‑1 (MQ=255)
                                                             ggATTTGTCATATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata        >  1:7396876/1‑71 (MQ=255)
                                                             ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGata        >  1:6745855/1‑71 (MQ=255)
                                                             ggATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTAatc        <  1:4418867/71‑2 (MQ=255)
                                                              gATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCAt                       >  1:6160229/1‑55 (MQ=255)
                                                               atatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg      >  1:4827057/1‑71 (MQ=255)
                                                                tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTa              >  1:4234647/1‑62 (MQ=255)
                                                                tatGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATg      >  1:1667073/1‑70 (MQ=255)
                                                                   gACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt   >  1:7170496/1‑70 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                          <  1:3238297/46‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                          <  1:5338410/46‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACt                          <  1:6383817/46‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                <  1:5745082/56‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                <  1:1866492/56‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTc                <  1:7408587/56‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTAtt            <  1:3400718/60‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt  >  1:4017711/1‑70 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt  >  1:6258590/1‑70 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt  <  1:646197/70‑1 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAtt  >  1:499143/1‑70 (MQ=255)
                                                                    aCAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGAt   >  1:2884403/1‑69 (MQ=255)
                                                                    |                                                                     
CGCAGTCTAGCCAGGATGGCGCGCCTGCCGAACCGCAGTAACACATTATCGACTGAACGCCGGATATGACAAATCCGGCGATTTGAACATACAACATAATCCCACCTTATTACTCATACCCTTCTATTGATATGGATT  >  AP009048/1659544‑1659681

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