Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	164,620	C→T	10.4%	intergenic (+86/‑110)	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,620	0	C	T	10.4%	237.7 / 17.3	88	intergenic (+86/‑110)	hrpB/mrcB	putative ATP‑dependent helicase/fused glycosyl transferase and transpeptidase
Reads supporting (aligned to +/- strand):  ref base C (42/35);  new base T (4/5);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 7.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-03

TGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAA  >  NC_000913/164356‑164893                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                     tGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTAACCCTCACCCTAACCCTCTCCCTCAAggggaag                                                                                                                                                                                                                                                                             >  2:442918/1‑270 (MQ=255)    tGGAGTTGCTTTAACCGTCCCAAAGGCCATTACAAATCACACGAGATTGACGCGACGGCTGGAAAGGGGCGTACCGTGGGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCGGGCCGGACGGCCCGGCAAATACTGCACCGACGCGACGGAAGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGGGGGAGAGTGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAgggt                                                                                                                                                                                                                                                                          <  4:446251/274‑2 (MQ=255)          tGCTTCCACCGTCCCAAAGGCCATTACAGATCACACCAGTGTTGAGCGACTTCTGGAAAGTGGCGTACCGTGAGGTGCAAAAAGGGATGAAAGGGCGTTATCACAAACACGTCTGGCCGGACGACCCGGCAAATACGGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCGTCTCCCTGTGGGGGAGGGCCGGGGTGAGGGCAGCAGCGCGCACGTTCACCCTCACCCTTTCCCTCTCCCTCAAGGGAGAGGGGACCGAt                                                                                                                                                                                                                                                                    <  4:259678/274‑1 (MQ=255)                                       aTCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATAAGCGCGCCCGGTCACACTCACCCTCACCCACTCCCTCAAGTGAGAGGGGACCGATCGGG‑ACAAAATGTTGAGA‑GATCtcttcctc                                                                                                                                                                                                                                    >  4:112615/1‑275 (MQ=255)                                                                      gAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAGCCACAAA‑TTTGAGA‑GATCTCTTCTCCTGTCTTGTCACCGA‑AGAAcctgtctcttat                                                                                                                                                                                                     >  4:133498/1‑263 (MQ=255)                                                                       aGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGCGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCt                                                                                                                                                                                                   >  2:23442/1‑275 (MQ=255)                                                                          gTACCGTGAGGTGCAAAAAGGGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGt                                                                                                                                                                                                                            <  2:133499/248‑1 (MQ=255)                                                                                tGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTg                                                                                                                                                                                                                                                    >  2:513762/1‑219 (MQ=255)                                                                                          aaaGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCCCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCCCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAC‑CTCAAATATTGAGA‑GATATCTTCTTCTGTCTGGTAACAGA‑CTAACAAAAAATCCGGCTTTTGCGCCTGAATATTGCg                                                                                                                                                                                >  4:398795/1‑275 (MQ=255)                                                                                               gaTGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCGCCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaa                                                                                                                                                                           >  2:209396/1‑275 (MQ=255)                                                                                                 tGANNGGGCGTTATCCCAAACATGGCTGGCCGGACGACCCGGCAAATACTGCANCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTgaga‑g                                                                                                                                                                                                                                               <  4:513757/206‑1 (MQ=255)                                                                                                     agtgCGTTAGCCCAAGCAGGTCTGGCCGGACGGCCAGGAAAATAAGGGACCGAGGCGAAGGATGAAAAAGTAGTCGTAAGGTTGGTTTTCTCCCTCTCCCGGGGGGGGAGGGCCGGCGTGAGGGCATCAGCGCGCACGTTCACTCTCACCCTAACCCTCGCCCTCAAGGGAGAGGGGACCGATCGAG‑GAAAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTACCAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAGAAGc                                                                                                                                                                      <  4:23442/271‑1 (MQ=255)                                                                                                             aTCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCAACCTAACCCCCTCCCTCAAGGGAGAGGGGACCGAACGAG‑CCAAAAGTTTGAGA‑GATATCTTCTTCTGTCTTGTTAAAGA‑AGAACAGACAACCGGGCTTTTGCGCCTGAATATTGAGGAAAAAAAACATGGCCcgg                                                                                                                                                             >  4:85601/1‑275 (MQ=255)                                                                                                                      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)                                                                                                                                     acCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGc                                                                                                                                     <  2:85602/275‑1 (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)                                                                                                                                                                           tATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCtt                                                                                                                                                                                                                              >  7:467874/1‑149 (MQ=255)                                                                                                                                                                            aTTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTgg                                                                                                                                         >  2:474755/1‑232 (MQ=255)                                                                                                                                                                                gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                                                                                                                                                             <  7:375870/99‑1 (MQ=255)                                                                                                                                                                                gTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTTGAGGGagag                                                                                                                                                                                                                                                                             >  8:375870/1‑99 (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)                                                                                                                                                                                             tctccttcccctCGTGGGAGAGGGGCGGCGGGAGGTCGTCCGCGCGCGCGTTCACCCGAACCCTAACCCTATCCCTCAAGGGGGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GGTATCTTCTTCTGTCTTGTAAAAGA‑AGAACAGAAAATCGGGCTTGTGCGCCTGAATAGTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGTAAAGGGAAACCGACGCGTCCGGTCAATCAAAAGGGAAGCCGTCGTCGTTACGAAGa                                                                             <  4:156042/275‑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                                                                                                                                                                                                                                                        <  5:390959/94‑1 (MQ=255)                                                                                                                                                                                                         tgtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAAt                                                                                                                                                                                                                                                        >  6:390959/1‑94 (MQ=255)                                                                                                                                                                                                          gtgGGGGAGGGACGGGGCTAGGGGTTCAGCGCGCACGTCCAACGTCCACCTAACCCATGCCCTAAATGGAGAGGGGAACGATCGAG‑AACAAATTTTGAGA‑GTTATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGTTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTac                                                                <  4:81426/275‑1 (MQ=255)                                                                                                                                                                                                          gtgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTgc                                                                                                                                                                                              <  6:270941/149‑1 (MQ=255)                                                                                                                                                                                                           tgGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTCCGAAGATGACgataatga                                                                 >  2:219723/1‑271 (MQ=255)                                                                                                                                                                                                            gggAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTacga                                                              <  2:249337/275‑1 (MQ=255)                                                                                                                                                                                                               agagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGCAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCCAGCCAATTGGAGGCAAAGGGAAACCGACGCGTCCGGTCAACCAAAAGGTCAGCCGGCGTCGTTACGAAGATGACGATGAttcagaaga                                                           >  2:61443/1‑268 (MQ=255)                                                                                                                                                                                                                gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTgag                                                                                                                                                                                                                                                  <  6:294698/93‑1 (MQ=255)                                                                                                                                                                                                                gagGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTgag                                                                                                                                                                                                                                                  >  5:294698/1‑93 (MQ=255)                                                                                                                                                                                                                 agGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCAGGTCAAACAAAAGGTAAGCAGTCGTCGTTACGAAGACGACGATACTTacgacagtt                                                         >  2:85314/1‑271 (MQ=255)                                                                                                                                                                                                                    gCCGGGGGGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTatg                                                      <  2:215487/275‑1 (MQ=255)                                                                                                                                                                                                                               ggCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCt                                                                                                                                                                                                   >  4:105758/1‑123 (MQ=255)                                                                                                                                                                                                                                gCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGaaaaa                                                                                                                                                                        <  6:140671/149‑1 (MQ=255)                                                                                                                                                                                                                                   tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCgg                                                                                                                                                                                                      >  8:466149/1‑116 (MQ=255)                                                                                                                                                                                                                                   tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCgg                                                                                                                                                                                                      <  7:466149/116‑1 (MQ=255)                                                                                                                                                                                                                                   tCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTg                                                                                                                                                                                               <  2:105759/123‑1 (MQ=255)                                                                                                                                                                                                                                    cAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCg                                                                                                                                                                                                                                                                      >  4:406163/1‑54 (MQ=255)                                                                                                                                                                                                                                     aGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATg                                                                                                                                                                   <  8:467874/149‑1 (MQ=255)                                                                                                                                                                                                                                        gcgcACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCg                                                                                                                                                                                                                                                                  <  2:406168/54‑1 (MQ=255)                                                                                                                                                                                                                                            aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGCAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                            >  7:260745/1‑149 (MQ=255)                                                                                                                                                                                                                                            aCGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                            <  5:234477/149‑1 (MQ=255)                                                                                                                                                                                                                                               ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                                                                                                                                                       >  5:74316/1‑135 (MQ=255)                                                                                                                                                                                                                                               ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAg                                                                                                                                                                       <  6:74316/135‑1 (MQ=255)                                                                                                                                                                                                                                               ttCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCGTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATg                            <  2:161412/274‑1 (MQ=255)                                                                                                                                                                                                                                                    cctcaccctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATAACcgc                                                                                                                                                    >  5:16406/1‑149 (MQ=255)                                                                                                                                                                                                                                                          cctAACCCTCTCCCTCAAGGGAGAGGGTACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGt                <  2:20756/275‑1 (MQ=255)                                                                                                                                                                                                                                                           ctAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCaa                                                                                                                                             >  5:264178/1‑149 (MQ=255)                                                                                                                                                                                                                                                              aCCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGAGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGAtt                                                                  >  1:53795/1‑221 (MQ=255)                                                                                                                                                                                                                                                               cccTCTCCCTTGAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                                                                                                                                                  <  7:236097/140‑1 (MQ=255)                                                                                                                                                                                                                                                               cccTCTCCCTTGAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGa                                                                                                                                                  >  8:236097/1‑140 (MQ=255)                                                                                                                                                                                                                                                                ccTCTCCCTTGAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACa                                                                                                                                                                                                               <  6:182159/78‑1 (MQ=255)                                                                                                                                                                                                                                                                ccTCTCCCTTGAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACa                                                                                                                                                                                                               >  5:182159/1‑78 (MQ=39)                                                                                                                                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGCGCCAATTGGa                                                                                                                                        <  8:260745/149‑1 (MQ=255)                                                                                                                                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa                                                                                                                                        >  7:356398/1‑149 (MQ=255)                                                                                                                                                                                                                                                                ccTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATACTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGa                                                                                                                                        >  8:458277/1‑149 (MQ=255)                                                                                                                                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                            >  7:469806/1‑127 (MQ=255)                                                                                                                                                                                                                                                                  tctcCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGa                                                                                                                                                            <  8:469806/127‑1 (MQ=255)                                                                                                                                                                                                                                                                    tcCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGC‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTACg                                                                                                                >  1:12496/1‑169 (MQ=255)                                                                                                                                                                                                                                                                    tcCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTAAGACGATTATGATGAATATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAg             >  1:63575/1‑268 (MQ=255)                                                                                                                                                                                                                                                                      ccTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa                                                                                                                                  >  7:23391/1‑149 (MQ=255)                                                                                                                                                                                                                                                                      ccTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCaaa                                                                                                                                  >  6:192461/1‑149 (MQ=255)                                                                                                                                                                                                                                                                        tCAAGGGAGAGGGGACCGATCGAG‑CACAACTTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCTGGCTTTTGCGCCTGAATATGGCGGAGAAAAAGCATGGCCGGGACTGACGGCGAGACAATTAGAAGCAAAGGGAGACAGACGCATACGGTAACACAAAAGGTACGCGGACTGCGTTACGAAAATTACGAAGATTCCGACGTGTATGCTACCTATATCGATGAAAAACGGATTCCGAGCAGAGGAAAgagataagcgag  >  4:498795/1‑264 (MQ=255)                                                                                                                                                                                                                                                                        tCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAgg                                                                                                                                >  8:419722/1‑149 (MQ=255)                                                                                                                                                                                                                                                                        tCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGggcaaaggcaa  <  2:398800/275‑1 (MQ=255)                                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                     TGCTGGAGTTGCTTTCACCTGCCCAAAGGCCATTACAAATCACACGAGATTTGAGCGACTTCTGGAAAGGAGCGTACCGTGAGGTGCAAAAAGAGATGAAAGGGCGTTATCCCAAACATGTCTGGCCGGACGACCCGGCAAATACTGCACCGACGCGACGGACGAAAAAGTATTCGTAAGGTTGGTTTTCTCCCTCTCCCTGTGGGAGAGGGCCGGGGTGAGGGCATCAGCGCGCACGTTCACCCTCACCCTAACCCTCTCCCTCAAGGGAGAGGGGACCGATCGAG‑CACAAATTTTGAGA‑GATATCTTCTTCTGTCTTGTAACAGA‑AGAACAGAAAATCGGGCTTTTGCGCCTGAATATTGCGGAGAAAAAGCATGGCCGGGAATGACCGCGAGCCAATTGGACGCAAAGGGAAACCGACGCGTCCGGTCAAACAAAAGGTAAGCCGTCGTCGTTACGAAGATGACGATGATTACGACGATTATGATGACTATGAGGATGAAGAACCGATGCCGCGCAAAGGTAAGGGCAAAGGCAA  >  NC_000913/164356‑164893

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