Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,543,047	T→G	intergenic (+414/‑68)	fimE → / → fimA	regulator for fimA/type 1 fimbriae major subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,543,047	0	T	G	97.9%	151.0 / ‑6.4	47	intergenic (+414/‑68)	fimE/fimA	regulator for fimA/type 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (30/16);  total (31/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATCAAACCGTTCAGTTAGGACAGGTTCG  >  NC_000913/4542801‑4543296                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          aattTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGaca                                                                                                                                                                                                                                                        >  1:394160/1‑250 (MQ=255)   tttCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGacaaa                                                                                                                                                                                                                                                      >  2:339106/1‑250 (MQ=255)        caTCACCTCCGCTATATGTAAAGCTAAAGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGACCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGacaaaacaa                                                                                                                                                                                                                                                  >  2:540487/1‑249 (MQ=255)                       atGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATACAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GCGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACCGAACGACTg                                                                                                                                                                                                                                 >  1:588364/1‑251 (MQ=255)                                       tCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAg                                                                                                                                                                                                                  <  1:624223/250‑1 (MQ=255)                                                 gACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGttttt                                                                                                                                                                                                       <  1:121206/251‑1 (MQ=255)                                                   cGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCGCTATGAGTCAAAATGGCCCCAAAGGTGTCATCGTTTGGGGGAAAACTGTGCA‑GTTTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTg                                                                                                                                                                                                     <  1:377766/251‑1 (MQ=255)                                                            cTCATTCTTCTCTCCAAAAACCAGCTCATGCAATATAAACATCTATAAATAACGATAACAATAGAATATTAAGCAAACAAATAAAATGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCACATGTTTCATCTTTTGGTGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTTACAAAACAAATTGTACAGAATGACTGCCCATGTCGCTTTAGTAATAGTTTTTTgaagtgaa                                                                                                                                                                                              >  2:14483/1‑244 (MQ=255)                                                                    tctctcCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAgcagca                                                                                                                                                                                       <  1:333011/248‑1 (MQ=255)                                                                         ccAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTGGAAAGGAAAGCAGCATGAAAAt                                                                                                                                                                                <  2:412717/250‑1 (MQ=255)                                                                          cAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTa                                                                                                                                                                              <  1:348008/251‑1 (MQ=255)                                                                                       aTGCAATATAACCATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGCAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAt                                                                                                                                                                 >  1:140738/1‑251 (MQ=255)                                                                                       aTGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAAt                                                                                                                                                                 >  1:260954/1‑251 (MQ=255)                                                                                        tGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATc                                                                                                                                                                >  1:281150/1‑251 (MQ=255)                                                                                            atatAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCAGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGCTTTAGAAATAGTTTTTTGACAGGACAGCAGCATGAAAATTAAAACTCTGGCAATCgttg                                                                                                                                                            >  1:249154/1‑251 (MQ=255)                                                                                                aaaCATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCGCCAAATGTTTCATCTATTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACTAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATATTGTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTc                                                                                                                                                         <  1:408391/250‑1 (MQ=255)                                                                                                        ataaataaAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGctc                                                                                                                                                <  2:112275/251‑1 (MQ=255)                                                                                                                        aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGt                                                                                                                                                                                                                          <  1:123626/161‑1 (MQ=255)                                                                                                                        aTAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGt                                                                                                                                                                                                                          >  2:123626/1‑161 (MQ=255)                                                                                                                                 tAAGCCACCAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCACCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGccct                                                                                                                       >  2:35900/1‑248 (MQ=255)                                                                                                                                     acaacaaATAAACGGAAAAAATTTGGCCGGGATGCTTTCCTCTATGAGTCAAAAGGGCCCCAAATGTGTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGGCGATGTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTTGGCTCTGGCCCTCAGTTCTTCAGCGGCTCTGGct                                                                                                                   <  2:80475/250‑2 (MQ=255)                                                                                                                                         caaATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAGCTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAACTAGTTTTTTGAAAGGAAAGCAGCATGACAATTAAATCTCTGGCAATCGTTGTTCTGTCGGCGCTGTCCCTCAGTTCTACAGCGGCTATGGCAGCTg                                                                                                               >  2:355351/1‑251 (MQ=255)                                                                                                                                            aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACa                                                                                                                                                                                                                                          >  1:354084/1‑125 (MQ=255)                                                                                                                                            aTAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACa                                                                                                                                                                                                                                          <  2:354084/125‑1 (MQ=255)                                                                                                                                                    aaaaaGTTTGTCTGAGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGTGACACATGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGACCGACTGCCCATGTCGATTTAGAAATAGTTTTATGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGGTGTTCTGTCGGCTCTGTCTCTCAGTTCTACAGCGTCTCTGGCCGCTTCCACGACGGtt                                                                                                    >  2:243056/1‑251 (MQ=255)                                                                                                                                                        aGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAAtg                                                                                                >  1:109266/1‑251 (MQ=255)                                                                                                                                                          tttGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAAtggt                                                                                              >  1:25715/1‑251 (MQ=255)                                                                                                                                                            tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGcc                                                                                                                   >  1:543378/1‑228 (MQ=255)                                                                                                                                                            tGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGGTCTACAGCGGCTCTGGcc                                                                                                                   <  2:543378/228‑1 (MQ=255)                                                                                                                                                                             tctATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTaa                                                                             >  1:20897/1‑249 (MQ=255)                                                                                                                                                                              ctATGAGTCACACTGGCCCCAAATGTTTCATCTTTTGGCGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGGCGGTTAATGGTGGGACCGTTCACTTTAAAgg                                                                          >  2:237755/1‑251 (MQ=255)                                                                                                                                                                                aTGAGTCAAAATGGCCCCAAATGTTTCATCTTTGGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCCCAGTTCTACAGCGGCTCTGGCCGCTGCCACTACGGTTAATGGTGGGACCGTTCACTTTAAAggg                                                                         <  2:459102/250‑1 (MQ=255)                                                                                                                                                                                        aaaTGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACCCGTTGCCCAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTAAGTTCTACAGCGGCTCTGGCCGCTGCACCGACGTTTAATGGTGGGACCGTTCACTTTAAAGGGGAAgttgtt                                                                >  2:250964/1‑251 (MQ=255)                                                                                                                                                                                              ccccAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAAcgcc                                                          >  1:245644/1‑251 (MQ=255)                                                                                                                                                                                                 cAAATGTTTCATCTATTGGGGGAAAACTGTGCAAGT‑TTGGCATTCAAACTCGTTGTCAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAAATAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCATTTCTTCAGCGGCTCTGGCCTCTGCCACGACTGTTAATGTTGGTACCGTTCGCTGTAAATGGGAAGTTTTTAACGCCGCt                                                       <  2:40437/251‑1 (MQ=255)                                                                                                                                                                                                    atgattCATCTTTTGGGGGAAAAATGTGCA‑GTGTTGGCAGTCAAACTCGTTGACTAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGATCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTgc                                                    >  2:392012/5‑251 (MQ=255)                                                                                                                                                                                                      gTTTCCTCTATTGGGGGAAACCTGTGCA‑GTGTTGGCCGTCTCATTCGTTGACACAACAAAAAGTACAGACCGACTGCCCCTGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTAGGCTCTGTCCCTCAGTTCTACCGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTACACGCCGCTTgcgc                                                  >  2:367497/1‑251 (MQ=255)                                                                                                                                                                                                          cATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGtt                                              >  2:612519/1‑251 (MQ=255)                                                                                                                                                                                                          cATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGtt                                              >  1:372899/1‑251 (MQ=255)                                                                                                                                                                                                                 tGGGGGAAAACTGTGCA‑TTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTGAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCGCAGTGCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTGAACGCCTCTTGCGCAGTTGATGCa                                        <  2:281150/250‑1 (MQ=255)                                                                                                                                                                                                                            tgtgCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCCCTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCCGGCTCTGTTGAt                            >  1:342023/1‑251 (MQ=255)                                                                                                                                                                                                                             gtgCA‑GTGTTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTATGTCGGCTCTGTCCATCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTACTGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTgctc                           >  2:205570/1‑248 (MQ=255)                                                                                                                                                                                                                                   gtgtTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATCAAAcc                      <  2:178725/251‑1 (MQ=255)                                                                                                                                                                                                                                    tgtTGGCAGTCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGGGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCCCTGTTGATCAAACCg                     >  2:455662/1‑251 (MQ=255)                                                                                                                                                                                                                                           aGTCAAACTCGCTGACAAAACAAAGTGTACAGACCGACTGCCCATGTCGATTTAGAAATCGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTATGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATCAAACCGTTCAGtt              >  2:606837/1‑251 (MQ=255)                                                                                                                                                                                                                                             tCAAACTCGTTGACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTCAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATCAAACCGTTCAGtt              >  2:240922/1‑249 (MQ=255)                                                                                                                                                                                                                                               aaaCTCGTTGACAAAACAAAGAGTACAGACAGACTGCCCATGCCGATTTAGACATAGCATTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTCATCAAACAGTTCAGTTAgg           >  2:120554/1‑250 (MQ=255)                                                                                                                                                                                                                                                       ttACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATCCAACCGTTCAGTTAGGACAGGTTCg  >  2:174805/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          AATTTCACATCACCTCCGCTATATGTAAAGCTAACGTTTCTGTGGCTCGACGCATCTTCCTCATTCTTCTCTCCAAAAACCACCTCATGCAATATAAACATCTATAAATAAAGATAACAATAGAATATTAAGCCAACAAATAAACTGAAAAAGTTTGTCCGCGATGCTTTCCTCTATGAGTCAAAATGGCCCCAAATGTTTCATCTTTTGGGGGAAAACTGTGCA‑GTGTTGGCAGTCAAACTCGTTTACAAAACAAAGTGTACAGAACGACTGCCCATGTCGATTTAGAAATAGTTTTTTGAAAGGAAAGCAGCATGAAAATTAAAACTCTGGCAATCGTTGTTCTGTCGGCTCTGTCCCTCAGTTCTACAGCGGCTCTGGCCGCTGCCACGACGGTTAATGGTGGGACCGTTCACTTTAAAGGGGAAGTTGTTAACGCCGCTTGCGCAGTTGATGCAGGCTCTGTTGATCAAACCGTTCAGTTAGGACAGGTTCG  >  NC_000913/4542801‑4543296

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