Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,592,188	C→T	H178Y (CAT→TAT)	tsr →	methyl‑accepting chemotaxis protein Tsr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,592,188	0	C	T	100.0%	227.1 / NA	69	H178Y (CAT→TAT)	tsr	methyl‑accepting chemotaxis protein Tsr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (35/34);  total (35/34)

TGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCG  >  NC_000913/4591949‑4592426                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                tGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGcc                                                                                                                                                                                                                                      >  1:102561/1‑251 (MQ=255) tGAAACAGGCGGAAAAAAACTGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGcc                                                                                                                                                                                                                                      >  2:153144/1‑251 (MQ=255)  gAAACAGGCGGAAAAAAACTGGGCAGAGTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGcc                                                                                                                                                                                                                                      <  1:428489/250‑1 (MQ=255)           ggAAAAAAACGGGGCAGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATa                                                                                                                                                                                                                            <  1:466775/251‑1 (MQ=255)                        gcagATTACGAAGCGTTGCCGCGTGCCCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTa                                                                                                                                                                                                               <  2:510207/248‑1 (MQ=255)                        gcagATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGcctcct                                                                                                                                                                                                                >  2:157156/4‑250 (MQ=255)                          agATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGAGATCGCCGTCAGCGATAACAATGCCTCCTACa                                                                                                                                                                                                             <  1:154294/250‑1 (MQ=255)                              tACGAAGCGTTTCCGCGTGACCCGCGTCATAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTAGTAGGTGCAGGCAAAATCAACGAGTTCTGTGAGCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATTATTTCGGCGTCAGCGATAGCAATGCCTCCTACAGcc                                                                                                                                                                                                          <  2:597998/250‑1 (MQ=255)                                cGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAgg                                                                                                                                                                                                       <  1:299895/251‑1 (MQ=255)                                 gAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGc                                                                                                                                                                                                      >  1:50476/1‑251 (MQ=255)                                   aGCGTTGCCGCGTGACCCGCGACAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCg                                                                                                                                                                                                     >  2:157704/1‑250 (MQ=255)                                                gACCCGCGTCAGCGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCCGGCGATGTGGATTCtggt                                                                                                                                                                                       >  1:131277/1‑251 (MQ=255)                                                         cAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGtga                                                                                                                                                                               <  1:79549/250‑1 (MQ=255)                                                              cACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACCATGCGCTGGCGGCGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCCGTATGTGGCTTCAATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAAACATGCCTCCTACAGCCAGGCGATGTGGATTCTGTGGGGCGTGGTGAtcg                                                                                                                                                                         >  2:228277/1‑251 (MQ=255)                                                                ccgccgCAGCGGCAGAGATCAAACGTAATTACGATATTGATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGAGTCTGGTGGGCGTGATGAtcgtc                                                                                                                                                                       <  1:238336/251‑1 (MQ=255)                                                                       aGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTgg                                                                                                                                                                <  2:472351/251‑1 (MQ=255)                                                                        gCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATTCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTTCTGGc                                                                                                                                                               <  2:324117/251‑1 (MQ=255)                                                                               agaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATc                                                                                                                                                        >  1:269891/1‑251 (MQ=255)                                                                                gaTCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATTGGGATTCTGGTGGGCGTGATGAGCGTCGTACTGGCGGTCATCt                                                                                                                                                       <  1:238875/251‑1 (MQ=255)                                                                                   cAAACGTAATTACGATATTTATCCCAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGAAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTGCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCg                                                                                                                                                    <  2:23744/251‑1 (MQ=255)                                                                                   cAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCtttc                                                                                                                                                    >  1:240062/1‑249 (MQ=255)                                                                                   cAAACGTAATTACGATAGTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGTGCGGGATGATCGTCGTACTGGCGGTCATCTTc                                                                                                                                                     <  2:56510/250‑1 (MQ=255)                                                                                           aTTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGCTTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCAttttacgcgtctgg                                                                                                                                            >  2:187570/1‑238 (MQ=255)                                                                                            ttACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGCTCGTCGTACTGGCGGGCATCTTCTCCCTCTggg                                                                                                                                           >  2:87113/1‑250 (MQ=255)                                                                                              aCGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGtt                                                                                                                                          <  1:256789/250‑1 (MQ=255)                                                                                                 atatTTATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGt                                                                                                                                      <  2:50476/251‑1 (MQ=255)                                                                                                     ttATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTTCTCTTCGCCGTCTgggttgggattt                                                                                                                                  >  2:472419/1‑241 (MQ=255)                                                                                                     ttATCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCTCCGCCTGGGTTGGTAtt                                                                                                                                   >  1:466731/1‑250 (MQ=255)                                                                                                        tCACAATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTCCCCGTCTGGTTTGGTAtttaa                                                                                                                                >  1:17598/1‑247 (MQ=255)                                                                                                           caATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGccc                                                                                                                            >  1:493372/1‑250 (MQ=255)                                                                                                           caATGCGCTGGCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGTTCGGGGTTAAAGccc                                                                                                                            >  1:231802/1‑250 (MQ=255)                                                                                                                cgcTGGAGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTCTCCGTCTGGGTTGGGATTAAAACCTCGCTg                                                                                                                       >  2:422240/1‑251 (MQ=255)                                                                                                                   gggCGGAGCTGAGCCAACTATTAGGTGCATGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGATCAAAACGATCGGCTCTATGATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTTGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTGCGGTATTAAAGCCTCGCTGGTa                                                                                                                    <  2:286478/250‑1 (MQ=255)                                                                                                                    ggCGGAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAg                                                                                                                   <  1:153144/251‑1 (MQ=255)                                                                                                                       ggAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGGCTGGTTCGGTATTAAAGCCTCGCCTGTAgcgc                                                                                                                >  1:385023/1‑251 (MQ=255)                                                                                                                       ggAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGGCATCTTTGCCGGCTGGGTCGGTATTTAAGACTCGCTGGtgtcgc                                                                                                                >  2:599254/1‑246 (MQ=255)                                                                                                                        gAGCTGATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGcc                                                                                                               <  2:170758/251‑1 (MQ=255)                                                                                                                            ggATCCAACTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATg                                                                                                           <  1:432888/250‑1 (MQ=255)                                                                                                                                  aaCTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGc                                                                                                     >  1:370288/1‑251 (MQ=255)                                                                                                                                  aaCTATTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTTGGTCTTTAAACCTCGCTGGTAAGGCCCATGAATCGc                                                                                                     >  2:140764/1‑251 (MQ=255)                                                                                                                                      attAGGTGCAGGGAAAATCAACGAGTTCTGTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACAGGGAGCAAAACGATCGGCTCTATGATAGCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCGGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGa                                                                                                 <  2:102561/250‑1 (MQ=255)                                                                                                                                          ggTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGa                                                                                             <  2:38589/251‑1 (MQ=255)                                                                                                                                               aGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCa                                                                                        <  2:201371/251‑1 (MQ=255)                                                                                                                                                   aaactCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCCTCGCCGTCTGGTTCCGTCTTAAAGGCACGCTGGGAGCGCCCATGTATCACCTGGTTgtaaccattcc                                                                                    >  1:277574/5‑241 (MQ=255)                                                                                                                                                   aaaaTCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCCTCTGGTTCGGTATTAAAACCTCGCTTGGAGCGCCCATGGATCGCCTGTTTggaagccttcc                                                                                    >  1:277582/1‑241 (MQ=255)                                                                                                                                                        cAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTAGGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCat                                                                                <  2:71479/250‑1 (MQ=255)                                                                                                                                                                   ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAggcgg                                                                     <  1:521816/250‑1 (MQ=255)                                                                                                                                                                   ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCGGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAggc                                                                       <  2:469369/248‑1 (MQ=255)                                                                                                                                                                   ttGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCCCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACCACATTCTTCATATTGCTggcggc                                                                    >  2:319479/1‑251 (MQ=255)                                                                                                                                                                      aTCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCACCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTTGGTGTTAAAGCCTCGCTGGGTGCGCCAATGAATCGCCTGATTGACAACATTTGTCATTATGCAGGGGGCGAt                                                                 >  1:515903/1‑251 (MQ=255)                                                                                                                                                                             gACCCAGGGATATCAGGCCGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTCTGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATTTTCGCCGTCTGGGTCGGGATTAAAGCCTCGCTGGTAGCGCCCATGAAACGCCTGCTTTACAGCAATCGTCATCTTGCAGGCGGCGacattgggt                                                          >  2:360720/1‑243 (MQ=255)                                                                                                                                                                              aCCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGGATTTAAGCCTCGCTGGTAGCGCCCATGAATCGCCTGATTGACAGCATTCGTCATTTTGCCGGCGGGGGTCTGGTGa                                                          >  1:395595/1‑250 (MQ=255)                                                                                                                                                                                     gATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTCGGTATTAAAGCCTCGCTGGGAGCGCCCATGAATCGCCTGATTGAAAGCATTCGTCATATTTCAGGCGGCGATCTGGTGAAACCCAtt                                                  >  1:304472/1‑251 (MQ=255)                                                                                                                                                                                        atCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGGTCGTCGTACTGGCGGTCATCTTCGCCCTCTGGGTCGGTATTAAAGCCTCGCTGGGAGCGCCCATGGATCACCAGCTTGGCAGAATTCGTCATCTTTCAGGCGGGGATATTGGGAAACCGATTg                                                 >  2:175725/1‑249 (MQ=255)                                                                                                                                                                                                ggTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGCTATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGCCATCTTCGCCGGTTGGTTCGGTATTACAGCCTCGCTGGTAGCGCCAATGAATCGACTGTTTGTCAACATTCGGCATATTTCCGGCGGCGATCTGGTGAAAACCATTGCGGTGG‑ATgg                                      >  1:587677/1‑251 (MQ=255)                                                                                                                                                                                                     cGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTa                                 <  2:583728/251‑1 (MQ=255)                                                                                                                                                                                                       agaAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCTCCGTCTGGGTCGGTATTAAAACCTCGCTGGTAGCGCCCATGAATCGCCTGATTTACAGCATTCCTCATATTGAAGGCGGCGATCTGGTGGAACCCATTTAgggggggtgctctaat                                >  1:539320/1‑237 (MQ=255)                                                                                                                                                                                                       agaAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTTTGGGTTGGTATTTAAGACTCGCTGGGAGCGCCCATGAATCGCCTGATTGACAGCATTCGTCGTATTTCCGGCGGCGGTCTGGTGGAACCGATTGAGGGGG‑ATGGGTCTAAt                               >  1:550513/1‑251 (MQ=255)                                                                                                                                                                                                          aGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTAc                                                                                                                                                                   <  1:482673/117‑1 (MQ=255)                                                                                                                                                                                                          aGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTAc                                                                                                                                                                   >  2:482673/1‑117 (MQ=255)                                                                                                                                                                                                             aGTATGTGGCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTTGCCGTCTGGGTCGGGGTTAAAACCTCGCGGGGAGGGGCAATTGATCGCCCGGGTTACAGCATTCTTCATATTGCAGGGGGGGCTCTGGTGGAACAGGATGAGGGGG‑GTGGGTCTAATgggagg                         >  2:11815/1‑246 (MQ=255)                                                                                                                                                                                                                    ggCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAAc                  <  2:179213/251‑1 (MQ=255)                                                                                                                                                                                                                    ggCTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCAGCTTCGCCGTCTGGTTCGGTATTAAAGCCGCGCTGGTAGCGCCAATGAATCGCCTGATTGAc                                                                                            <  2:314309/178‑1 (MQ=255)                                                                                                                                                                                                                      cTTACATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGGTTGGTATTAAAACCGCGCCGGTAGCGCCAATAGATCCCCTGATTTGCCGCATTCGTCATATTGCAGGCGGCGGTCCGGGGAAACCCATTGAGGGGGGATG‑CTCTAATGAGATGTGGccaccg                >  2:101687/1‑246 (MQ=255)                                                                                                                                                                                                                         aCATGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCGTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAACTGGCa             <  2:281010/251‑1 (MQ=255)                                                                                                                                                                                                                            tGGAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTTGGCGTGATGATCGTCGTACTGGCGGTCCTCTTCGCCTTCTGTTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCATGCGGCGATCTGGTGAAAcc                                                      <  2:79782/208‑1 (MQ=255)                                                                                                                                                                                                                             ggAGCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAACTGGCagag          <  2:478317/250‑1 (MQ=255)                                                                                                                                                                                                                                gCAAAACGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGtt      <  1:422240/251‑1 (MQ=255)                                                                                                                                                                                                                                    aaCGATCGGCTCTATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTgcg  <  1:230588/251‑1 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                TGAAACAGGCGGAAAAAAACTGGGCGGATTACGAAGCGTTGCCGCGTGACCCGCGTCAGAGCACCGCCGCAGCGGCAGAGATCAAACGTAATTACGATATTTATCACAATGCGCTGGCGGAGCTGATCCAACTGTTAGGTGCAGGCAAAATCAACGAGTTCTTTGATCAGCCGACCCAGGGATATCAGGACGGTTTCGAGAAGCAGTATGTGGCTTACATGGAGCAAAACGATCGGCTCCATGATATCGCCGTCAGCGATAACAATGCCTCCTACAGCCAGGCGATGTGGATTCTGGTGGGCGTGATGATCGTCGTACTGGCGGTCATCTTCGCCGTCTGGTTCGGTATTAAAGCCTCGCTGGTAGCGCCAATGAATCGCCTGATTGACAGCATTCGTCATATTGCAGGCGGCGATCTGGTGAAACCGATTGAGGTGG‑ATGGCTCTAATGAGATGGGGCAACTGGCAGAGAGTTTGCG  >  NC_000913/4591949‑4592426

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