Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	402,391	Δ723 bp	coding (645‑1367/1416 nt)	phoA →	alkaline phosphatase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	402391	403113	723	34 [0]	[0] 34	phoA	alkaline phosphatase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 402390	0 (0.000)	34 (0.620)	31/482	0.7	100%	coding (644/1416 nt)	phoA	alkaline phosphatase
	?	NC_000913	403114 =	0 (0.000)					coding (1368/1416 nt)	phoA	alkaline phosphatase

TGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAC  >  NC_000913/402870‑403354                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 gtgtcaaaacctataacggcgcgctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGAcc                                                                                                                                                                                                                                            <  2:152121‑M2/7‑1 (MQ=255)     caaatcctataacggcgcgctggtcgtcgatattcacgaaaaagatctcccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttcttccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgtccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGAtc                                                                                                                                                                                                                                        <  1:82920‑M2/11‑1 (MQ=255)             ataacggcgcgctgggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTac                                                                                                                                                                                                                                >  1:191841‑M2/233‑251 (MQ=255)                          gggcgtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGc                                                                                                                                                                                                                     <  1:375258‑M2/30‑1 (MQ=255)                              gtcgatattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGctct                                                                                                                                                                                                               >  1:179416‑M2/216‑251 (MQ=255)                                  atattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTgggg                                                                                                                                                                                                           >  1:12940‑M2/212‑251 (MQ=255)                                  atattcacgaaaaagatcacccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTgggg                                                                                                                                                                                                           >  1:519468‑M2/212‑251 (MQ=255)                                       cacgaaaaagatcccccaacgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGa                                                                                                                                                                                                       <  2:376745‑M2/44‑1 (MQ=255)                                            aaaagatcacccaacgattctggaaatggcacaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGCCCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATaa                                                                                                                                                                                                  >  2:338477‑M2/202‑250 (MQ=255)                                                          cgattctggaaatggcaaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAg                                                                                                                                                                                   <  2:314315‑M2/64‑1 (MQ=255)                                                                          aaaagccgcaggtctggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCg                                                                                                                                                                   >  1:486081‑M2/172‑251 (MQ=255)                                                                                         ggcgaccggtaacgtttctaccgcagagttgcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCGCGTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCGGCCGGGGGGTTTTTTTg                                                                                                                                                      <  2:83442‑M2/93‑1 (MQ=255)                                                                                                                       gcaggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACg                                                                                                                      <  1:555918‑M2/125‑1 (MQ=255)                                                                                                                        caggatgccacgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACgg                                                                                                                     >  1:191576‑M2/126‑251 (MQ=255)                                                                                                                                  cgcccgctgcgctggtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGc                                                                                                           <  1:181545‑M2/136‑1 (MQ=255)                                                                                                                                           cgctggtggcacatgtgaccgcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTa                                                                                                  <  2:114681‑M2/145‑1 (MQ=255)                                                                                                                                                gtggcacatgtgacctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAAc                                                                                              >  1:106872‑M2/102‑250 (MQ=255)                                                                                                                                                            acctcgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCAGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACTTGAAGAGGGATg                                                                                 >  2:491037‑M2/90‑251 (MQ=255)                                                                                                                                                                cgcgcaaatgctacggtccgagcgcgaccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTg                                                                              <  2:15388‑M2/165‑1 (MQ=255)                                                                                                                                                                           tacggtccgagcgcgaccagtgaaaaatggccgggtaacgctctggaaaaaggcagaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCAGGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTTCTATGAAAATa                                                                   <  2:59316‑M2/176‑1 (MQ=255)                                                                                                                                                                                           ccagtgaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCtt                                                  >  1:73417‑M2/59‑251 (MQ=255)                                                                                                                                                                                              gggaaaaatggccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACAGGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCt                                               <  1:335791‑M2/196‑1 (MQ=255)                                                                                                                                                                                                gaaaaatgtccgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCtt                                              <  2:558091‑M2/197‑1 (MQ=255)                                                                                                                                                                                                  aaaatgtccgggtcacgatctggacaaaggcggcaaaggatagattacagaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATAACGGGCGCATAAGCTCATGGTTAAAACATGACGAGGGATGGTGCTAAGAAAATAACATTACTGGTTACCTTGcttttt                                           >  2:389642‑M2/52‑250 (MQ=255)                                                                                                                                                                                                   aaatgtccgtgtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCGCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATGTTCGCTGCCGGGTGGTTTTTTGGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCAGACGCTCATGGTTAAAGCATTAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCg                                          <  1:82704‑M2/201‑1 (MQ=255)                                                                                                                                                                                                         ccgggtaacgctcttgaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTgg                                    <  1:461624‑M2/207‑1 (MQ=255)                                                                                                                                                                                                          cgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgc                                                                                                                   <  1:498341‑M2/128‑1 (MQ=255)                                                                                                                                                                                                          cgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGgc                                                                                                                   >  2:498341‑M2/44‑171 (MQ=255)                                                                                                                                                                                                          cgggtaacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGt                                   >  1:307919‑M2/44‑251 (MQ=255)                                                                                                                                                                                                              taacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttt                               >  1:37565‑M2/40‑251 (MQ=255)                                                                                                                                                                                                              taacgctctggaaaaaggcggaaaaggatcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGttttt                               >  2:19152‑M2/40‑251 (MQ=255)                                                                                                                                                                                                                                        gatcgattaccggCCAGACCGATCTCGTCTACGCCATGAAAGCCGCTCTGGGGCTGAAATAAACCCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGATGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCgagag     <  1:96685‑M2/238‑1 (MQ=255)                                                                                                                                                                                                                                         atcgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGCCTTAATGGCAGATCACGGGCGCATACACTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTATTTTTAACCACCGTCGGCGCTGCCgag       >  1:490881‑M2/13‑248 (MQ=255)                                                                                                                                                                                                                                           cgattaccgaCCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAc  <  1:428995‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 TGGAATTCGCTAAAAAGGAGGGTAACACGCTGGTCATAGTCACCGCTGATCACGCCCACGCCAGCCAGATTGTTGCGCCGGATACCAAAGCTCCGGGCCTCACCCAGGCGCTAAATACCAAAGATGGCGCAGTGATGGTGATGAGTTACGGGAACTCCGAAGAGGATTCACAAGAACATACCGGCAGTCAGTTGCGTATTGCGGCGTATGGCCCGCATGCCGCCAATGTTGTTGGACTGACCGACCAGACCGATCTCTTCTACACCATGAAAGCCGCTCTGGGGCTGAAATAAAACCGCGCCCGGCAGTGAATTTTCGCTGCCGGGTGGTTTTTTTGCTGTTAGCAACCAGACTTAATGGCAGATCACGGGCGCATACGCTCATGGTTAAAACATGAAGAGGGATGGTGCTATGAAAATAACATTACTGGTTACCTTGCTTTTCGGTCTGGTTTTTTTAACCACCGTCGGCGCTGCCGAGAGAAC  >  NC_000913/402870‑403354

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