Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,269,672 (A)7→8 coding (259/714 nt) aphA → acid phosphatase/phosphotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,269,6651.A95.7% 149.4 / ‑0.3 46coding (252/714 nt)aphAacid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (20/24);  total (21/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

CACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGC‑‑AAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTA  >  NC_000913/4269544‑4269794
                                                                                                                          |                                                                                                                                  
cACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGcc                                                                                                                  >  1:313231/1‑140 (MQ=255)
  cccATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAAACCTTCTCGCCa                                                                                                                 <  2:1003414/140‑1 (MQ=255)
   ccATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGa                                                                                                               <  2:925895/140‑1 (MQ=255)
    cattcattGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGaa                                                                                                              >  1:112877/1‑140 (MQ=255)
      ttcattGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAg                                                                                                            <  1:612324/140‑1 (MQ=255)
                ttCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTAt                                                                                                  <  1:891961/140‑1 (MQ=255)
                ttCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTAt                                                                                                  >  2:582936/1‑140 (MQ=255)
                 tCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATc                                                                                                 >  2:11812/1‑140 (MQ=255)
                     tCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                             <  1:11812/140‑1 (MQ=255)
                     tCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                             <  2:677678/140‑1 (MQ=255)
                       gCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                             <  2:768361/138‑1 (MQ=255)
                       gCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                             >  1:768361/1‑138 (MQ=255)
                       gCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaaaa                                                                                           >  1:747896/1‑140 (MQ=255)
                        cacaAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaaaaa                                                                                          >  2:1017771/1‑140 (MQ=255)
                        cacaAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaaaaa                                                                                          >  1:311136/1‑140 (MQ=255)
                               tGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCtgtg                                                                                   <  2:311136/140‑1 (MQ=255)
                                        ccTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGaa                                                                          >  1:456525/1‑140 (MQ=255)
                                         cTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGaaa                                                                         >  1:717591/1‑140 (MQ=255)
                                                  cGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACa                                                                <  1:695632/140‑1 (MQ=255)
                                                           aTGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTggg                                                       <  2:112877/140‑1 (MQ=255)
                                                             ggCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGAt                                                     <  1:508431/140‑1 (MQ=255)
                                                                  tGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCt                                                          <  1:102130/130‑1 (MQ=255)
                                                                  tGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAAtt                                                <  2:313231/140‑1 (MQ=255)
                                                                  tGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCt                                                          >  2:102130/1‑130 (MQ=255)
                                                                     gggttGATATCGATGACACGGTACTTTTTTTCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAg                                             >  2:95588/4‑140 (MQ=255)
                                                                               cGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGGAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCaaaa                                   >  1:206531/1‑140 (MQ=255)
                                                                                    acacGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGt                              <  1:592891/140‑1 (MQ=255)
                                                                                      acGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGtcg                            >  2:609433/1‑140 (MQ=255)
                                                                                          tACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGtcgctcg                        >  1:673258/1‑140 (MQ=255)
                                                                                          tACTTTTTTCCAGTCCGGGCTTCTGGCGCGGC‑‑AAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGtcgctcgc                       >  1:547053/1‑140 (MQ=255)
                                                                                           aCTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGtcgctcgc                       >  2:449486/1‑140 (MQ=255)
                                                                                             ttttttCCAGTCCGGGCTTCTGGCGCGGC‑‑AAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAg                    <  2:547053/140‑1 (MQ=255)
                                                                                                  tCCAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtga                <  2:673258/140‑1 (MQ=255)
                                                                                                   ccAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgat               <  1:563323/140‑1 (MQ=255)
                                                                                                    cAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgatt              <  2:747896/140‑1 (MQ=255)
                                                                                                    cAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgatt              <  2:717591/140‑1 (MQ=255)
                                                                                                    cAGTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgatt              >  2:535377/1‑140 (MQ=255)
                                                                                                      gTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgattga            <  2:456525/140‑1 (MQ=255)
                                                                                                      gTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgattga            <  1:609433/140‑1 (MQ=255)
                                                                                                      gTCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgattga            <  1:582936/140‑1 (MQ=255)
                                                                                                       tCCGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgattgat           >  2:764524/1‑140 (MQ=255)
                                                                                                        ccGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGata          <  1:764524/140‑1 (MQ=255)
                                                                                                        ccGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGata          <  2:1074855/140‑1 (MQ=255)
                                                                                                         cGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgat               >  2:366397/1‑134 (MQ=255)
                                                                                                         cGGGCTTCTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgat               <  1:366397/134‑1 (MQ=255)
                                                                                                                cTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATTTGAAAAATCCTGTGTTCTTGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTa  <  2:603413/140‑1 (MQ=255)
                                                                                                                cTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTa  <  2:467970/140‑1 (MQ=255)
                                                                                                                cTGGCGCGGCAAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTa  >  1:651948/1‑140 (MQ=255)
                                                                                                                          |                                                                                                                                  
CACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGC‑‑AAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTA  >  NC_000913/4269544‑4269794

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