Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,269,672 (A)7→6 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,6660A.97.4% 148.4 / ‑2.7 38coding (253/714 nt)aphAacid phosphatase/phosphotransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base . (17/20);  minor base C (1/0);  total (18/20)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTACGCCGCG  >  NC_000913/4269534‑4269801
                                                                                                                                    |                                                                                                                                       
gCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGGTGACACGGGACTATATTACAGTCCGTGCTTAAGGCGCGGCaaaaaacc                                                                                                                                  >  2:407506/1‑138 (MQ=255)
      cAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTGCCAGTCCTGGCGTCTGGCGCGCCAAAAAACCTTCTCg                                                                                                                           >  1:388041/1‑140 (MQ=255)
              cattcattGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAg                                                                                                                   <  2:798661/140‑1 (MQ=255)
                ttcattGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCCGACTTCTGGGGCGGCACAAAACCGTCTCGCCAGAAAGCg                                                                                                                 >  1:351824/1‑140 (MQ=255)
                   attGGGTTTCGGTCGCCAAAATTGAAAATAGGCACGCAGAGCGTCCGCGAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAg                                                                                                              <  2:388041/140‑1 (MQ=255)
                      gggTTTCGGTCGCACAAGTTGACAATAGCCTCGCACGGCATCCGCCAACGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGAtt                                                                                                           <  2:276968/140‑1 (MQ=255)
                       ggTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGGTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTa                                                                                                          >  2:160990/1‑140 (MQ=255)
                        gTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTAt                                                                                                         <  2:17420/140‑1 (MQ=255)
                         tttCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATc                                                                                                        >  1:580824/1‑140 (MQ=255)
                           tCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTg                                                                                                      >  1:192035/1‑140 (MQ=255)
                             ggTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                                    <  2:1190/140‑1 (MQ=255)
                                  cacaAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                                    >  2:514554/1‑135 (MQ=255)
                                  cacaAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGaa                                                                                                    <  1:514554/135‑1 (MQ=255)
                                              aTAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTggg                                                                                   <  1:160990/140‑1 (MQ=255)
                                               tAGCCTCGCAGGGCGTCTGCCACTGGCGGTGCGGTGTGATATCGATGATACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGa                                                                                  <  2:351824/140‑1 (MQ=255)
                                                aGCCTCGCAGAGCGTACGCCAAAGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGaa                                                                                 <  1:350493/140‑1 (MQ=255)
                                                  ccTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAATATTATCTGGAAAAGCCTGTGTTCTGGGaaaa                                                                               >  2:306756/1‑140 (MQ=255)
                                                           gCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACaa                                                                      <  2:192035/140‑1 (MQ=255)
                                                                cGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCt                                                                 >  1:131002/1‑140 (MQ=255)
                                                                   cAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTggg                                                              <  2:472172/140‑1 (MQ=255)
                                                                   cAATGGCGGTGGGGTATGATATCGATGACGCGGTACTGTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTggg                                                              <  1:407506/140‑1 (MQ=255)
                                                                             ggggTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAg                                                    <  2:131002/140‑1 (MQ=255)
                                                                                  ttGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTc                                               >  2:601341/1‑140 (MQ=255)
                                                                                        tCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAag                                         >  2:442338/1‑140 (MQ=255)
                                                                                         cGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAAtt                                                       >  2:855878/1‑125 (MQ=255)
                                                                                         cGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAAtt                                                       <  1:855878/125‑1 (MQ=255)
                                                                                            tGACACGGTACTTTTTTCCAGTCCGCGCTCCTGGCGCTGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGt                                     <  1:331064/140‑1 (MQ=255)
                                                                                              acacGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGtcg                                   <  1:306756/140‑1 (MQ=255)
                                                                                                    tACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGcc                             >  1:160740/1‑140 (MQ=255)
                                                                                                               aGTCCGGGCTTCTGGCGCGGCACAAAACCTTCTCGCCAGTAAGCGAAGCTTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCtgattgat                  <  2:231613/140‑1 (MQ=255)
                                                                                                                gTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAACAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGata                 >  1:810390/1‑140 (MQ=255)
                                                                                                                gTCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTGAGCATTCGAAAGGAGGTCGCTAGCCAGCTGATTGata                 >  1:231613/1‑140 (MQ=255)
                                                                                                                 tCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGata                 >  1:659512/1‑139 (MQ=255)
                                                                                                                 tCCGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGata                 <  2:659512/139‑1 (MQ=255)
                                                                                                                  ccGGGCTTCTGGCGCGGCAAAACACTTTCGCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATg               <  1:326449/140‑1 (MQ=255)
                                                                                                                   cGGGCTTCTGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGTATGAATTCAGCAGTGCTAAAGAGGTAGCTCAACAGCTGATTGATATGc              >  2:391604/1‑140 (MQ=255)
                                                                                                                           tGGCGCGGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTGAGCATTCCTAAAAAGGTCGCTCGCCCGCTGATTGATATGCATGTAcct      >  2:326449/1‑138 (MQ=255)
                                                                                                                              cgcgGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTAcgccgc   >  2:652857/1‑140 (MQ=255)
                                                                                                                               gcgGCAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTACGCcgcg  <  2:580824/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                       
GCTGAACAGGCACCCATTCATTGGGTTTCGGTCGCACAAATTGAAAATAGCCTCGCAGGGCGTCCGCCAATGGCGGTGGGGTTTGATATCGATGACACGGTACTTTTTTCCAGTCCGGGCTTCTGGCGCGGCAAAAAAACCTTCTCGCCAGAAAGCGAAGATTATCTGAAAAATCCTGTGTTCTGGGAAAAAATGAACAATGGCTGGGATGAATTCAGCATTCCAAAAGAGGTCGCTCGCCAGCTGATTGATATGCATGTACGCCGCG  >  NC_000913/4269534‑4269801

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: