Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,993,531	T→C	R793R (CGT→CGC)	cyaA →	adenylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,993,531	0	T	C	98.6%	239.4 / NA	69	R793R (CGT→CGC)	cyaA	adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (34/34);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCG  >  NC_000913/3993290‑3993776                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      aCAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTGCCGGCGGTGGTGGACGGCGTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATAGCACCACTGCGAAGGCAGCAAAGAGGAGTTGGTACGTGACGTCAGTCGCTGCTACTCGTCATCGCATGACCGCTTGACCTAc                                                                                                                                                                                                                                              <  2:222523/251‑1 (MQ=255)         cACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGGGGGGGACGGCTTTGCCAGCGAAGGGATCATCCCGTTCTTTTTCGGAGAAA‑CGCAAGACGAGAATGGCCTTAAAATCTACATTCTCGACGAAAACAACCCGGGTGAGGGATATCACCCCTGCGGAGGCAGCAAAAAAGAGCTGGGACGTGACGTCAGTCGCCTCTACTCGTCATCGCATTACCGCTTTTCCTACCGCTCaa                                                                                                                                                                                                                                       >  2:75734/1‑250 (MQ=255)             gCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGGGGAGGGCATTGCCAGCGCAGGGATCAACCAGTTCCTTTTCGCAGAAA‑CGGAAAAAGAGAAAGGATTTAATATCTTCATACTCGACGAAAACAACCGGGGTGAAGGAAATAACCACTGGGGAGGCAGCAAAAAGGAGCTGGTACGCGACGTCAGTCGCTTCTACTCGTCATCGCAAGACAGCTTTACCAACGGCTCCAGcatt                                                                                                                                                                                                                                  >  1:196392/1‑248 (MQ=255)               cTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGGCGGCTTTTCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAAACCCAA‑ACAAGAATGGCTTTAATATCTAACTTCTCCACGAAAGCAAACGGGGTTAGGGATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTAAAGGACGTCAAGCGCTTTTACTCCTCATCGCATGACCCCGTTACCTACGGCTCCAGCTtcat                                                                                                                                                                                                                                >  1:221496/1‑251 (MQ=255)                           gTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTAc                                                                                                                                                                                                                                              <  1:136426/225‑1 (MQ=255)                             tGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTg                                                                                                                                                                                                                   <  1:45170/250‑1 (MQ=255)                                                     ggcggTGGTGGACGGCTTTGCCAGCGAAGGGATCCTCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCa                                                                                                                                                                                                                        >  1:509226/1‑221 (MQ=255)                                                     ggcggTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTGACCTACGGCTCAAGCTTCATCAACTTCa                                                                                                                                                                                                                        <  2:509226/221‑1 (MQ=255)                                                       cggTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGtt                                                                                                                                                                                        >  1:225454/1‑251 (MQ=255)                                                        ggtggtGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTg                                                                                                                                                                                       >  2:93010/1‑251 (MQ=255)                                                               gACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGAATATCACCACTGCGAAGGCGGCCCTGAGGAGCTGGAACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTAACTACGGCTCAAGCTTAATCAACTTCAACCTGCCGCATATCACTCTGATTGTGAAGGTTGATGGTct                                                                                                                                                                                >  2:277272/1‑250 (MQ=255)                                                                aCGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGAAGGTCGt                                                                                                                                                                               <  1:142115/251‑1 (MQ=255)                                                                 cGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGACGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGACGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTg                                                                                                                                                                              >  2:39826/1‑251 (MQ=255)                                                                        gCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGt                                                                                                                                                                       <  2:375071/251‑1 (MQ=255)                                                                            gCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTTAACAGGTGAtt                                                                                                                                                                   <  2:470667/251‑1 (MQ=255)                                                                            gCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGAtt                                                                                                                                                                   >  1:470667/1‑251 (MQ=255)                                                                                           aattCTTTTTCGAAGAAA‑CGCAAGACGAGAATGACTTTAATATCTACATGCCCGACGAAAGCAACCTGGTTGAGGTATATCACCACTGCCAAGGCAGTAAAGATGATCTTGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAGGACCGCTTTACCTTCGTCTCAAGCTTCTTCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATTATCGTGAATAGTTGATTCCTTTGCGCACaa                                                                                                                                                      <  2:498342/247‑1 (MQ=255)                                                                                               cTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCt                                                                                                                                                 >  1:486682/1‑250 (MQ=255)                                                                                                    tCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCgg                                                                                                                                            <  2:179528/250‑1 (MQ=255)                                                                                                     cGAAGATA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGt                                                                                                                                           >  1:144370/1‑250 (MQ=255)                                                                                                     cGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGt                                                                                                                                           >  1:268491/1‑250 (MQ=255)                                                                                                          aaa‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGCCGAAACCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCACAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATg                                                                                                                                     >  1:261566/1‑251 (MQ=255)                                                                                                          aaa‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATg                                                                                                                                     >  1:127497/1‑251 (MQ=255)                                                                                                            a‑CGCAAGACGAGAATGGCTTTAATTTTTACATTCTCGACGAAAGCAGCCGGGTTGAGGTATATCACCATTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGTTTCGACTCTTCATCGCATGACCGCTTTACCTACGGCGCAAGCTTCATCAACGTCAACCTGCCGCAGTTCTATCAGATGGTGAATGTTGATGTGCGTGAACAGGGGATGCCGTTGCGCACAAAATCTATCGGTAACATgcc                                                                                                                                   <  2:324692/251‑1 (MQ=255)                                                                                                                      gagaATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCa                                                                                                                           >  1:246684/1‑250 (MQ=255)                                                                                                                         aaTGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGAAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCGCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGAGGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCa                                                                                                                       <  2:489812/251‑1 (MQ=255)                                                                                                                                 tAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACg                                                                                                               >  1:184580/1‑251 (MQ=255)                                                                                                                                  aaTATCTACATTCTCGACGAAAGCAGCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGCCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGTACAGGTGATTCCGTTCCGCACAAAATCTATCGGAAACATGCCGCCTGCCAATCAGGATCACg                                                                                                               >  1:224715/1‑250 (MQ=255)                                                                                                                                  aaTATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACg                                                                                                               <  1:39826/250‑1 (MQ=255)                                                                                                                                         aCATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAc                                                                                                           <  1:75734/247‑1 (MQ=255)                                                                                                                                          cATTCTCGACGAAAGCAACCCGGTTGAGGAATCTCACCCCTGCGAAGGCAGCCAAGAGAAGCTGGTACGTGACGTCAGTCGCTTCTACTCGCCATCGCATGACCGCTTTACCTACGGCCCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATAcgccgc                                                                                                      >  2:211434/1‑251 (MQ=255)                                                                                                                                             tcccgacgaAAGCAACCGGGTTGAGGTATATCTCCACTGCGAAGGCAGCAAAGAGGAGCTGTTACGTTACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTTCCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTa                                                                                                    <  2:374720/247‑1 (MQ=255)                                                                                                                                              ctcGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTAtt                                                                                                  >  2:241285/1‑251 (MQ=255)                                                                                                                                                cgacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTa                                                                                                 <  2:144370/250‑1 (MQ=255)                                                                                                                                                cgacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTa                                                                                                 <  2:486682/250‑1 (MQ=255)                                                                                                                                                cgacgaAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGGTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTAt                                                                                                   <  2:82674/248‑1 (MQ=255)                                                                                                                                                 gacgaACGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCCAAGAGGAGCTGGTACGTGACGACCGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCCCGATACGCCGCTAttcc                                                                                                >  2:370970/1‑248 (MQ=255)                                                                                                                                                 cgtcctgtcttgcGGGTTGAGGTAGATCCCCCCTGCGATGGCAGCAAAGAGGAGCTGGTACGTGGCGTCAGTCTCTTCTCCTCGTCATCGCATGACCGCTTTGCCTACGGCTCAAGCTTTAGCAACTTCAACCTGCCGCAGTGCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACTAAATCTATCGGTAACATGCCGCCTGCCAATCTGGATCACGATACGCCGCTATTACa                                                                                               <  2:59871/239‑1 (MQ=255)                                                                                                                                                   tccgagcgaCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGGCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACa                                                                                               <  2:59861/241‑1 (MQ=255)                                                                                                                                                    gaAAGCTACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGc                                                                                             >  1:152693/1‑250 (MQ=255)                                                                                                                                                        gCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCGGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGGCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACAGGCCGCCTGCCAATCAGGATCACGGTACGCCGCTATTACATCAatat                                                                                        <  1:207646/251‑1 (MQ=255)                                                                                                                                                           acagggTAGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACACGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGAACACGATACGCCGCTATTACa                                                                                               >  1:59871/4‑241 (MQ=255)                                                                                                                                                           aCCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACa                                                                                               >  1:59861/1‑241 (MQ=255)                                                                                                                                                                   gAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgat                                                                              <  1:190119/250‑1 (MQ=255)                                                                                                                                                                    aGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatg                                                                             >  2:146811/1‑250 (MQ=255)                                                                                                                                                                     ggTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGtgatga                                                                            >  1:331881/1‑250 (MQ=255)                                                                                                                                                                          atacccacTGCGAAGGTAGCAAAGAGGAGCTGGTACGGGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCGCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGATATTACAGCAATATTTGTCGTGAGGATCGTg                                                                       <  2:196392/246‑1 (MQ=255)                                                                                                                                                                              ccacTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTTCGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGAGACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGa                                                                  <  2:223115/251‑1 (MQ=255)                                                                                                                                                                                      aaGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAgg                                                          >  1:486534/1‑251 (MQ=255)                                                                                                                                                                                          cagcaAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACGTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTGCGTGATGAACGTGCCTGAAAGCGAGGCTTa                                                      <  2:492309/251‑1 (MQ=255)                                                                                                                                                                                            gcaAAGAGGAGCTGGTACGTGATGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATc                                                    <  1:188551/251‑1 (MQ=255)                                                                                                                                                                                             caAAGAGTAGCTGGTCCGTGCCGTCAGGCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGGGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCTCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGATCGTGCCGTAAAGCGAGGCTTATcc                                                   <  2:447793/251‑1 (MQ=255)                                                                                                                                                                                             caAAGAGGAGCTGGTACGTGACGGCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTTATCAACTTCAACCTGCCGCAGTTCTATCAGATGGTGAAGGTTGAGGGTCGTGAACAGGTGAGTCCGTTCCGCACAAAATCTTTCGGTAACATGCCGCCTGCCAATCAGGATCCCGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCTAGGCTTAtac                                                   <  2:138972/251‑3 (MQ=255)                                                                                                                                                                                                  gggagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCAGGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGAGTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCAt                                              <  2:152693/250‑1 (MQ=255)                                                                                                                                                                                                    gagCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCAGCATGc                                            <  1:146811/251‑1 (MQ=255)                                                                                                                                                                                                         ggTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGGTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAAt                                         <  2:197953/249‑1 (MQ=255)                                                                                                                                                                                                         ggTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTAACGCTTCATAAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTAGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGACAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCAGGCATGCAATCt                                       >  1:442552/1‑251 (MQ=255)                                                                                                                                                                                                           tACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaa                                                                                                                                                     <  1:192852/139‑1 (MQ=255)                                                                                                                                                                                                           tACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACaaa                                                                                                                                                     >  2:192852/1‑139 (MQ=255)                                                                                                                                                                                                                     cAGTCGCTTCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCCTGCAATCTTAGCGGAAactg                           >  2:165429/1‑251 (MQ=255)                                                                                                                                                                                                                             tCTACTCTTCATCGAATGACCGCTTTACCTGCGGCTCAAGCTTCATCAACTTCAACCGGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGGGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                     <  2:216800/249‑1 (MQ=255)                                                                                                                                                                                                                             tCTACTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGtt                     <  1:93010/249‑1 (MQ=255)                                                                                                                                                                                                                                aCTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCa                  <  2:415583/249‑1 (MQ=255)                                                                                                                                                                                                                                 cTCGTCATCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCAccc               >  1:92386/1‑251 (MQ=255)                                                                                                                                                                                                                                        tCGCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCctgct        >  1:229339/1‑251 (MQ=255)                                                                                                                                                                                                                                          gCATGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCtgctgc      >  2:269319/1‑251 (MQ=255)                                                                                                                                                                                                                                            aTGACCGCTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGt    >  1:297641/1‑251 (MQ=255)                                                                                                                                                                                                                                              gACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCGGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGAGCACGAGACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAAGCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTc   <  1:165429/250‑1 (MQ=255)                                                                                                                                                                                                                                              gACCGCTTTTCCTACGGCTCAAGCTTAATCTACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGACTGCCCATCCGGATCACGATCCGCCGCTATTACAGCAATATATTTCGTGATGAACGTGCCGGAAAGCGAGGATTCTCAGGCATACAATCTTATCGGACACTAACAGTTTCACCCACctcatgcttcg  >  2:279615/1‑242 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                      ACAAACTGCACGGCCTGTCAGTGCAGGTTGAAACCAATCACGTCAAATTACCGGCGGTGGTGGACGGCTTTGCCAGCGAAGGGATCATCCAGTTCTTTTTCGAAGAAA‑CGCAAGACGAGAATGGCTTTAATATCTACATTCTCGACGAAAGCAACCGGGTTGAGGTATATCACCACTGCGAAGGCAGCAAAGAGGAGCTGGTACGTGACGTCAGTCGCTTCTACTCGTCATCGCATGACCGTTTTACCTACGGCTCAAGCTTCATCAACTTCAACCTGCCGCAGTTCTATCAGATTGTGAAGGTTGATGGTCGTGAACAGGTGATTCCGTTCCGCACAAAATCTATCGGTAACATGCCGCCTGCCAATCAGGATCACGATACGCCGCTATTACAGCAATATTTTTCGTGATGAACGTGCCGGAAAGCGAGGCTTATCCGGCATGCAATCTTAGCGGAAACTGACTGTTTCACCCGCCTGCTGCGTCG  >  NC_000913/3993290‑3993776

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