Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	624,598	0	A	T	26.0%	108.9 / 12.7	74	intergenic (‑88/‑287)	fepB/entC	ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base A (19/35);  new base T (0/19);  total (19/55)
Fisher's exact test for biased strand distribution p-value = 1.71e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624459‑624737                                                                                                                                            |                                                                                                                                            gCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGgaga                                                                                                                                             <  2:533746/140‑1 (MQ=255)   aataGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTTTTAATAAGGTTAAGGGCGAAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAgg                                                                                                                                           >  2:1558849/4‑140 (MQ=255)      aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                        <  1:1069623/140‑1 (MQ=255)      aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                        <  1:1112770/140‑1 (MQ=255)      aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                        <  2:846890/140‑1 (MQ=255)      aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                        <  1:134338/140‑1 (MQ=255)      aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGtt                                                                                                                                        <  1:362838/140‑1 (MQ=255)          cTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTTggg                                                                                                                                    <  2:1441221/140‑1 (MQ=255)           ttttAATTGACGGGCGTTGCGGTTGAGGGGGGCTTGTCTCAAAAAACAGATTCATGTTTTTAATAAGGTTAAGGGCGTATTGAAAAATTCCAAAAAGCGAAAAAACGGTCCCCTCCCCCGTTTGGGTTGTGGGTTTggtt                                                                                                                                   <  2:599661/138‑3 (MQ=255)            cttAATAGAAGGGCGTTGCGGTCCAGCGGGGCGAGTCACACAACACAGCTCCCTGTTATTAATCAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTg                                                                                                                                  <  2:883243/139‑1 (MQ=38)             ttattaGATGGGGGTTGGGGTTGTGGGGGGGGTGTCTCAAAAATCAGATTCTTGTTTTTAATATGTTTAAGGGCTTATTTAAAAATTCGACAACGCGCACAATCCGTCCCCTCGCCCTTTTGGGTTGTGGGTTTGGGTgt                                                                                                                                 <  2:258112/140‑2 (MQ=255)                  aGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGTGGa                                                                                                                            >  2:1524907/1‑140 (MQ=255)                  aGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                            >  2:984457/1‑140 (MQ=255)                   gAAGGGCGTTGCGGTAGAGGGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                           <  2:1518792/140‑1 (MQ=255)                   gAAGGGCGGTGCGGTAGTCGGGGGCGATTCTCACAAACCAGCTTCCTTTTATTAATATGGTTAAGGGCGTACTGACAAATACGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGAGGGTTTGGGTGAGGGGaa                                                                                                                           <  2:1518840/140‑1 (MQ=255)                            tGCGGTAGAGCGGGGCGAGTCTCACCAATCCGCTTCCTGTTATTAATAAGGTTAAGGGCGTACTGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTAAGGGGAGAGGGTTAGCGAGAGGGGACCAGCCAGCa                                                                                                                  >  2:1545668/1‑140 (MQ=255)                                   ttggggggggAGTATAAAAAATCAGCTACCTTTTTTTAATATGGTTAAGGGCTTATTGAAAATTTCGAAAAAGCGAACAATCCGTCCCCTCGCCCCTTTGGGTAGTGGGTTTGGGTGCGGGGAACAGCCAGCACTGGTGc                                                                                                           <  1:1392570/136‑1 (MQ=19)                                   gagCGGGGCGAGTCACACAAATCAGCTTCCTTTTATTAATAAGGTTATGGGCTTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGc                                                                                                           <  2:178760/140‑1 (MQ=255)                                          gCGTGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACAtt                                                                                                    <  2:704808/140‑1 (MQ=255)                                           cGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTa                                                                                                   >  2:176130/1‑140 (MQ=255)                                                ctcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:1548178/140‑1 (MQ=255)                                                ctcACAAAACAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:1548690/140‑1 (MQ=255)                                                ctaacaATTCAGCTTCCTGTTTTTATTATGTTTATGGGCTTAATGAAATATTCGAAAATGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:532030/137‑1 (MQ=21)                                                ctaaaaaaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:620157/135‑1 (MQ=255)                                                cacacaAAACAGCTTCCTGTAATTAATACGGTTATGGGCGTAATGACAAAGTCCACAACGCGCACAAGCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  1:1456633/138‑1 (MQ=25)                                                cacaaaaaTCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATTAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGTTGAGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:982050/135‑1 (MQ=38)                                                atcACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCt                                                                                              <  2:1490828/139‑1 (MQ=255)                                                   acaAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAc                                                                                           >  2:956370/1‑140 (MQ=255)                                                      aaccAGCTCCCTGTAATTAAAATGGTTATGGGCTAAATCACAATTACGACAAACCGCAAACCCCGACCCCCCCCCCCTTTGGGTCGTGGTTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCAcccc                                                                                        <  1:205060/137‑1 (MQ=255)                                                            cttcttGTTATTAATACGTTTACGGGCTTAATGAAAAATTCGACAAAGCGCAAAATCGGTCCCCTCGCCCCTTTGGGGTGTGGTTTTGGGTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  1:513933/140‑1 (MQ=25)                                                            cTTCGTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:512052/140‑1 (MQ=255)                                                            cTTCCTGTTTTTAATAAGGTTATGGGCGTAATGACAAATTCCACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:512279/140‑1 (MQ=255)                                                            cTTCCTGTTATTAATACGGTTATGGGCGTAATGACAAATTCGACAACGCGCACAATCCGTCCCCTCGCCCCTTTGGGGTGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:370553/140‑1 (MQ=255)                                                            cTTCCTGTTATTAATAAGGTTACGGGCGTAATGACAAATCCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTGGGGGAGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  2:24520/140‑1 (MQ=255)                                                            cTTCATGTTATTAATAAGGTTATGGGCGTAATGAAAATTTCGACAAAGCGCACAATCCGTCCCCTCGCCCATTTGGGGTGTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCt                                                                                  <  1:1608636/140‑1 (MQ=255)                                                              tcgtGTTATTAATTAGGTAAAGGGCGTAATGACAAATTCGAAAAAGCGCACAATCCGGCCCCTCGCCCCTTTGGGGAGTGGGTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  2:1013887/137‑1 (MQ=37)                                                              tCCTGTTATTAATCCGGTTAAGGGCGTAATTACAAATTCGACAAAGCGCACAAACCGTCCCAACGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCa                                                                                <  2:1013779/140‑1 (MQ=25)                                                                      ttaATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGaa                                                                        >  1:886124/1‑140 (MQ=255)                                                                          tatggTTAAGGGCGTAATGACAAATTCGACAACGCGCACAATCCGACCCCTCGCCCCTTTGGGTAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGga                                                                    <  2:1099706/137‑1 (MQ=255)                                                                                 aaGGGCGTAATGACAAATTCGACAATGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAggggg                                                             <  2:1468833/140‑1 (MQ=255)                                                                                      cGTAATGACAAATTCGACAAAGCGCACAACCCGACCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGaa                                                        <  1:701199/140‑1 (MQ=255)                                                                                        tAATGACAAATTCGACAAAGCGCACATTCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGTGCATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  2:154290/1‑140 (MQ=255)                                                                                        tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:376865/1‑140 (MQ=255)                                                                                        tAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACg                                                      >  1:511470/1‑140 (MQ=255)                                                                                           ttacacATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1617626/138‑1 (MQ=255)                                                                                           tccaaaaTTCCAAAACGCGCACAATCCGCCCCCCCGCCCCTTGGGGGAGTGGGTTAGGTTGCGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1352189/136‑1 (MQ=255)                                                                                           tcaCAAATTCGACAACGCCCACAATCCGTCCCCACGCCCATTTGGGGCGAGGGTTTGGGGGTGGGCAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1069571/138‑1 (MQ=255)                                                                                           tGACAAATTCGACAACGCGCACAATCCGCCCCCTCGCCCCTTTGGGGAGAGGGTTCGGGTGAGGGGACCCGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:1352127/140‑1 (MQ=25)                                                                                           tGAAAAATTGCACAATGCGAACAATCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGgcg                                                   <  1:468459/140‑1 (MQ=21)                                                                                             aaaatttCGACAACGCGCACAACCCGTCCCCTCGCCCCTTTGGGGTGTGGGTTTGGGGGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                 <  1:848776/135‑1 (MQ=25)                                                                                             aaaaaTTCGACAAAGCGCACAATCCGTCCCCTCGCCCATTTGGGGTGAGGGTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCa                                                 <  1:848724/138‑1 (MQ=255)                                                                                               aatttCGAATAAGCGCCCATTCCGTCCCCCCCCCCATTGGGGGTGCTGGTTTGGGTTAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCGCAGCCCTCACCCTGGAAGGGAGAGGGGGCTGAACGGCGCAgg                                               <  1:887752/137‑1 (MQ=255)                                                                                                 ctaccacacAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTCGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGAc                                             <  1:1607374/135‑1 (MQ=255)                                                                                                  ttCGACAATGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  2:805235/140‑1 (MQ=255)                                                                                                  ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  2:1460655/140‑1 (MQ=255)                                                                                                  ttCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACa                                            <  2:642476/140‑1 (MQ=255)                                                                                                    cGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATc                                          >  1:1582738/1‑140 (MQ=255)                                                                                                        aaaGCGCACACTCCTTCCCCTCGCCCCTTTGGGTATAGGGTTTGGGTGAGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                      <  1:282957/140‑1 (MQ=255)                                                                                                           gcgcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                                   >  2:1153302/1‑140 (MQ=255)                                                                                                           gcgcACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTgc                                   >  2:1150442/1‑140 (MQ=255)                                                                                                                      gcccccTCCCCCTTTTGGGTAGTGGTTTAGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                        <  1:150564/138‑1 (MQ=255)                                                                                                                      gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                        <  2:1528804/140‑1 (MQ=255)                                                                                                                      gTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGaa                        >  1:1339326/1‑140 (MQ=255)                                                                                                                         ccccccccctTTTGGGGAGTGGTTCAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     <  1:858299/134‑1 (MQ=255)                                                                                                                         ccccCGCCCCTTTGGGTATTGGGTTTGGTTTTGGGTAAGAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     <  1:859921/140‑1 (MQ=255)                                                                                                                         cccTCGCCCCTTTGGGGAGTGGTTTTGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     <  1:994043/140‑1 (MQ=255)                                                                                                                         cccTCGCCCCTTGGGTGAGTGTGTTTGGGTGTGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                     <  1:994190/140‑1 (MQ=255)                                                                                                                           cTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCAt                   >  1:1209588/1‑140 (MQ=255)                                                                                                                            tCGCCCCTTTGGGGCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATc                  <  2:1310556/140‑1 (MQ=255)                                                                                                                             cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                 >  1:434300/1‑140 (MQ=255)                                                                                                                             cGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCa                 >  1:519683/1‑140 (MQ=255)                                                                                                                               ccccTTTGGGTTGTGGGTTTGGGTGCGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCAACCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaat               <  1:1062021/140‑1 (MQ=255)                                                                                                                                   tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg           >  2:395888/1‑140 (MQ=255)                                                                                                                                   tttGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATg           >  1:366930/1‑140 (MQ=255)                                                                                                                                    ttGGGGTGAGGGTTTGGGTGTGGGTAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGc          <  1:1352030/140‑1 (MQ=255)                                                                                                                                          agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa    >  1:1595751/1‑140 (MQ=255)                                                                                                                                          agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa    >  2:878597/1‑140 (MQ=255)                                                                                                                                          agagGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCa    >  2:879064/1‑140 (MQ=255)                                                                                                                                            agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  2:168042/1‑140 (MQ=255)                                                                                                                                            agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:1246664/1‑140 (MQ=255)                                                                                                                                            agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:1238695/1‑140 (MQ=255)                                                                                                                                            agGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCAtt  >  1:763336/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                            GCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATT  >  NC_000913/624459‑624737

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