Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,540,114	A→T	intergenic (‑613/‑843)	nanC ← / → fimB	N‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,540,114	0	A	T	100.0%	120.7 / NA	40	intergenic (‑613/‑843)	nanC/fimB	N‑acetylneuraminate outer membrane channel/Type 1 fimbriae regulatory protein FimB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/22);  total (18/22)

AGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGA  >  NC_000913/4539869‑4540349                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGt                                                                                                                                                                                                                                        <  1:458403/251‑1 (MQ=255)           gCCATAGCTTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAACACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAgg                                                                                                                                                                                                                              >  2:385777/1‑251 (MQ=255)           gCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAg                                                                                                                                                                                                                               <  1:12415/250‑1 (MQ=255)            ccATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAggg                                                                                                                                                                                                                             <  2:245371/251‑1 (MQ=255)            ccATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAggg                                                                                                                                                                                                                             >  1:571887/1‑251 (MQ=255)                aGATTGCGAGATTCCCTTGCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAATAACTAAAGATCTTAGTTTAACTATGTGTTTTATAAAGAATTGGGTTAGAGTCTAAACAAGGGGAg                                                                                                                                                                                                                          <  2:369833/250‑1 (MQ=255)                     gCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGc                                                                                                                                                                                                                    <  2:143740/251‑1 (MQ=255)                                     cATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGt                                                                                                                                                                                                       <  1:524248/248‑1 (MQ=255)                                       ttGACTATCTCGGTGCACATTTATACCTGTTATACCAGATCCAAACTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGACCACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCt                                                                                                                                                                                                  >  2:408425/1‑251 (MQ=255)                                             aTCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTg                                                                                                                                                                                             <  1:320467/250‑1 (MQ=255)                                                      cACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGTGAGCTGTGCAAGCTAACTTAGTGAGCTTGGTGaaatccagtg                                                                                                                                                                                   <  2:95247/251‑8 (MQ=255)                                                             ataCCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTAccc                                                                                                                                                                            <  1:385777/251‑1 (MQ=255)                                                                 cTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGcc                                                                                                                                                                         >  1:465358/1‑250 (MQ=255)                                                                                    aTCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGATCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATtct                                                                                                                                                     <  1:11700/251‑1 (MQ=255)                                                                                      cACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCTATCTTTTTCTAAGAGGAAGATGCCATGTGAATCCATACAAACACTTGCGGTGGTCTTCAAAAACTAAAGAGCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTTCAAGCTAACTCAGTGAGCTGGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtca                                                                                                                                                   <  1:556182/251‑1 (MQ=255)                                                                                         gaaaTCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGAGGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCtcatca                                                                                                                                                <  2:180524/249‑1 (MQ=255)                                                                                           aaTCCATACAACAAAACCAGATTTGCAATTCTTGTCACAAAAAATGCCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACAATTGCGGTGGTCTTCAAAACAGAAAGTGCTAATTTTAACTATTTGTTTTAAAAATAATTTGTTTAGAGTCAAAACAAGGGGAGCTTTGCATGCTAACTCAGTGAGCTTGGTGAAAACCAGAGTTCTCCCGCCATCAGTCTGAGCATAAGTCTcatcatt                                                                                                                                              <  2:589416/251‑2 (MQ=255)                                                                                                  acaacaAAACCAGATTTGCAATTCGGGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATg                                                                                                                                       <  1:75306/251‑1 (MQ=255)                                                                                                    aacaaAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGtt                                                                                                                                     >  2:290681/1‑251 (MQ=255)                                                                                                       aaaaCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAATAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGAGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGtttac                                                                                                                                  >  2:358312/1‑249 (MQ=255)                                                                                                          aCCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTgg                                                                                                                               >  1:210613/1‑251 (MQ=255)                                                                                                           ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCATACGAACACTTGCGGTGGTCTCCAAAAACTAAAGATCTTAGGTTACCTAGTTGTTTTATAAAGAAGTTGTTTATAGTCTAAACAAGGGGAGCATGGCAAGCTaa                                                                                                                                                                                                             <  2:264663/172‑1 (MQ=255)                                                                                                           ccAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCCAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTaa                                                                                                                                                                                                             >  1:264663/1‑172 (MQ=255)                                                                                                                  tGCAAGTCGTGGCTCAAAATAGGTGGATCTTTTACAAAGAGAAAGATGCCATGTGAAGCCAGATGAAAACTTGCGGTGGTCATAAAAAACTAAAGATGTTAGTTGAACTATTTGTTTAATAAATAATTTGTTTAAAGTCTAAACAAGGGGAGGTTGCCAATCTAACTCAGTGATCTTGGTGAAAATCAGTGTTGTCCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTgg                                                                                                                        <  2:371373/250‑1 (MQ=255)                                                                                                                                 aaaaTATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACt                                                                                                        >  1:70530/1‑251 (MQ=255)                                                                                                                                   aaTATGTCGATTTTTTTCTAGGAGGAAGATGCCATGTGCAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATAATGCAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTccattt                                                                                                      >  2:217044/1‑246 (MQ=255)                                                                                                                                                ttttCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTAAAAAACTAAAGATCTTAGTTTAACTATTTTTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTGGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATc                                                                                           <  1:87622/249‑1 (MQ=255)                                                                                                                                                ttttCTAAGAGGAAGATGCCATGTGAAGCCAGACAAACACTTGCGGTGTTCTTCAAAAACTAAAGATCTTAGTTTAACTATTGGTTTTATAAATAATTTGTTTAGAGTCAAAACAAGGGTAGCGTTGCAAGCGAACTCAGTGAGCTTGGTGAAAATCAGGGTTTCCCCGCCATCATGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGATTGTAACTGGTCACTTCTGAAGTCg                                                                                              <  2:558776/246‑1 (MQ=255)                                                                                                                                                        gagGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCt                                                                                 >  2:396330/1‑251 (MQ=255)                                                                                                                                                         agGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCtt                                                                                <  1:216578/251‑1 (MQ=255)                                                                                                                                                             aGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTTGTCACTTCTGAAGTCGATCTGGAGAGGCttgttt                                                                            >  1:480385/1‑250 (MQ=255)                                                                                                                                                                             cAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAg                                                            <  2:229232/251‑1 (MQ=255)                                                                                                                                                                             cAGACGAACACTTGCGGTGGTCTTAAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAg                                                            <  2:36510/251‑1 (MQ=255)                                                                                                                                                                                aCGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGgat                                                         <  1:548317/251‑1 (MQ=255)                                                                                                                                                                                                               aGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACCAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGCAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCTGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGc                          >  2:575140/1‑251 (MQ=255)                                                                                                                                                                                                               aGATCTTAGTTTAACTATTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGAAGCTTTGCATGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTAACTTCTGAAGTCGATCTGGAGAGGTTTGTTGATGTTTGTGGTTTCAGGATGATGTTTCACTTTGGTTGTTTGCCGttatgc                          >  2:160032/1‑246 (MQ=255)                                                                                                                                                                                                                               aTTTGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCtgtg          >  1:540181/1‑251 (MQ=255)                                                                                                                                                                                                                                tttGTTTTATAAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCtgtg          >  1:64919/1‑250 (MQ=255)                                                                                                                                                                                                                                       tataAATAATTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATAAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGGTGGTGTTTTCAGGATGATGTTTCACTTGGTTTGTTTGCCGTATCGCCCGGCGAATGGCCGTGATTAAgg   >  1:267192/1‑250 (MQ=255)                                                                                                                                                                                                                                        ataaataaTTTGTTTAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGa  >  1:143304/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            AGTGATGTTTGCCATAGATTGCGAGATGCCCTGTCTCATTGACTATCTTGGTGCACCTTTATACCTGTTATACCAGATCAAAAATCACGCAATCCATACAACAAAACCAGATTTGCAATTCGTGTCACAAAATATGTCGATCTTTTTCTAAGAGGAAGATGCCATGTGAAGCCAGACGAACACTTGCGGTGGTCTTCAAAAACTAAAGATCTTAGTTTAACTATTTGTTTTATAAATAATTTATTAAGAGTCTAAACAAGGGGAGCTTTGCAAGCTAACTCAGTGAGCTTGGTGAAAATCAGTGTTTACCCGCCATCAGGCTGAGCATAATTCTCATCATGAAATATGTTTCCTGGTTTGTGGCTTGTAACTGGTCACTTCTGAAGTCGATCTGGAGAGGCTTGTTGATGTTGGTGTTTTCAGGATGATGTTTCACTTAGTTTGTTTGCCGTATCGCCCGGCGAATGGCTGTGATTGAGGA  >  NC_000913/4539869‑4540349

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: ‑