Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	624,596	0	A	T	29.4%	108.7 / 10.8	68	intergenic (‑86/‑289)	fepB/entC	ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (26/22);  new base T (0/20);  total (26/42)
Fisher's exact test for biased strand distribution p-value = 6.71e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCT  >  NC_000913/624457‑624733                                                                                                                                            |                                                                                                                                          aaGCTAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:830492/140‑1 (MQ=255) aaGCAACAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:129144/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTCATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:1442914/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:134815/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:630824/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:639060/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:640480/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:492748/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:841055/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:107913/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:1442851/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:1442839/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:134737/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:1316833/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:1101388/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:1314103/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:1481481/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:1500019/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  2:895122/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:1119242/140‑1 (MQ=255) aaGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:1076898/140‑1 (MQ=255)  aGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgag                                                                                                                                          <  2:382724/140‑1 (MQ=255)   gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTCGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGGTATTAATAAGGTTAAGGGCGTTATGACACATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgagt                                                                                                                                         <  2:294181/140‑2 (MQ=255)   gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGATATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                         <  2:1467800/140‑1 (MQ=255)               ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           >  2:255233/1‑126 (MQ=255)               ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           >  2:254491/1‑126 (MQ=255)               ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:254491/126‑1 (MQ=255)               ttAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGga                                                                                                                                           <  1:255233/126‑1 (MQ=255)                    tgatggggGTTGGGTTAGTGGGGGGCGTGTCTCAAAAATCAGCTTCCTTTTATTATTTTGTTTTAGGGCTTATTGAAAATTTCGAAAAAGCGCACATTCCGTCCCCTCGCCCCTTTGGGTTGTGGGTTTGGTTGTGGGGa                                                                                                                        <  2:187551/136‑1 (MQ=255)                    ctaaGGGGTTTGCGGTAGAGCGGGGCTACTCTCACAAATCAGCTACCTGTTATTAAAATGGTTACGGGCGAAATGAAAATTTCGACAAAGCCCAAATTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTATGGTTAGGGGa                                                                                                                        <  2:838624/138‑1 (MQ=25)                    aGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                        >  2:222633/1‑140 (MQ=255)                     gatgggCGTTGCGGTTGAGGGGGGCGTGTCTAAAAATTCTGCTTCCTTTTATTAATAAGGTTAAGGGCGTAATGACAAATTGGACAAAGCGCAAAATCGGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGTGGGGaa                                                                                                                       <  2:151257/137‑1 (MQ=25)                         ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                   <  2:540425/140‑1 (MQ=255)                         ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagc                                                                                                                   >  1:863666/1‑140 (MQ=255)                         ggCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCGGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGGTAGGGTGAGGGGAAcagc                                                                                                                   >  1:269751/1‑140 (MQ=255)                             ttGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccagc                                                                                                               >  2:1042286/1‑140 (MQ=255)                             ttGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAAcagccagc                                                                                                               >  2:884123/1‑140 (MQ=255)                              ggCGGTATAGCGGGGCTAGTCACAAAAATCAGCTTCCTGTTATTAATACGGTTAAGGGCTTAATGAAAAATTCGACAAAGCGCAAAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGTGGGGAACAGCCAGCa                                                                                                              <  1:1472024/139‑1 (MQ=255)                                     gagCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGc                                                                                                       >  2:520329/1‑140 (MQ=255)                                     gagCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGc                                                                                                       >  2:518519/1‑140 (MQ=255)                                      tgCGGGGCTTTTCTCAAAAATCAGCTTCTTTTTATTAATATGGTTATGGGCTTAATGACAAATTCGACAATGCGCAAAATCGGTCCCCTCGCCCGTTTGGGTTGTGGTTTCGGGTGTGGGGAACAGCCAGCACTGGTGCg                                                                                                      <  1:1504510/139‑1 (MQ=255)                                          gggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACa                                                                                                  >  2:1116958/1‑140 (MQ=255)                                          gggCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCACGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACa                                                                                                  >  1:1255958/1‑140 (MQ=255)                                             cgtgTCTCAAAAATCAGCTTCCTTTTATTAATAAGGTTAAGGGCGTAATGACAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTa                                                                                               <  2:1441680/137‑1 (MQ=38)                                             cGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTa                                                                                               >  1:1486495/1‑140 (MQ=255)                                                  ctcACAAATCAGCTTCCGGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                          <  2:527338/140‑1 (MQ=255)                                                  atcACAAATCTGCTTCCTTTTATTAATATGTTTTTGGGCGTAATGAAAATTTCGACAATGCGCAAATTCCTTCCCCTCGCCCTTTTGGGTAGTGGTTTTGGGTTTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                          <  1:1428363/139‑1 (MQ=255)                                                  atcAAAAATCCGATTCTTTTTATTATTATGGTTATGGTCTTATTTACAATTTCGAATAAGCGCACAATCCGTCCCCTCGCCCTTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                          <  1:1284946/139‑1 (MQ=255)                                                        aaTCCGCTTCATTTTATTAATACGTTTATGGGCGTAATGAAAAATACGACAAAGCGCAAAATCCGTCCCCACGCCCCTTTGGGGTGAGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAcccc                                                                                    <  1:643967/140‑1 (MQ=37)                                                          tCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAg                                                                                  >  2:220397/1‑140 (MQ=255)                                                            tgCTTCCTTTTTTTAATATGTTAATGGGCTTAATTACAATTCCGAAAACGCCAACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGcc                                                                                <  1:432818/139‑1 (MQ=16)                                                              gttCCTGTTATTAATAAGTTTAAGGGCGTAATGAAAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  2:802851/139‑1 (MQ=37)                                                              cttcttTTTATTAATACGGTTATGGGCGTAATTACAAATTCCACAACGCGCACAATCCGTCCCCAGGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  2:822116/140‑1 (MQ=25)                                                              cTTCGTTTTTTTAATAAGTTTATGGGCTTAATGACAATTTCTAAAACGCGCACAATCCGTCCCCCCGCCCGTTGGGGTTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  1:1181406/140‑1 (MQ=255)                                                              cTTCCTTTTTTTAATATGGTAATGGGCTTAAGCACAAATCCTCCAATGCCCACATCCCGCCCCCCCGCCCCTTGGGGGATCGGTTTCGGGTGCGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  1:325482/140‑1 (MQ=255)                                                              cTTCCTTTTATTAATACGTTTAAGGGGTTATTGACAATTTCTCCAAAGCGAACAACCCTACCCCCCGCCCCTTTGGGGAGAGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  1:519885/140‑1 (MQ=18)                                                              cTTCCTGTTTTTAATATGGTTCCGGGCGTAATGACAAATTCGACAAAGCGCTCAAACCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  2:135477/140‑1 (MQ=25)                                                              cTTCCTGTTATTAATACGGTTACGGGCGTAATGACAATTTCGAAAAAGCCCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              <  1:739840/140‑1 (MQ=25)                                                              cTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                              >  1:1030286/1‑140 (MQ=255)                                                              cTTCCTGGGATTAAGAAGGGTCAGGGCGTAAGGCCACATTCGACAAGGCGCACAATCCGTCCCCTCGCCCCTGTGGGGAGTGGGGTAGGGGGCGGGGACCCGCCAGCAGTGGTGCGAAGAGTAACCCTCACCCCAGCCCt                                                                              >  2:457626/1‑140 (MQ=255)                                                                  cTGTTATTAATAAGGTTAAGGGCGAAATGACAACTGCGACAAAGCGCACAATCCGGCCCCTCGCCCCTTTGGGGTGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCAcc                                                                          <  2:271788/140‑1 (MQ=37)                                                                       ataaatatGGTTTCGGGCTTAATCACATATTCGACAACGGGCACATCCCGTCCCCTCGCCCCTTGGGGGTGCGGGTTGGGGTGTGGGGAACCGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGa                                                                     <  1:119263/137‑1 (MQ=255)                                                                             ctggatAAGGGCGAATATACAAATTCGACAATGAGAACTAGCCGTCCCATCGGCCCATTGGGGTGAGGGTTAGGTTGTGGGTAACAGGCATAACTGTTTCGAACATAACCCCTCACCCCAGCCCTCACCCTGTAAGGgag                                                               <  2:614799/135‑1 (MQ=255)                                                                             aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                               >  2:230780/1‑140 (MQ=255)                                                                             aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                               >  2:229539/1‑140 (MQ=255)                                                                             aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                               >  2:739097/1‑140 (MQ=255)                                                                             aaGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGCCCCCACGCCCCATTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgag                                                               <  2:616141/140‑1 (MQ=255)                                                                             aaGGTAAAGGGCGTAATGACAAATCCCAAAAACCGCACAACCCGTCCCCTCTCACCGTTGGGGAGTTGGTTACCGTGTGCGGAACAGCCAGCACTGGTGCGATCCCTCCCCCTCACCCCAGCCCACACCCTGGAAGGgag                                                               <  2:616096/140‑1 (MQ=255)                                                                              tgtttACGGGCTTATTTAAAAATACCACAAAGCGAAAAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGAGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                              <  1:942062/137‑1 (MQ=18)                                                                              tggTTACGGGCTTAATTAAAAATACCAAAAAGCGCAAAACCCGCCCCCTCCCCCATTGGGGGTGTGGTTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGgaga                                                              <  1:1191532/139‑1 (MQ=255)                                                                                   acgggCGTACGGCCAAATTCGACAAAGCGCCCACGCCGGCCCCGCGCCCCTTGGGGGTGAGGGTTCGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggggg                                                         <  1:1444027/138‑1 (MQ=255)                                                                                        cGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                    <  1:167727/140‑1 (MQ=255)                                                                                          tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                  >  1:704947/1‑140 (MQ=255)                                                                                          tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                  >  1:610811/1‑140 (MQ=255)                                                                                          tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                  >  1:610839/1‑140 (MQ=255)                                                                                             tGACAAATTCGACAACGCGCACAATCCTTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                               <  2:265766/140‑1 (MQ=255)                                                                                             tGAAAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCGTTGGGGGTGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                               <  1:1161989/140‑1 (MQ=255)                                                                                               caaaaTTCCACAACGCCCAAAAACCGTCCCCACGCCCCTTTGGGGTGAGGGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                             <  1:874986/138‑1 (MQ=38)                                                                                               aCAAATTCTACAAAGCGAAAAATCCGTCCCCTCGCCCATTTGGGTTGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                             <  1:966387/140‑1 (MQ=255)                                                                                               aCAAATTCTAAAAAGCGCACAATCCGTCCCCTCGCCCATTTGGGGAGAGGGTTCGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                             <  1:1297926/140‑1 (MQ=255)                                                                                               aCAAATTCGACAAAGCGCAAAATCCTACCCCTCCCCCCTCTGGGTAGTGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                             <  1:328412/140‑1 (MQ=255)                                                                                                 aatttCGACAACGCGCCCAATCGGGCCCCGCGCCCGTTGGGGGTGTGGGTTGGGGTGGGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAgg                                           <  1:937676/137‑1 (MQ=255)                                                                                                  aaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGa                                          >  2:1495271/1‑140 (MQ=255)                                                                                                    ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGGAGAGGGTTAGGGTTAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                        <  1:556298/140‑1 (MQ=255)                                                                                                         cAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                   >  1:193409/1‑140 (MQ=255)                                                                                                         cAAACCGCAACACCCGGCCCCCCGCCCCTTTGGGGATCGGGTTTTGGTGGGGGTAACAGCCACCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcaca                                   <  1:1121860/140‑1 (MQ=255)                                                                                                                        gTCCCCTCGCCCCTTTGGGGAGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCTTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                    <  1:771652/140‑1 (MQ=255)                                                                                                                        gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                    >  1:288040/1‑140 (MQ=255)                                                                                                                        gTCACCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                    <  2:1093929/140‑1 (MQ=255)                                                                                                                           cccTCGCCCCTTTGGGGAGAGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCACACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                 <  2:112475/140‑1 (MQ=255)                                                                                                                           cccTCGCCCCTTGGGTGTGTGGGTTAGGTTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                 <  1:987038/140‑1 (MQ=255)                                                                                                                               cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa             >  2:15323/1‑140 (MQ=255)                                                                                                                               cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa             >  2:979386/1‑140 (MQ=255)                                                                                                                                gCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaa            >  2:1306250/1‑140 (MQ=255)                                                                                                                                gCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaa            >  2:1379129/1‑140 (MQ=255)                                                                                                                                          ggAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCCCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTtct  <  1:691132/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                          AAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCT  >  NC_000913/624457‑624733

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