Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 402,739 G→A L331L (TTG→TTA phoA → alkaline phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913402,7390GA94.4% 42.2 / NA 18L331L (TTG→TTAphoAalkaline phosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (9/8);  minor base C (1/0);  total (10/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGG  >  NC_000913/402572‑402965
                                                                                                                                                                       |                                                                                                                                                                                                                                  
ccccTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCTAATCCTTGTgg                                                                                                                                                 <  1:399259/251‑1 (MQ=255)
                               tGCCAGTGCGCTGGCTAGGACCGCAAGCAACGAACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGAtctc                                                                                                                        >  1:169688/1‑245 (MQ=255)
                                         cTGGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAAAGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCCTTTTCCTGCaa                                                                                                                                                                                                   >  1:505043/1‑160 (MQ=255)
                                               aGGACCGAATGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATTAGAAATGTTTTTTCTTGCAAGTTGAAGGTGCGCCAATCGATGATCCGGATCATGCGGCGAATCCGTGTGGGCAAATTCGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGgc                                                                                                    <  1:380887/249‑1 (MQ=255)
                                                  aCCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCACATCCGCATCGTAATGACCGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGCGAAAGGCTTTTTCCTGTAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTgg                                                                                               >  1:497777/1‑251 (MQ=255)
                                                              gTACCATGGCAATATCGATGAGCACGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGAATTCGCTaaa                                                                                    >  1:173889/1‑250 (MQ=255)
                                                               tACCATGGCAATATCGATAAGCCCGCAGGCACCTGTAAGCAAAATCCGCAACGTCATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGACATTGAAttattg                                                                                                                                                                                                                                 >  1:460300/1‑108 (MQ=255)
                                                                aCCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGGCCGACAAAGCCATTGAATTATTGAGTAAAATTGAGAAAGGCTTTTTCCTGCAAGGTGAAGGt                                                                                                                                                                                          <  1:181362/146‑1 (MQ=255)
                                                                 ccATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGAATTCGCTAAAAagn                                                                                <  1:149387/251‑2 (MQ=255)
                                                                                      cGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTCTACCAACCCTGGCGCAGATGACCGCAAAACCCCTTGAATTCTTGAGTAAAATTGAGAAAGGCTTTTTCCTGCAAGCTGACGGTGCGTCAATCGCTAAACAGGATCATGCTGCGAATCATTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGACTTCGCTAAAAaggag                                                                              >  1:472475/1‑232 (MQ=255)
                                                                                        cAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTCCCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGAAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTc                                                         >  1:361510/1‑251 (MQ=255)
                                                                                                 gTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCgn                                                                                                             <  1:305293/190‑2 (MQ=255)
                                                                                                 gTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGANCATNNNNNGann                                                                                                                                                         >  2:305293/1‑144 (MQ=255)
                                                                                                      ccAACTCCGCAACGTAATGACAGTGTACCAACCCTGGCGCTGATGACCGACAAATACATTGAATTATTTAGCAAAAATGAGAAAGGCTTTTTCCTGCAAGTGGAAGGTGCGTCAATCGATAAACAGTATCATGCATCGAATCCTTGTGGGCAAATTGGCGAGACGTTCGGTCTCGATGAACCCGTACAACGGGCGCGGTAAGTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCa                                                        <  1:322274/238‑1 (MQ=255)
                                                                                                                          cAGTGTACCAACCCTGGCGCAGATGACCGACACAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCACGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGAATTAGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGCttgttg                       >  1:168753/1‑251 (MQ=255)
                                                                                                                                       cTGGCGCAGATGACCGACAAAGCCATTGAATTATTTAGTAAAAATGAGAAAGTCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGGGGGCAAATTGGCGAGACTGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGAATTCGCTAAAAATGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATAc             <  1:343925/248‑1 (MQ=255)
                                                                                                                                               naTGACCGACAAAGCCATTGAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCACGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTGCAACGGGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCgg  >  1:404359/2‑251 (MQ=255)
                                                                                                                                                                  naaTTATTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGa                                                                                                                                                                      >  1:297998/2‑68 (MQ=255)
                                                                                                                                                                  gAATTATTGAGTAAAAATGAGAAAGGCTTTTTCCGGCAAGTTGAAGGTGCGTCAATCGATAAACAGGa                                                                                                                                                                      <  2:297998/68‑1 (MQ=255)
                                                                                                                                                                       |                                                                                                                                                                                                                                  
CCCCTGCTTGGCCTGTTTGCTGACGGCAATATGCCAGTGCGCTGGCTAGGACCGAAAGCAACGTACCATGGCAATATCGATAAGCCCGCAGTCACCTGTACGCCAAATCCGCAACGTAATGACAGTGTACCAACCCTGGCGCAGATGACCGACAAAGCCATTGAATTGTTGAGTAAAAATGAGAAAGGCTTTTTCCTGCAAGTTGAAGGTGCGTCAATCGATAAACAGGATCATGCTGCGAATCCTTGTGGGCAAATTGGCGAGACGGTCGATCTCGATGAAGCCGTACAACGGGCGCTGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGG  >  NC_000913/402572‑402965

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