Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,872,074 T→C 8.0% intergenic (+122/‑25) yeaH → / → yeaI conserved hypothetical protein/predicted diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,872,0740TC 8.0% 136.2 / 4.2 50intergenic (+122/‑25)yeaH/yeaIconserved hypothetical protein/predicted diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (24/22);  new base C (2/2);  total (26/24)
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.46e-01

ATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATT  >  AP009048/1872005‑1872133
                                                                     |                                                           
aTCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                            >  1:2872004/1‑71 (MQ=255)
  ctcATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                                           >  1:4521287/1‑70 (MQ=255)
     aTTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGCTCCa                                                         <  1:4512818/69‑1 (MQ=255)
      tttATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGa                                                       >  1:4225803/1‑70 (MQ=255)
       ttattaCATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaa                                                      >  1:438578/1‑70 (MQ=255)
        tattaCATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaaaa                                                    >  1:4013548/1‑71 (MQ=255)
          ttaCATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAAc                                                   <  1:1141602/70‑1 (MQ=255)
             catacatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                                <  1:2875700/70‑1 (MQ=255)
             catacatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                                <  1:4495134/70‑1 (MQ=255)
               tacatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                            >  1:2920532/1‑56 (MQ=255)
               tacatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                            >  1:1629394/1‑56 (MQ=255)
               tacatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                             <  1:3734994/71‑1 (MQ=255)
                acatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                            >  1:2035174/1‑71 (MQ=255)
                acatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                            >  1:2040400/1‑71 (MQ=255)
                acatTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                            >  1:1209206/1‑71 (MQ=255)
                  atTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCg                                          >  1:688489/1‑71 (MQ=255)
                  atTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGCTCCAGAAAACAATCATTCg                                          >  1:2389113/1‑71 (MQ=255)
                     aGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCgg                                         >  1:211224/1‑69 (MQ=255)
                     aGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGa                                        >  1:2464/1‑70 (MQ=255)
                     aGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGa                                        <  1:2660987/70‑1 (MQ=255)
                         tACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTcc                                                          >  1:3549348/1‑48 (MQ=255)
                          aCAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACaa                                                 <  1:1965428/56‑1 (MQ=255)
                          aCAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                                    <  1:1994461/69‑1 (MQ=255)
                          aCAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                                   >  1:70335/1‑70 (MQ=255)
                          aCAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                                   <  1:2283944/70‑1 (MQ=255)
                              tCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACa                                                  <  1:3408991/51‑1 (MQ=255)
                               cGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                               >  1:3021272/1‑69 (MQ=255)
                                 cTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                                <  1:4014380/50‑1 (MQ=255)
                                    tAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                          <  1:257536/69‑1 (MQ=255)
                                            aCAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaa                                      >  1:434156/1‑49 (MQ=255)
                                              aGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                             <  1:1622517/40‑1 (MQ=255)
                                              aGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc              <  1:3316630/71‑1 (MQ=255)
                                               gCATAACCTTTACATAATTTAGTTCCAGAAAACaa                                                 <  1:2365986/35‑1 (MQ=255)
                                                 aTAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCg                                          >  1:1316150/1‑40 (MQ=255)
                                                 aTAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                                    >  1:960543/1‑46 (MQ=255)
                                                  tAACCTTTACATAATTTAGTTCCAGAAAACAATCACTCGGaaaaa                                    >  1:3219450/1‑45 (MQ=255)
                                                    aCCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                               >  1:3606133/1‑48 (MQ=255)
                                                    aCCTTTACATAATTTAGCTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTa                  >  1:3679668/1‑61 (MQ=255)
                                                     ccTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc              <  1:393305/64‑1 (MQ=255)
                                                       tttACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCa            <  1:784663/64‑1 (MQ=255)
                                                       tttACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCa            <  1:544254/64‑1 (MQ=255)
                                                       tttACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGt      >  1:1578016/1‑70 (MQ=255)
                                                         tACATAATTTAGCTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGtt     <  1:1103721/69‑1 (MQ=255)
                                                          aCATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                               >  1:1233738/1‑42 (MQ=255)
                                                          aCATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTcc             >  1:3349038/1‑60 (MQ=255)
                                                          aCATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCa            >  1:1092428/1‑61 (MQ=255)
                                                          aCATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGtt     <  1:396207/68‑1 (MQ=255)
                                                          aCATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTa    <  1:27880/69‑1 (MQ=255)
                                                            aTAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAac                       >  1:887843/1‑48 (MQ=255)
                                                            aTAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGt                   <  1:1448972/52‑1 (MQ=255)
                                                            aTAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTAtt  <  1:4267523/69‑1 (MQ=255)
                                                             tAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa                        >  1:4596433/1‑46 (MQ=255)
                                                             tAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc              <  1:2164122/56‑1 (MQ=255)
                                                                     |                                                           
ATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCATGGGGTTATT  >  AP009048/1872005‑1872133

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: