Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	92,687	A→G	P425P (CCA→CCG)	ftsI →	peptidoglycan DD‑transpeptidase FtsI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	92,687	0	A	G	96.8%	189.7 / ‑2.2	62	P425P (CCA→CCG)	ftsI	peptidoglycan DD‑transpeptidase FtsI
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (33/27);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

CGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTCGGGCCGGACGGTCGCTA  >  NC_000913/92448‑92932                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         cGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAgc                                                                                                                                                                                                                                               <  1:251837/251‑1 (MQ=255)     aagcggAATTAACCCTGACGGGGGAATTACAGAAGTCGGTTAACTTCGGTTTGTCCTGTGTGACGTTGTCGATGCCGCCCTCTGCGTTAGTGTAGGCTTACTACCGTTTTGGACTGTGAAAAGCGACCAATTGGGGGTTGGTCGGAGAACTCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCAGTTTGTCTTTCTGCTACGGGCTAATGGTTACACCGGTACAGTTAgcgc                                                                                                                                                                                                                                             <  2:163090/244‑1 (MQ=255)           aaTTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGTTAACACCGTTACAGTTAGCGCGAGTCTAc                                                                                                                                                                                                                                     <  1:390107/251‑1 (MQ=255)           aaTTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGCCTGACATAGAGCGGGCCACCTTCTCTTTCGGCTACGGGCTAATGATAACACCGTTCAAGTTAGCGCGAGTCtcc                                                                                                                                                                                                                                     >  2:394048/1‑249 (MQ=255)                 cccTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTAAAAAACAACGGTGGTCTGACATAGAGAGGGCAACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAgcctccccaact                                                                                                                                                                                                                               >  1:416168/1‑240 (MQ=255)                  ccTGACCGGGATATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCCTTCACGATGCCGTCCTCAGCGTTAGCAGATCCTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCCAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGACTAATGGTAACACCGTTACAGTTAGCGCGAATCTACGCAACt                                                                                                                                                                                                                               >  1:58602/1‑250 (MQ=255)                         ggggTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGc                                                                                                                                                                                                                         >  1:354930/1‑249 (MQ=255)                                      aGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTAt                                                                                                                                                                                                          <  1:192030/251‑1 (MQ=255)                                         cGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGc                                                                                                                                                                                                       <  1:336829/251‑1 (MQ=255)                                         cGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGGTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCg                                                                                                                                                                                                        <  2:443005/250‑1 (MQ=255)                                           aGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGccc                                                                                                                                                                                                     <  1:183243/251‑1 (MQ=255)                                            gTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCa                                                                                                                                                                                                    <  2:408864/251‑1 (MQ=255)                                             tAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCAc                                                                                                                                                                                                   >  1:20532/1‑251 (MQ=255)                                                   cGGTGTTTTCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCg                                                                                                                                                                                              <  1:315359/250‑1 (MQ=255)                                                     gtgtTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTCGCGGCAGCTACGGCATTTATCGCCCACTGTCGAtt                                                                                                                                                                                           <  2:190779/251‑1 (MQ=255)                                                              aGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGtt                                                                                                                                                                                  >  1:290921/1‑251 (MQ=255)                                                                gcggCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGGTGGGCGGAGAACGCAGTGGCTTATATCCGCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGa                                                                                                                                                                                <  2:126677/249‑1 (MQ=255)                                                                     gTTAGCGATGCCGTCCTCAGCGTTAGGAGATAATTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGAcccct                                                                                                                                                                           <  2:244778/251‑2 (MQ=255)                                                                       tAGCTATGCCGTCCTCAGCGTTAGTATATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGAcccc‑cc                                                                                                                                                                         <  2:321024/250‑1 (MQ=255)                                                                            aTGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAACGCGACCACTTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACACAAAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGAAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGACAACTACGGCATTTATCGCCCCATGTCGATTACCAACGTT‑ACCCCACCGGTTcc                                                                                                                                                                   >  2:381951/1‑251 (MQ=255)                                                                             tGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTccc                                                                                                                                                                  >  1:158866/1‑251 (MQ=255)                                                                                         cGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTc                                                                                                                                                      >  1:186795/1‑251 (MQ=255)                                                                                                    acttactCACTTTTTGGACTGGGAAAAGCGACCAATTTGGGGTGGGTCGGAGAACGCAGTGGCATATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCTATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAAt                                                                                                                                            <  2:262076/250‑1 (MQ=255)                                                                                                                tttGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATACTTTGTCCGCAc                                                                                                                               >  2:243068/1‑251 (MQ=255)                                                                                                                    gACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGTGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTAGATTACCAAAGTTGACCCA‑CCGCGTCCCGGTGCACGTGTCTTCCCGGAACCCATTGGCCGCACggtg                                                                                                                           >  2:43160/1‑251 (MQ=255)                                                                                                                    gACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGCTTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACcgtg                                                                                                                           >  2:75782/1‑251 (MQ=255)                                                                                                                       tgagAAAAGAGACCAATTTGGGGGTGGTCGGAGCACGAAGTGGCTTATATCCTCAAACACAACGGGGGTCTGACATAGCGAGGGCCCACTTCTCTTTCGCCTACGGGCTAACGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGACGCTTACACCAGTTGCCCCA‑CCGGTTCCCGGTGAACGTGTCTTCCCGGACTCCATTggacgaccagtag                                                                                                                          >  2:198866/4‑237 (MQ=255)                                                                                                                                cGACCAATTTGGGGTTGGTCGGCGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGACCGCACCGTGGGGCATatgatt                                                                                                               >  2:347006/1‑250 (MQ=255)                                                                                                                                cGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGTGCATatgatg                                                                                                               >  2:129521/1‑251 (MQ=255)                                                                                                                                   ccAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCAAATGATGGa                                                                                                             >  2:361780/1‑250 (MQ=255)                                                                                                                                    cAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACAGAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGaaa                                                                                                           <  2:351896/251‑1 (MQ=255)                                                                                                                                     aaTTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGaa                                                                                                            >  1:174636/1‑249 (MQ=255)                                                                                                                                          ggggTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTg                                                                                                      <  1:129521/250‑1 (MQ=255)                                                                                                                                               tGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCt                                                                                                 <  1:449784/250‑1 (MQ=255)                                                                                                                                                   cGGAGAACGCAGTGGCTTATGTCCTCAAAAATAACGGTGGTCTGACATGGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGTGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTAcc                                                                                              >  1:146391/1‑249 (MQ=255)                                                                                                                                                        aaCGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTAAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGgg‑cggc                                                                                      >  2:414171/1‑251 (MQ=255)                                                                                                                                                           gaagtagCTTATGTCCTCAACAACAACGGTGGACTGACATAGAGAGGGCCACATTCTCTTTCGGCTACGAGCTAATGGTCACACCGTTACAGTTAGAGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCCAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGACCGCACCGAGGTGCACATGATGGCAAACGAGGCGCGACCAgcgggggac                                                                                    >  2:452830/7‑243 (MQ=255)                                                                                                                                                                    tataTCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGATACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGa                                                                            >  1:113239/1‑249 (MQ=255)                                                                                                                                                                         ccTCAAAAACAACGGTGGTCTGAAATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAggcgg                                                                      <  1:96474/250‑1 (MQ=255)                                                                                                                                                                          cacAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCGTCCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCGCCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGgcggcg                                                                    <  2:158866/249‑1 (MQ=255)                                                                                                                                                                              aaaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTa                                                                <  1:41040/251‑1 (MQ=255)                                                                                                                                                                               aaaCAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCCACTATCGGCAGCTACGGCCTTTATCGCCCCCTGTCGCTTACCAAAGTTCACCCC‑CCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGTCCGCCCCGTGGGGCATTTGATGGAAAGCGGGGCGCTACCGGG‑ATGCGGCGGCGTGAcggcggggctta                                                                >  2:90969/1‑244 (MQ=255)                                                                                                                                                                                    acGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGACCGTGTCTTCCCGGAATCCATTGACCGCACCGTGGTGCCTATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCt                                                          >  2:446989/1‑251 (MQ=255)                                                                                                                                                                                         ggtCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCtatcgt                                                     >  1:364789/1‑251 (MQ=255)                                                                                                                                                                                                     atagGGCCACCCTCTCTTTCGGCTACGGGCTAATGGTAACCCCGTTACAGTTAGCGCGAGTCTACGCAACTATCGTCAGCTACGGCATTTAGCGCCCACTGGCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaa                                         <  2:352400/249‑1 (MQ=255)                                                                                                                                                                                                      gagGGCCACCTTCTCTTTCGGCTACGCGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTGCCCGGTGAACGTGTCTGCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTaaaa                                        <  1:19411/251‑1 (MQ=255)                                                                                                                                                                                                          gCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTATACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGCCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCgg                                    >  2:406496/1‑251 (MQ=255)                                                                                                                                                                                                           ccACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGAGAGTCTAAGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCACGGTGAAAGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGAGCCTATGATGGAAAGCCTGGCGCT‑CCAGGGAGGCGGCGGCGTGAAGGCGGCGATTAAAGGAGACCGTATCACCCTTAAAAcagg                                    >  2:133663/1‑247 (MQ=255)                                                                                                                                                                                                             cccTTCTCTTTCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTAc                                 <  1:414171/250‑1 (MQ=255)                                                                                                                                                                                                                     tttCGGCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCtt                                                                                                                                                    >  1:237392/1‑129 (MQ=255)                                                                                                                                                                                                                          gCTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTc                   <  2:362128/251‑1 (MQ=255)                                                                                                                                                                                                                           cTACGGGCTACTGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAGGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAGAGGCTATCGTATCGCCCTTAAAACCGGTACC‑GCGAAACAGGTCg                  >  2:4484/1‑251 (MQ=255)                                                                                                                                                                                                                           cTACGGGCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTCg                  >  1:314107/1‑251 (MQ=255)                                                                                                                                                                                                                               gggCTAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTCGGGcc              <  1:363148/251‑1 (MQ=255)                                                                                                                                                                                                                               gggCTAAGTGTAACACCGGTACTGTTAGCGCGAGGCTACGCAACTATCGTCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTCGGGcc              <  2:20532/251‑1 (MQ=255)                                                                                                                                                                                                                                 gCTAATGGTAACAACGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTCTCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCACGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGAATAAAGGCTATCGTATCGCCATTAAAACCGGTACA‑GCGAAAAAGGTCGGGCCgg            >  1:62683/1‑251 (MQ=255)                                                                                                                                                                                                                                   tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGCTTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGGGGTGCATATGATGGAAAGAGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTAGCGTATCGCCATTAAAACCGG‑ACCAGCGAAAAAGGGCGGGCCGGAc          >  2:1803/1‑251 (MQ=255)                                                                                                                                                                                                                                   tAATGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGCATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTCAAACCGGTACC‑GCGGAACAGGGCGGGCCgtcc          >  2:206291/1‑248 (MQ=255)                                                                                                                                                                                                                                      tGGTAACACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCACACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGGGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTGCCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGCACC‑GCGAAAAAgggcgggcagggaggg       >  2:175186/1‑243 (MQ=255)                                                                                                                                                                                                                                          aaCACCGTTTCAGTTAGCGCGAGTCTACGCAACTACCGGCAGCTACGGCATTTATCGCCCACTGTCGATTTCCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTACC‑GCGAAAAAGGTCGGGCCGGACGGTCGCt   <  2:354930/251‑1 (MQ=255)                                                                                                                                                                                                                                          aaCACCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGTGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTCGGGCCGGACGGTCGc    <  1:160286/250‑1 (MQ=255)                                                                                                                                                                                                                                           acacCGTTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAATTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACCGTGGTGCATATGATGGAAAGCGGGGCGCTACCGGG‑CGGCGGCGGCGTGAAGGCTGCGATTAAAGGCTATCGTATCGCCATTAAACCCGGTACC‑GCGAAAAAGGTAGGGCCCGACgggccctc  >  2:179159/1‑245 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                         CGCTACAGCGAATTAACCCTGACCGGGGTATTACAGAAGTCGAGTAACGTCGGTGTTTCCAAGCTGGCGTTAGCGATGCCGTCCTCAGCGTTAGTAGATACTTACTCACGTTTTGGACTGGGAAAAGCGACCAATTTGGGGTTGGTCGGAGAACGCAGTGGCTTATATCCTCAAAAACAACGGTGGTCTGACATAGAGAGGGCCACCTTCTCTTTCGGCTACGGGCTAATGGTAACACCATTACAGTTAGCGCGAGTCTACGCAACTATCGGCAGCTACGGCATTTATCGCCCACTGTCGATTACCAAAGTTGACCCC‑CCGGTTCCCGGTGAACGTGTCTTCCCGGAATCCATTGTCCGCACTGTGGTGCATATGATGGAAAGCGTGGCGCTACCAGG‑CGGCGGCGGCGTGAAGGCGGCGATTAAAGGCTATCGTATCGCCATTAAAACCGGTACC‑GCGAAAAAGGTCGGGCCGGACGGTCGCTA  >  NC_000913/92448‑92932

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