Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,400,336	C→A	R17R (CGG→AGG)	hfq →	RNA‑binding protein Hfq

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,400,336	0	C	A	98.9%	330.5 / ‑2.7	93	R17R (CGG→AGG)	hfq	RNA‑binding protein Hfq
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (48/43);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-01

GAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACT  >  NC_000913/4400201‑4400470                                                                                                                                        |                                                                                                                                       gAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGa                                                                                                                                    <  1:665383/140‑1 (MQ=255)        tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGtt                                                                                                                             >  2:657937/1‑140 (MQ=255)        tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGtt                                                                                                                             >  2:1085382/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATGAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        >  1:16487/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        <  1:413893/140‑1 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        >  2:839183/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        >  2:690030/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        >  2:653964/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        >  1:637822/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        >  2:1205573/1‑140 (MQ=255)             ttGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGt                                                                                                                        <  1:976265/140‑1 (MQ=255)                agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTc                                                                                                                     >  2:532042/1‑140 (MQ=255)                agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTc                                                                                                                     >  1:289808/1‑140 (MQ=255)                  aCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCta                                                                                                                   >  2:1251682/1‑140 (MQ=255)                  aCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCta                                                                                                                   >  2:937849/1‑140 (MQ=255)                    gTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATGAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCtatt                                                                                                                 <  2:16487/140‑1 (MQ=255)                     tATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCtattt                                                                                                                >  2:604146/1‑140 (MQ=255)                     tATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCtattt                                                                                                                >  1:428647/1‑140 (MQ=255)                       tCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCtatttat                                                                                                              >  1:690517/1‑140 (MQ=255)                       tCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCtatttat                                                                                                              >  1:1046834/1‑140 (MQ=255)                          tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTg                                                                                                           <  2:1008893/140‑1 (MQ=255)                          tGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTg                                                                                                           <  2:560569/140‑1 (MQ=255)                               aaTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGTTCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                                                                                      <  2:595823/140‑1 (MQ=255)                               aaTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                                                                                      <  1:839183/140‑1 (MQ=255)                               aaTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                                                                                      <  1:552711/140‑1 (MQ=255)                               aaTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                                                                                      <  2:893449/140‑1 (MQ=255)                               aaTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACCTGTTCCAGTTTCTATTTATTTGGTGaa                                                                                                      <  2:835195/140‑1 (MQ=255)                                aTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                                                                                     >  2:985250/1‑140 (MQ=255)                                                  ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGT‑AGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                  <  2:524301/140‑1 (MQ=255)                                                   gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                  <  1:985250/140‑1 (MQ=255)                                                    ttCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCa                                                                                 >  2:1075184/1‑140 (MQ=255)                                                          gTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGa                                                                           <  2:954791/140‑1 (MQ=255)                                                          gTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGa                                                                           <  1:937849/140‑1 (MQ=255)                                                             cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc                                                                        >  1:195422/1‑140 (MQ=255)                                                             cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc                                                                        >  2:1148810/1‑140 (MQ=255)                                                                   aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                  >  2:945705/1‑140 (MQ=255)                                                                   aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                  >  1:93791/1‑140 (MQ=255)                                                                   aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                  >  1:858227/1‑140 (MQ=255)                                                                   aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                  >  1:460226/1‑140 (MQ=255)                                                                   aGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                  >  1:318686/1‑140 (MQ=255)                                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 <  2:389555/140‑1 (MQ=255)                                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 <  1:757566/140‑1 (MQ=255)                                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 <  1:737812/140‑1 (MQ=255)                                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 <  1:653970/140‑1 (MQ=255)                                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 <  1:276371/140‑1 (MQ=255)                                                                    gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGACCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                 <  2:637822/140‑1 (MQ=255)                                                                     cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGGTCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                <  2:68383/140‑1 (MQ=255)                                                                     cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                >  2:160352/1‑140 (MQ=255)                                                                     cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                <  2:858227/140‑1 (MQ=255)                                                                     cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                >  1:210880/1‑140 (MQ=255)                                                                      atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCa                                                               <  2:428647/140‑1 (MQ=255)                                                                      atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCa                                                               <  1:532042/140‑1 (MQ=255)                                                                         taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGtt                                                            >  2:173443/1‑140 (MQ=255)                                                                              aaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGa                                                       >  2:532467/1‑140 (MQ=255)                                                                              aaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGa                                                       >  1:159952/1‑140 (MQ=255)                                                                               aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                      >  2:594057/1‑140 (MQ=255)                                                                               aaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGAt                                                      >  1:672183/1‑140 (MQ=255)                                                                                 agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                    <  2:256144/140‑1 (MQ=255)                                                                                 agagagAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATcc                                                    <  2:318686/140‑1 (MQ=255)                                                                                     agaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                <  1:1251682/140‑1 (MQ=255)                                                                                       aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGa                                              >  1:266057/1‑140 (MQ=255)                                                                                       aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGa                                              >  1:1021503/1‑140 (MQ=255)                                                                                                 gggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGt                                    >  1:805333/1‑140 (MQ=255)                                                                                                 gggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGt                                    <  2:289808/140‑1 (MQ=255)                                                                                                  ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTc                                   <  2:672183/140‑1 (MQ=255)                                                                                                   gCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCa                                  >  1:735636/1‑140 (MQ=255)                                                                                                     aaTCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGc                                >  2:52065/1‑140 (MQ=255)                                                                                                     aaTCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGc                                >  1:2811/1‑140 (MQ=255)                                                                                                      aTCTTTACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGcc                               >  1:850681/1‑140 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  1:960762/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  1:1205573/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  2:690517/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  2:321935/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  1:653964/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  2:1166977/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  1:657937/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  1:160352/140‑1 (MQ=255)                                                                                                         tttACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGa                            <  2:383211/140‑1 (MQ=255)                                                                                                           tACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCGGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                          >  2:497596/1‑140 (MQ=255)                                                                                                           tACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAg                                 <  2:602329/133‑1 (MQ=255)                                                                                                           tACAAGATCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAg                                 >  1:602329/1‑133 (MQ=255)                                                                                                                 aTCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                    >  1:1111610/1‑140 (MQ=255)                                                                                                                 aTCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                    >  1:653516/1‑140 (MQ=255)                                                                                                                 aTCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                    >  1:802525/1‑140 (MQ=255)                                                                                                                 aTCCGTTCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTAc                    >  2:338707/1‑140 (MQ=255)                                                                                                                       tCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAc              <  1:52065/140‑1 (MQ=255)                                                                                                                       tCCTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAc              <  1:338707/140‑1 (MQ=255)                                                                                                                         cTGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgc            <  2:1111610/140‑1 (MQ=255)                                                                                                                          tGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcg           <  1:1148810/140‑1 (MQ=255)                                                                                                                          tGAACGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcg           <  1:690030/140‑1 (MQ=255)                                                                                                                            aaCGCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGAt         >  1:1118160/1‑140 (MQ=255)                                                                                                                               gCACTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc      <  2:735636/140‑1 (MQ=255)                                                                                                                                 aCTGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTACGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTa    >  1:794735/1‑140 (MQ=255)                                                                                                                                   tGCGTAGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt  >  1:26153/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       GAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACT  >  NC_000913/4400201‑4400470

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