Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	164,557	G→C	5.6%	intergenic (+23/‑173)	hrpB → / → mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	164,557	0	G	C	5.6%	273.7 / 6.2	81	intergenic (+23/‑173)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base G (34/42);  new base C (2/3);  total (36/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-02

TGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTAC  >  NC_000913/164284‑164831                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                      tGGCTGTGAGAATGCAAGGGATGTTTGGCGAGGCCACCAATCCGACGATTGTCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtgt                                                                                                                                                                                                                                                                                      <  2:252590/275‑1 (MQ=255)       gtggaATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgaga                                                                                                                                                                                                                                                                                <  2:338/272‑1 (MQ=255)                                            gacggTCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                                                                                                                           <  2:266468/275‑1 (MQ=255)                                            gacgaGCGCCCATGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                                                                                                                           <  2:68635/275‑1 (MQ=255)                                             acgaTCGCCCGGGGGCGCGTGCCGCTGGTGCTGGTGTGGCTTTCACCTGCCCAAAGGCCATGACAAATCACACAGGATTGGAGCGACTTCTAGAAATGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                                                                                                                           <  2:254905/274‑1 (MQ=255)                                              cgaTCGGCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAAGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                                                                                                                           <  2:168194/273‑1 (MQ=255)                                                         gggCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGACCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                                                                                                              <  2:67149/275‑1 (MQ=255)                                                                         tGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCTCACCCTAACCCTCTCCCTCAAggggaag                                                                                                                                                                                                               >  2:442918/1‑270 (MQ=255)                                                                            tGGAGTTGCTTTAACCGTCCCAAAGGCCATTACAAATCACACGAGATTGACGCGACGGCTGGAAAGGGGCGTACCGTGGGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGGCCCGGCAAATACTGCACCGACGCGACGGAAGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGTGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggt                                                                                                                                                                                                            <  4:446251/274‑2 (MQ=255)                                                                                  tGCTTCCACCGTCCCAAAGGCCATTACAGATCACACCAGTGTTGAGCGACTTCTGGAAAGTGGCGTACCGTGAGGTGCAAAAAGGGATGAAAGGGCGTTATCACAAACACGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCGTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTTCACCCTCACCCTTTCCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                                                                                                                      <  4:259678/274‑1 (MQ=255)                                                                                       tCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGGGCGTAACGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACAAGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtg                                                                                                                                                                                                                                                                                       <  3:78218/188‑1 (MQ=255)                                                                                                               aTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATAAGCGCGCCCGGTCACACTCACCCTCACCCACTCCCTCAAGTGAGAGGGGACCGATCGGG‑ACAAAATGTTGAGA‑GATCtcttcctc                                                                                                                                                                      >  4:112615/1‑275 (MQ=255)                                                                                                                                 gACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTgg                                                                                                                                                                                                                                                                                    <  7:117571/149‑1 (MQ=255)                                                                                                                                 gACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTgg                                                                                                                                                                                                                                                                                    >  8:117571/1‑149 (MQ=255)                                                                                                                                   cTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGga                                                                                                                                                                                                                                                                                  <  5:294856/149‑1 (MQ=255)                                                                                                                                    ttCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGgag                                                                                                                                                                                                                                                                                 >  6:475593/1‑149 (MQ=255)                                                                                                                                        ggAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                                                                                                                                                             <  5:66304/149‑1 (MQ=255)                                                                                                                                         gAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggga                                                                                                                                                                                                                                                                            <  6:357498/149‑2 (MQ=255)                                                                                                                                             ggAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                                                                                                                                        <  7:409133/149‑1 (MQ=255)                                                                                                                                              gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCACAAA‑TTTGAGA‑GATCTCTTCTCCTGTCTTGTCACCGA‑AGAAcctgtctcttat                                                                                                                                       >  4:133498/1‑263 (MQ=255)                                                                                                                                               aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGt                                                                                                                                                                                                                                                                      <  5:448757/149‑1 (MQ=255)                                                                                                                                               aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGCGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCt                                                                                                                                     >  2:23442/1‑275 (MQ=255)                                                                                                                                                  gTACCGTGAGGTGCAAAAAGGGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGt                                                                                                                                                              <  2:133499/248‑1 (MQ=255)                                                                                                                                                  gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                                                                                                                                   <  5:475593/149‑1 (MQ=255)                                                                                                                                                  gTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAg                                                                                                                                                                                                                                                                   <  8:322965/149‑1 (MQ=255)                                                                                                                                                        tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTg                                                                                                                                                                                      >  2:513762/1‑219 (MQ=255)                                                                                                                                                               cAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAgcgcgc                                                                                                                                                                                                                                                      <  7:187911/149‑1 (MQ=255)                                                                                                                                                                 aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                                                                                                                                                             >  7:103809/1‑124 (MQ=255)                                                                                                                                                                 aaaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCACAGGGAGAggg                                                                                                                                                                                                                                                                             <  8:103809/124‑1 (MQ=255)                                                                                                                                                                  aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAC‑CTCAAATATTGAGA‑GATATCTTCTTCTGTCTGGTAACAGA‑CTAACAAAAAATCCGGCTTTTGCGCCTGAATATTGCg                                                                                                                  >  4:398795/1‑275 (MQ=255)                                                                                                                                                                   aaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCtg                                                                                                                                                                                                                                                                                       <  2:328806/112‑1 (MQ=255)                                                                                                                                                                       gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                                                                                                             >  2:209396/1‑275 (MQ=255)                                                                                                                                                                        aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                                                                                                                             <  8:151513/149‑1 (MQ=255)                                                                                                                                                                        aTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                                                                                                                             >  7:151513/1‑149 (MQ=255)                                                                                                                                                                         tGANNGGGCGTTATCCCAAACATGGCTGGCCGGACGACCCGGCAAATACTGCANCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTgaga‑g                                                                                                                                                                                 <  4:513757/206‑1 (MQ=255)                                                                                                                                                                         tGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccc                                                                                                                                                                                                                                            >  5:419271/1‑149 (MQ=255)                                                                                                                                                                          gAAAGGGCGTTTTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccct                                                                                                                                                                                                                                           >  5:353147/1‑149 (MQ=255)                                                                                                                                                                             agtgCGTTAGCCCAAGCAGGTCTGGCCGGACGGCCAGGAAAATAAGGGACCGAGGCGAAGGATGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCGGGGGGGGAGGGCCGGCGTGAGGGCATCAGCGCGCACGTTCACTCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAG‑GAAAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTACCAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAGAAGc                                                                                                        <  4:23442/271‑1 (MQ=255)                                                                                                                                                                              gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                                                                                                                                        >  7:8465/1‑116 (MQ=255)                                                                                                                                                                              gggCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCggg                                                                                                                                                                                                                                                                        <  8:8465/116‑1 (MQ=255)                                                                                                                                                                                  gTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                                                                                                                                   >  5:364/1‑149 (MQ=255)                                                                                                                                                                                     aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGtt                                                                                                                                                                                                                                                 >  4:518350/1‑132 (MQ=255)                                                                                                                                                                                     aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCAACCTAACCCCCTCCCTCAAGGGAGAGGGGACCGAACGAG‑CCAAAAGTTTGAGA‑GATATCTTCTTCTGTCTTGTTAAAGA‑AGAACAGACAACCGGGCTTTTGCGCCTGAATATTGAGGAAAAAAAACATGGCCcgg                                                                                               >  4:85601/1‑275 (MQ=255)                                                                                                                                                                                         cAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcacc                                                                                                                                                                                                                                             <  2:518355/132‑1 (MQ=255)                                                                                                                                                                                         cAAACATGTCGTGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACCAAAAAGTATTCGCAAGGCTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCtctc                                                                                                                                                                                                                            <  7:345930/149‑1 (MQ=38)                                                                                                                                                                                              aTGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACcg                                                                                       <  2:112616/274‑1 (MQ=255)                                                                                                                                                                                               tGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCaa                                                                                                                                                                                                                      <  8:315993/149‑1 (MQ=255)                                                                                                                                                                                                 tCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGCCGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCACCTCCCTGTGGGCGATGGCCTGGGTGAGGGCATCAGCGCGCCAGCTCATCCTCACCCTAACCCTCCCCCTAAAGGGAGAGGGGGCCGACCGAG‑CACAAA‑CTTGATACGCTCTCTTCTTCTGTCGTTTACCAGACACAAC‑GAAAATAGTGTTTGTGAGCCTGCATATTGCGGGGAAAAAGCAGGACCGGTAATGGccaggac                                                                                   >  4:193408/1‑270 (MQ=255)                                                                                                                                                                                                   tGGCCGGACTACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTAGGAAGGGTGGTTCTCTCCCTCTCCCTCTGGTAGAGAGCCGGGGAGAGGGCAGCAGCGCGAACGTTAACCCTTACCCTAACGCTACCTCTCAAGGGAGAGGGGACCCATGCAG‑CACAACATTTGAGA‑GTTATCTTCTTCTGTCCTGTAACAGA‑CGAACAGACAATCCGCGTTTTCCGCTCGAATCGGGCGGTAAAAAAGCACGCCCGGGACGTACCCCGAGcc                                                                                 >  4:302312/1‑275 (MQ=255)                                                                                                                                                                                                       cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                                                                                              <  6:364/149‑1 (MQ=255)                                                                                                                                                                                                       cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                                                                                               <  8:150318/148‑1 (MQ=255)                                                                                                                                                                                                       cGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGagag                                                                                                                                                                                                               >  7:150318/1‑148 (MQ=255)                                                                                                                                                                                                         gacgacCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAAc                                                                                                                                                                                                                                  >  4:30934/1‑127 (MQ=255)                                                                                                                                                                                                             acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCctc                                                                                                                                                                                                                              <  2:30934/127‑1 (MQ=255)                                                                                                                                                                                                             acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc                                                                       <  2:85602/275‑1 (MQ=255)                                                                                                                                                                                                                             aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                                                                                                                         <  8:56871/116‑1 (MQ=255)                                                                                                                                                                                                                             aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCt                                                                                                                                                                                                                         >  7:56871/1‑116 (MQ=255)                                                                                                                                                                                                                             aCCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAAtttt                                                                                                                                                                                       >  7:191583/1‑149 (MQ=255)                                                                                                                                                                                                                                 acgcgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCACGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTcc                                                   >  2:81427/1‑275 (MQ=255)                                                                                                                                                                                                                                    cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGC‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATTTTGCGGAGAAAAAGCATTGCCGGGAATGACCGCGAGACAAGTGGACGCAAAGGGAAACCGACCCGTCCggg                                                >  4:215484/1‑274 (MQ=255)                                                                                                                                                                                                                                    cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑Gat                                                                                                                                                                               <  6:419271/149‑1 (MQ=255)                                                                                                                                                                                                                                    cgacgGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCGGTGGGGGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGGGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGt                                                <  2:302316/275‑1 (MQ=255)                                                                                                                                                                                                                                     gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                                                                                              <  8:307252/119‑1 (MQ=255)                                                                                                                                                                                                                                     gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGAGAgg                                                                                                                                                                                                              >  7:307252/1‑119 (MQ=255)                                                                                                                                                                                                                                     gacggacgAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTc                                               <  2:416390/275‑1 (MQ=255)                                                                                                                                                                                                                                               aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                                                                                                                                   <  6:95062/88‑1 (MQ=255)                                                                                                                                                                                                                                               aaaGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTaa                                                                                                                                                                                                                                   >  5:95062/1‑88 (MQ=255)                                                                                                                                                                                                                                                   tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCtt                                                                                                                                                                >  7:467874/1‑149 (MQ=255)                                                                                                                                                                                                                                                    aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg                                                                           >  2:474755/1‑232 (MQ=255)                                                                                                                                                                                                                                                        gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                                                                                               >  8:375870/1‑99 (MQ=255)                                                                                                                                                                                                                                                        gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                                                                                               <  7:375870/99‑1 (MQ=255)                                                                                                                                                                                                                                                        gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAAc                                                                                                                                                           >  5:140671/1‑149 (MQ=255)                                                                                                                                                                                                                                                         tAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCcgtc                            <  2:193410/274‑1 (MQ=255)                                                                                                                                                                                                                                                           aGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACaga                                                                                                                                                        <  6:353147/149‑1 (MQ=255)                                                                                                                                                                                                                                                              ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTga                                                                                                                                                                                     >  6:47830/1‑118 (MQ=255)                                                                                                                                                                                                                                                              ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTga                                                                                                                                                                                     <  5:47830/118‑1 (MQ=255)                                                                                                                                                                                                                                                              ttggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACaga‑aga                                                                                                                                                    >  5:270941/1‑149 (MQ=255)                                                                                                                                                                                                                                                                ggttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAAc                                                                                                                                                  <  8:191583/149‑1 (MQ=255)                                                                                                                                                                                                                                                                 gttTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctc                                                                                                                                                                                                                                          >  2:537375/1‑63 (MQ=255)                                                                                                                                                                                                                                                                     tctccttcccctCGTGGGAGAGGGGCGGCGGGAGGTCGTCCGCGCGCGCGTTCACCCGAACCCTAACCCTATCCCTCAAGGGGGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GGTATCTTCTTCTGTCTTGTAAAAGA‑AGAACAGAAAATCGGGCTTGTGCGCCTGAATAGTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGTAAAGGGAAACCGACGCGTCCGGTCAATCAAAAGGGAAGCCGTCGTCGTTACGAAGa               <  4:156042/275‑1 (MQ=255)                                                                                                                                                                                                                                                                     tctccctctcccTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTtcaccctcaccc                                                                                                                                                                                                                                      <  4:537370/63‑1 (MQ=255)                                                                                                                                                                                                                                                                            ctccctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAAc                                                                                                                                                           >  2:523507/1‑129 (MQ=255)                                                                                                                                                                                                                                                                                ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACaga‑a                                                                                                                                                      <  4:523502/129‑1 (MQ=255)                                                                                                                                                                                                                                                                                ctGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCtttt                                                                                                                                  >  5:202333/1‑149 (MQ=255)                                                                                                                                                                                                                                                                                 tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAAt                                                                                                                                                                                          >  6:390959/1‑94 (MQ=255)                                                                                                                                                                                                                                                                                 tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAAt                                                                                                                                                                                          <  5:390959/94‑1 (MQ=255)                                                                                                                                                                                                                                                                                  gtgGGGGAGGGACGGGGCTAGGGGTTCAGCGCGCACGTCCAACGTCCACCTAACCCATGCCCTAAATGGAGAGGGGAACGATCGAG‑AACAAATTTTGAGA‑GTTATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGTTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac  <  4:81426/275‑1 (MQ=255)                                                                                                                                                                                                                                                                                  gtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTgc                                                                                                                                <  6:270941/149‑1 (MQ=255)                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                      TGGCGGTGAGAATGCAAGAGATGTTTGGCGAGGCCACCAATCCGACGATCGCCCAGGGGCGCGTGCCGCTGGTGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTAC  >  NC_000913/164284‑164831

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