Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,634,421 T→A 5.7% intergenic (‑422/‑332) ydfJ ← / → ydfK predicted transporter/predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,634,4210TA 5.7% 195.3 / 3.6 70intergenic (‑422/‑332)ydfJ/ydfKpredicted transporter/predicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (37/29);  new base A (2/2);  total (39/31)
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.85e-01

CCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGC  >  AP009048/1634352‑1634486
                                                                     |                                                                 
ccTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGtc                                                                  <  1:968438/71‑1 (MQ=34)
 cTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGtct                                                                 >  1:6165124/1‑71 (MQ=34)
   tGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGtctc                                                                >  1:1156863/1‑70 (MQ=35)
    gATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTa                                                              >  1:2136698/1‑71 (MQ=35)
        tttttctttttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGACTCTACac                                                           <  1:2075235/70‑2 (MQ=255)
              tttttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATGAGATAg                                                    >  1:590980/1‑71 (MQ=255)
              tttttcAGATAAAAACTGGTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAg                                                    >  1:7159736/1‑71 (MQ=255)
                tttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTctct                                                               >  1:5955414/1‑58 (MQ=31)
                tttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGtt                                                  >  1:1269330/1‑71 (MQ=35)
                tttcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGtt                                                  >  1:1046694/1‑71 (MQ=35)
                  tcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTg                                                >  1:367526/1‑71 (MQ=34)
                  tcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTg                                                >  1:3512238/1‑71 (MQ=34)
                  tcAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTg                                                >  1:2421049/1‑71 (MQ=34)
                   cAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTg                                                <  1:5917985/70‑1 (MQ=35)
                     gATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             <  1:3054039/71‑1 (MQ=34)
                     gATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             <  1:4363383/71‑1 (MQ=34)
                     gATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             <  1:3063060/71‑1 (MQ=34)
                     gATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             <  1:5314806/71‑1 (MQ=34)
                       tAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGAt                                                      <  1:2995601/60‑1 (MQ=35)
                       tAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGtt                                           >  1:5016317/1‑71 (MQ=34)
                       tAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGtt                                           >  1:1526163/1‑71 (MQ=34)
                            aCTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGc                                      >  1:6550882/1‑71 (MQ=38)
                            aCTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGc                                      >  1:7580588/1‑71 (MQ=38)
                            aCTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGc                                      <  1:7567486/71‑1 (MQ=38)
                                 tATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCa                                              >  1:5828848/1‑58 (MQ=35)
                                       tGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTg                             <  1:4132818/69‑1 (MQ=40)
                                       tGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTg                             <  1:5740498/69‑1 (MQ=40)
                                        gTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGa                                                       <  1:3584516/42‑1 (MQ=38)
                                        gTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             >  1:7710688/1‑52 (MQ=35)
                                        gTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTg                                          >  1:6489718/1‑55 (MQ=35)
                                         tataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCt                                     <  1:3271005/59‑1 (MQ=255)
                                         tataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCt                                     <  1:4804965/59‑1 (MQ=255)
                                         tataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAg                                 <  1:1384925/63‑1 (MQ=40)
                                          ataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCt                               <  1:5236554/64‑1 (MQ=40)
                                          ataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCa                        >  1:1304294/1‑71 (MQ=39)
                                          ataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCa                        >  1:3901986/1‑71 (MQ=39)
                                          ataCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCa                        >  1:1313478/1‑71 (MQ=39)
                                           taCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCtt                                    >  1:7061813/1‑58 (MQ=255)
                                           taCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCtt                                    >  1:256481/1‑58 (MQ=255)
                                           taCTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCt                               <  1:2127867/63‑1 (MQ=40)
                                             cTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATaa                     <  1:4256262/71‑1 (MQ=255)
                                             cTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATa                      <  1:6844599/70‑1 (MQ=255)
                                              ttttAAACCCAATCCGTGTAGAGTCTCTACATAAGATAg                                                    >  1:5949754/1‑39 (MQ=40)
                                              ttttAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGt                                                   >  1:119543/1‑40 (MQ=40)
                                              ttttAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGt                                                   >  1:2252510/1‑40 (MQ=40)
                                              ttttAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCtt                              <  1:6516683/61‑1 (MQ=255)
                                                ttAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCtt                                    <  1:6486269/53‑1 (MQ=255)
                                                ttAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTgcgc                          >  1:6958429/1‑63 (MQ=255)
                                                ttAAACACAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATa                      >  1:6983182/1‑67 (MQ=255)
                                                 tAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGc                                >  1:7587216/1‑56 (MQ=255)
                                                 tAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCg                 >  1:4809917/1‑71 (MQ=255)
                                                    aCCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             >  1:5028117/1‑40 (MQ=40)
                                                    aCCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAg                                             >  1:3672732/1‑40 (MQ=40)
                                                    aCCCAATCAGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCt                               >  1:61142/1‑54 (MQ=255)
                                                       ctaTCCGTTTAGAGTCTCTACATAAGATAGTTTGCTGTTGCCGCTTCa                                  <  1:5692034/46‑1 (MQ=14)
                                                        aaTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGc                                      <  1:3499455/43‑1 (MQ=255)
                                                        aaTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGAtt          >  1:3620760/1‑71 (MQ=255)
                                                        aaTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGAtt          >  1:127820/1‑71 (MQ=255)
                                                         aTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGtt                                           <  1:6653760/37‑1 (MQ=40)
                                                         aTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGtt                                           <  1:997799/37‑1 (MQ=40)
                                                         aTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGtt                                           <  1:4869016/37‑1 (MQ=40)
                                                         aTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGAttt         <  1:3772201/71‑1 (MQ=39)
                                                         aTCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGAttt         <  1:1190062/71‑1 (MQ=39)
                                                           ccGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGc                                         >  1:1872866/1‑37 (MQ=40)
                                                            tatgtAGAGACTCTACACAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCt              >  1:3676618/3‑63 (MQ=14)
                                                            tatgtAGAGACTCTACACAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGa            >  1:2783092/3‑65 (MQ=14)
                                                             gtgtAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgc     >  1:639476/1‑71 (MQ=39)
                                                             gtgtAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgc     >  1:3896120/1‑71 (MQ=39)
                                                             atgtAGAGACTCTACACAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGAttt         <  1:3172607/66‑1 (MQ=14)
                                                              tgtAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATaaa                    <  1:2182785/55‑1 (MQ=255)
                                                              tgtAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATaaa                    <  1:299006/55‑1 (MQ=255)
                                                                 agagTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgctgc  >  1:7307516/1‑70 (MQ=255)
                                                                  gagTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTtgc     <  1:6214554/66‑1 (MQ=40)
                                                                     tctcTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATa                      <  1:540802/46‑1 (MQ=255)
                                                                     |                                                                 
CCTTGATATTTTTCTTTTTCAGATAAAAACTGTTATCTATGTATACTTTTAAACCCAATCCGTGTAGAGTCTCTACATAAGATAGTTTGCAGTTGCCGCTTCAGCTTGCGCCATAAACCGCCTGATTTTTGCTGC  >  AP009048/1634352‑1634486

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: