Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,770	C→T	intergenic (+94/‑146)	rho → / → rfe	transcription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,770	0	C	T	100.0%	183.0 / NA	56	intergenic (+94/‑146)	rho/rfe	transcription termination factor Rho/UDP‑N‑acetylglucosamine‑‑undecaprenyl‑phosphate N‑acetylglucosaminephosphotransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/35);  total (21/35)

GACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTA  >  NC_000913/3967532‑3968014                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      gACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTg                                                                                                                                                                                                                                           <  1:391892/251‑1 (MQ=255)         aagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGAGGAGTTTTGCTTTTATATCTGTAATATTAATGCCGCGCTGGGGATGTTAGGAAAGTTCTTGGAATTTGCTGGCATGttt                                                                                                                                                                                                                                   >  2:478471/1‑250 (MQ=255)                  aGAAAATGTGGATCCTGCGCAAAATCATTCACCCGAGGTGCGAAATCGAGGCAATGGAATTCCACATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGTAACGCTCATAAATTTGTCGTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGATTTTAGATCTGTAATCTTAATGCCGCGCTGGCTATGTTAGGAAAATTCTTGAAATTTGCTGGCATTTTATGCAATTTg                                                                                                                                                                                                                          <  2:588148/251‑1 (MQ=255)                   gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGTAATTCCTCATTAAGAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTAGGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATAGCTTTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGc                                                                                                                                                                                                                         <  1:307449/251‑1 (MQ=255)                            gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATATGCA‑TATCAAAt                                                                                                                                                                                                               <  1:109900/251‑1 (MQ=255)                               ccTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATgg                                                                                                                                                                                                             <  2:227231/250‑1 (MQ=255)                                 tGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGCCGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAACCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGACTTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGtt                                                                                                                                                                                                           >  2:331129/1‑250 (MQ=255)                                 tGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGtt                                                                                                                                                                                                           >  2:113213/1‑250 (MQ=255)                                   cgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCAAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAAt                                                                                                                                                                                                        >  1:166124/1‑251 (MQ=255)                                      aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAAGTACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTTTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAAATCTTGGAATTTGCTGGCGTGTTATGCAATTTGCA‑TATCAAATGGTTAAtttt                                                                                                                                                                                                     <  2:34066/251‑1 (MQ=255)                                               cACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAgg                                                                                                                                                                                             >  1:189932/1‑250 (MQ=255)                                                        ggCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGt                                                                                                                                                                                   >  2:530899/1‑251 (MQ=255)                                                           gAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TTTCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTg                                                                                                                                                                                >  2:117670/1‑251 (MQ=255)                                                             aaTCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGa                                                                                                                                                                              <  1:628303/251‑1 (MQ=255)                                                              aTCGATGCAATGGAATTCCTCATTAAGAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGaa                                                                                                                                                                             <  2:16225/251‑1 (MQ=255)                                                                     cAATTGAATTCCTTATTAATAAACTGGCAATGACCAAGACCAATGACGATTTATTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCGTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGTATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTGTGCACAGGACTGGTGGGTTTGGAATGGAcgt                                                                                                                                                                      <  2:237202/251‑3 (MQ=255)                                                                               ccTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGa                                                                                                                                                            <  1:39375/251‑1 (MQ=255)                                                                                      aataaACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGTAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAgg                                                                                                                                                     <  1:83599/251‑1 (MQ=255)                                                                                                cAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGt                                                                                                                                           <  1:486121/251‑1 (MQ=255)                                                                                                 aaTGACCAAGACCAATGACGATTTCTTCGAACTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGAGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGCTGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGtt                                                                                                                                          >  2:238177/1‑251 (MQ=255)                                                                                                    gaccaagaccaaTGGCGATGTCTTCGAAATGATGAAACGCTCATAAAGTTGTCTTATGCCAAAAACGCCACTTGTTGACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTATGAAAATTCTTGGAATTTGCTGGCATGTTAGGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACTGACTTTCCCTTCTTAAGAAAGGTCTTCGGGGTtata                                                                                                                                       <  2:402398/251‑1 (MQ=255)                                                                                                      ccaagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACt                                                                                                                                     <  2:559426/251‑1 (MQ=255)                                                                                                        aagaccaaTGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTc                                                                                                                                   <  2:156022/251‑1 (MQ=255)                                                                                                                       ttcttcGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTtct                                                                                                                    >  2:587317/1‑251 (MQ=255)                                                                                                                              aaaTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGagag                                                                                                             <  1:21988/251‑1 (MQ=255)                                                                                                                               aaTGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAgc                                                                                                            >  1:156408/1‑251 (MQ=255)                                                                                                                                  gatgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAgcatg                                                                                                         <  1:131524/251‑1 (MQ=255)                                                                                                                                    tgaAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGAGGGTTTGGAACGGACTTTCCCTTCTGAATACAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAgcatgc                                                                                                        >  2:11888/1‑250 (MQ=255)                                                                                                                                          gCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGa                                                                                                 <  1:139702/251‑1 (MQ=255)                                                                                                                                              aTAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAAttt                                                                                             <  2:189932/251‑1 (MQ=255)                                                                                                                                                              ccAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAAGGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAAGTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGGCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACt                                                                             <  2:562244/251‑1 (MQ=255)                                                                                                                                                                 aaaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGACTTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGCATTTACTGCCAGTGAGTACTGAt                                                                          >  1:629802/1‑251 (MQ=255)                                                                                                                                                                 aaaaCGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTGTGCACAGGACTGGTGGGTGTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATTCTTCTGCTAATAATTTTCTCTGAGAGCATTCATTGTGAATTTACTTACAGTGAGTACTGAt                                                                          <  2:380585/251‑1 (MQ=255)                                                                                                                                                                  aaaCGCCACGTGATTCCGTGGCGTTATGCTTTTCTATCTGTAATCTTAATGACGCGCTGGGGATGTTAGGAAATTTCTCGGAATTTGCTGGCATGTTATGCAATATGCA‑TATCAATTGGTTAATTTTTGCACAGGCCTGGTGGGTTTGGAACGGAATTTCCCTTCTGACTAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGAtc                                                                         >  2:606745/1‑251 (MQ=255)                                                                                                                                                                         acgtgtttccgagtcgcttcacTTTTATACCTGAAATCTTCATGCAGCCCTGGCGATGTTCGGAATTTTCTTACAATTTGCCGGCCTGCTATGAAATTTGAA‑TATCAACTGGTAAATTTTTGCACAGGACTCGTGGGCTTGGAACGGTCTTTCACTTCTGAATAAAGGTATTAGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGACTTTACTGACAGTGAGTACTGAtctcctc                                                                    >  2:620420/22‑246 (MQ=255)                                                                                                                                                                         aCGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCCTGTTATGCAATTTGCA‑TATCAGATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCtcatc                                                                    >  1:514462/1‑249 (MQ=255)                                                                                                                                                                             gtTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTAtt                                                              <  1:331129/251‑1 (MQ=255)                                                                                                                                                                               ttACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTCCTGATCTCATCAGTAttt                                                             <  1:603856/250‑1 (MQ=255)                                                                                                                                                                                        gTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCAc                                                   <  1:116558/251‑1 (MQ=255)                                                                                                                                                                                            tGCTTTTATATCTGAAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGac                                                >  1:161886/1‑250 (MQ=255)                                                                                                                                                                                                tttATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCGCTGAGAGCATGCATTGTGAATTTACTGACAGGGAGTACTGATCTCATCAGTATTTTTTTATTCACGACActg                                            <  2:629802/250‑1 (MQ=255)                                                                                                                                                                                                  tataTCTTTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAAGTTTCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACActgtt                                          <  2:167483/250‑1 (MQ=255)                                                                                                                                                                                                   atatCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACActgtttc                                        >  1:73419/1‑251 (MQ=255)                                                                                                                                                                                                   atatCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACActgttt                                         <  1:103557/250‑1 (MQ=255)                                                                                                                                                                                                         gtaagctcagtGCCGCGGTGGCTATGTTAGGAAAGTTCTTGGAATTTGCTGGCATGTTATGCAATGTGCAGTA‑CAAATGGTTAATTTTTGCACAGGACTGGTGTGTTTGGAGCTGACTTTCCCTTCTGAATAAATGTCTTCGTGGTTATACTTCTGCTAATAAGTTTCTCGGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGAGCGCATCAGTATGTTGTTATTCGCGACCctgtttctgttgt                                  <  2:3579/241‑3 (MQ=255)                                                                                                                                                                                                                 aaTGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAACGGTCTTCGTGGTTATACTTCGGCTAATAATTTTCTCTGAGAGAATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGt                          >  1:514461/1‑251 (MQ=255)                                                                                                                                                                                                                 aaTGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCg                           <  2:156408/250‑1 (MQ=255)                                                                                                                                                                                                                  aTGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTTATTTTTGCACATGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTa                         <  1:145956/251‑1 (MQ=255)                                                                                                                                                                                                                   tGCCGCGCTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGCCTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTATTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTa                         >  1:258409/1‑250 (MQ=255)                                                                                                                                                                                                                      cgcgCTGCAGATGTTAGGTAAATTCTTGGAATTTGCTGGCCTGCTATGCAATTTGCA‑TAACACCTGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATCCTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGACTTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAggg                     >  2:543601/1‑250 (MQ=255)                                                                                                                                                                                                                       gcgcTGGCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTg                    <  1:113213/251‑1 (MQ=255)                                                                                                                                                                                                                            ggCGATGTTAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCa                 <  2:244578/249‑1 (MQ=255)                                                                                                                                                                                                                                    tAGGAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTGTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTAACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTc        <  2:453197/250‑1 (MQ=255)                                                                                                                                                                                                                                      ggAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTTAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTTTTTCTGTTGTGTGCCCGTAAGGTGGCACAAAAAGTCgg      <  1:186417/250‑1 (MQ=255)                                                                                                                                                                                                                                      ggAAAATTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGCTAATTTTTGCACAGGCCTGGTGGGTTTGGAAAGGAATTTCCCTTCTGAATAAAGGACTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGCATTTACTGACAGTGAGTACTGATCTCATCAGCATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAACAAGTCGGt     >  1:266701/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaaTTCTTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACTCTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTa  <  1:355793/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      GACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGTGTTTACGTGGCGTTTTGCTTTTATATCTGTAATCTTAATGCCGCGCTGGCGATGTTAGGAAAATTCCTGGAATTTGCTGGCATGTTATGCAATTTGCA‑TATCAAATGGTTAATTTTTGCACAGGACTGGTGGGTTTGGAACGGACTTTCCCTTCTGAATAAAGGTCTTCGTGGTTATACTTCTGCTAATAATTTTCTCTGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTA  >  NC_000913/3967532‑3968014

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