Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,170,593 G→C 5.3% intergenic (‑184/‑57) ycfQ ← / → ycfR predicted DNA‑binding transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,170,5930GC 5.3% 235.2 / 2.1 75intergenic (‑184/‑57)ycfQ/ycfRpredicted DNA‑binding transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (33/38);  new base C (2/2);  total (35/40)
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.98e-01

TTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAA  >  AP009048/1170525‑1170663
                                                                    |                                                                      
tttttAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGtt                                                                      <  1:3135187/71‑1 (MQ=255)
  tttAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTa                                                                     <  1:4063473/70‑1 (MQ=255)
  tttAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTa                                                                     <  1:387684/70‑1 (MQ=255)
   ttAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTat                                                                    <  1:2194379/70‑1 (MQ=255)
    tAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatc                                                                   <  1:2488095/70‑1 (MQ=255)
     aaaTTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcg                                                                  <  1:3149691/70‑1 (MQ=255)
      aattaattAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcga                                                                 <  1:5676155/70‑1 (MQ=255)
       attaattaATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTat                                                                    <  1:7348689/66‑1 (MQ=255)
         taattaatGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTatcg                                                                  >  1:1229146/1‑66 (MQ=255)
              aatGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAg                                                         >  1:4881778/1‑70 (MQ=255)
               atGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAg                                                         <  1:4101531/69‑1 (MQ=255)
               atGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAg                                                         <  1:3237445/69‑1 (MQ=255)
               atGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTACTTaacga                                                                 <  1:5351092/61‑1 (MQ=255)
                    tGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGt                                                                       <  1:919782/50‑1 (MQ=255)
                     gTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgc                                                  <  1:5931813/70‑1 (MQ=255)
                     gTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgc                                                  <  1:5941224/70‑1 (MQ=255)
                     gTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCAATTACTTAACGATCGATAACTAAatga                                                  <  1:1035511/70‑2 (MQ=21)
                      ttATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGt                                                                       <  1:6579769/48‑1 (MQ=255)
                      ttATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTat                                                                    <  1:3052339/51‑1 (MQ=255)
                        ataAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTa                                                                     >  1:6579614/1‑48 (MQ=255)
                        ataAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTaa                                                      >  1:4281145/1‑63 (MQ=255)
                        ataAAAAATATCTTGTATGTGATCAAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttgc                                              >  1:5248320/1‑71 (MQ=255)
                         taAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttgc                                              <  1:5560951/70‑1 (MQ=255)
                         taAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAtt                                                    >  1:907005/1‑64 (MQ=255)
                         taAAAAATATCTTGTATGTGATACAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAtt                                                    >  1:3414634/1‑64 (MQ=255)
                           aaaaaTATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttgc                                              <  1:2201447/68‑1 (MQ=255)
                           aaaaaTATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttgc                                              <  1:3767466/68‑1 (MQ=255)
                             aaaTATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACg                                          >  1:5909634/1‑70 (MQ=255)
                             aaaTATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGATTGCGACg                                          >  1:5667662/1‑70 (MQ=255)
                               atatCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGt                                         <  1:3123376/69‑1 (MQ=255)
                               atatCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGt                                         <  1:4932635/69‑1 (MQ=255)
                                  tCTTGTATGTGATCCAGATCACATCTATCATTTACTTaacgatc                                                               >  1:4859465/1‑44 (MQ=38)
                                  tCTTGTATGTGATACAGATCACATCTATCATTTACTTaacgatc                                                               >  1:1739649/1‑44 (MQ=38)
                                      gTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtca                                       <  1:4534993/64‑1 (MQ=255)
                                      gTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtca                                       <  1:4967009/64‑1 (MQ=255)
                                      gTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtca                                       <  1:2169628/64‑1 (MQ=255)
                                       tATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtcattcat                                  >  1:2095393/1‑68 (MQ=255)
                                       tATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtca                                       >  1:5199868/1‑63 (MQ=255)
                                       tATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtca                                       >  1:6674403/1‑63 (MQ=255)
                                       tATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCt                                <  1:3126658/70‑1 (MQ=255)
                                        aTGTGATCTGGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtcattcat                                  <  1:2954131/67‑1 (MQ=255)
                                        aTGTGATCTGGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtcattcat                                  <  1:3305297/67‑1 (MQ=255)
                                        aTGTGATCTGGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGtcattcat                                  <  1:358778/67‑1 (MQ=255)
                                        aTGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttgc                                              <  1:7636301/55‑1 (MQ=255)
                                        aTGTGATCCAGATCACATCTAACATTTAGTTATCGATCGTTAAGTAAttgcttgc                                              <  1:6480123/55‑1 (MQ=255)
                                           tgATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCAt                            >  1:1857077/1‑70 (MQ=255)
                                           tgATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCAt                            <  1:7028264/70‑1 (MQ=255)
                                            gATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCAt                            >  1:1629033/1‑69 (MQ=255)
                                              tCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttg                                               >  1:5106739/1‑48 (MQ=255)
                                              tCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgcttg                                               >  1:4382584/1‑48 (MQ=255)
                                               ccAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGGGACGTCATTCATCTGCATAAGGc                       >  1:7545293/1‑71 (MQ=255)
                                               ccAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGc                       >  1:7038198/1‑71 (MQ=255)
                                                  gATCACATCTATCATTTAGTTATCGATCGTTAAGTAAttgctt                                                <  1:1984099/43‑1 (MQ=255)
                                                    tCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACta                  >  1:6672162/1‑71 (MQ=255)
                                                     cacaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATc                                 >  1:2891139/1‑55 (MQ=255)
                                                     cacaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATc                                 >  1:6158526/1‑55 (MQ=255)
                                                     cacaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATc                                 >  1:4235245/1‑55 (MQ=255)
                                                      acaTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACtat                 >  1:3539652/1‑70 (MQ=255)
                                                        aTCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACtattat              >  1:2222127/1‑71 (MQ=255)
                                                         tCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGc                              <  1:6515661/54‑1 (MQ=255)
                                                         tCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGc                              <  1:486546/54‑1 (MQ=255)
                                                         tCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACtagtat              >  1:299461/1‑70 (MQ=255)
                                                          cTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCAt                            >  1:3587773/1‑55 (MQ=255)
                                                          cTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCAt                            >  1:3868356/1‑55 (MQ=255)
                                                          cTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACtatt                >  1:1522035/1‑67 (MQ=255)
                                                           tATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCAc                    <  1:890851/62‑1 (MQ=255)
                                                           tATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCAc                    <  1:4217361/62‑1 (MQ=255)
                                                           tATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGa            >  1:6654024/1‑70 (MQ=255)
                                                             tCATTTAGTTATCGATCGTTAAGTAATTTCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGaaa          >  1:3915835/1‑70 (MQ=255)
                                                             tCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGaaa          >  1:689804/1‑70 (MQ=255)
                                                             tCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGaaa          <  1:7462366/70‑1 (MQ=255)
                                                                 ttAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACg      >  1:4960508/1‑70 (MQ=255)
                                                                 ttAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACg      >  1:5258517/1‑70 (MQ=255)
                                                                   aGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTa    <  1:2997798/70‑1 (MQ=255)
                                                                   aGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTa    <  1:4230979/70‑1 (MQ=255)
                                                                    gTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACg      <  1:5609991/67‑1 (MQ=255)
                                                                    gTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTaaa  <  1:4600971/71‑1 (MQ=255)
                                                                    gTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTaa   <  1:766437/70‑1 (MQ=255)
                                                                    |                                                                      
TTTTTAAATTAATTAATGATTGTTATAAAAAATATCTTGTATGTGATCCAGATCACATCTATCATTTAGTTATCGATCGTTAAGTAATTGCTTGCGACGTCATTCATCTGCATAAGGCCACTATTATGAAAAACGTAAA  >  AP009048/1170525‑1170663

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: