Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,443	0	A	C	40.3%	131.1 / 19.3	67	Y40S (TAC→TCC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base A (39/1);  new base C (0/27);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 6.82e-18
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592                                                                                                                                                       |                                                                                                                                                      tCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAc                                                                                                                                                      >  1:872979/1‑151 (MQ=255) tCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAc                                                                                                                                                      >  2:577991/1‑151 (MQ=255)  cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTa                                                                                                                                                       <  1:97173/149‑1 (MQ=255)       tAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgcc                                                                                                                                                >  1:200394/1‑151 (MQ=255)        aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgcc                                                                                                                                                <  2:364308/150‑1 (MQ=255)                  acGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcg                                                                                                                                     >  1:783425/1‑151 (MQ=255)                   cGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                    >  1:308392/1‑151 (MQ=255)                     gagaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCAc                                                                                                                                  >  1:485441/1‑151 (MQ=255)                       gaACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  2:414043/1‑150 (MQ=255)                        aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  1:733939/1‑149 (MQ=255)                         ccctaaaaggccacaCCCTGCCGTTATTTTTCGTATAAAGACCCTTCTTTTGCTTAC‑TTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGTTATTTTTATCAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGc                                                                                                                              <  2:795230/140‑1 (MQ=255)                           cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             >  2:731012/1‑150 (MQ=255)                              cAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  1:367296/1‑148 (MQ=255)                               aaTGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  1:328031/1‑147 (MQ=255)                               aaTGCCACACCCTGCCTTTTTTATTCGAAAAAAGACCCTTTGTTTGCTTAC‑GTTAATCGAAATTAAGAAAAAAAAAATAAACACCCGTAAGAATTTTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                                                                                         <  2:491638/150‑1 (MQ=255)                                 tGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                      >  2:191915/1‑151 (MQ=255)                                     acaccccgcctttattatTCTAAAAAAGAACCTTCTTTTGCGTTC‑GTTTATCGGAATAAAGGAAAAAAAAATAAACACCCTTCAGAATATTTTTATAAATGTGACGCTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                                   <  2:978425/138‑1 (MQ=255)                                       acccTGCCGTTATTATTCGTAAAAAGACCCTTCGTTTGCTTAC‑GTTTATCGAAATAAAGGAAGAAAAAATAAACACCCGTCAGGTTTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  <  2:726174/147‑1 (MQ=255)                                        aCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               >  1:553843/1‑151 (MQ=255)                                             cctttttttttCGTAAAAAGACCCTTCGTTTGCTTAT‑TTTTATGGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGATTTTGTTATAAATGTGCCGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                          <  1:927493/142‑1 (MQ=255)                                             cctttattatTAATAAAAAAACCTTTCTTTTGCGATC‑GTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGTTTTTTTTAACAAGGTAACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             <  1:249705/144‑1 (MQ=255)                                             ccacTATTTTTGGTAAAAAAATCTTTCTTTTGCGTTC‑GTTTTATGAAATAAGGGAAAAAAAAATAAACACCCGTCAGGATATTGTTAAAAGGGGGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             <  2:998783/145‑1 (MQ=255)                                              cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                         >  2:748665/1‑151 (MQ=255)                                                 tattatTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      >  1:658744/1‑151 (MQ=255)                                                 tattatTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      >  2:971560/1‑151 (MQ=255)                                                       tcataaaaaGAACCTT‑GTTTGGGTTCTTTTTATGGAAATAAAGAAAGAAAAAATAAAAACCCGTAGGAATTTTGTTATAAATGTGCCTCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                                 <  2:743021/144‑1 (MQ=255)                                                        cgaataAAAACCCTTCGTTTGCGTAC‑TTTTATGGGAATTAAGAAAAAAAAAATAAACACCCGTCAGGATTTTTTTACCATTGTGACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:958064/148‑1 (MQ=255)                                                             aaaGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:605955/1‑150 (MQ=255)                                                             aaaGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:610899/1‑150 (MQ=255)                                                                   ctttCTTTTGCTTCC‑GTTTACCGGATTAAAGGAAAAAAAAATAAACACCCGTCAGTATTTTGTTAAAATGGTGACGCCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      <  1:954854/147‑1 (MQ=255)                                                                        gTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               >  2:938083/1‑151 (MQ=255)                                                                        gTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAAGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                               >  1:1005041/1‑151 (MQ=255)                                                                           ttggtttgttttATGGAAATTAAAGAAAAAAAAATAAACACCCCGAAGGTTTTTGTTATAAATGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                                <  1:219175/139‑1 (MQ=255)                                                                             cGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  2:257277/1‑151 (MQ=255)                                                                                 gtttttACCGAAATAAGGAAAGAAAAAATAAACACCCGAAGGGATGTTTTTATTAATGGGACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                        <  1:377658/147‑1 (MQ=255)                                                                                  tgtttacagaaaTAAAGGAAGAAAAAATCAAAACCTGTCGGGATTTTGTTAACATTGTTGTCACCCTTTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                        <  2:814245/139‑1 (MQ=255)                                                                                   gTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCAt                                                                     >  2:867833/1‑151 (MQ=255)                                                                                    tttATCGAAATTAAGGAAAAAAAAATTAACACCCGTCAGGATATTTTTATAAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                                      <  2:239773/149‑1 (MQ=255)                                                                                        tCGAAATAAAGAAAAAAAAAATTAACAACCGTCAGAATATTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 <  2:455297/150‑1 (MQ=255)                                                                                          gaaaTAAAGAAAAAAAAAATTAACACCCGTAAGGATATTTTTATAAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                                <  2:828981/147‑1 (MQ=255)                                                                                              taaaGGAAAAAAAAATAAACACCCGTAAGGATTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             <  2:535038/146‑1 (MQ=255)                                                                                               tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         >  2:286256/1‑151 (MQ=255)                                                                                               tAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         >  1:524831/1‑151 (MQ=255)                                                                                                 aGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                       >  1:201157/1‑151 (MQ=255)                                                                                                  ggAAAAAAAAATCACCACCCGGCAGTATATTGTTATAAATGTGACATTCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATaa                                                          <  2:49945/147‑1 (MQ=255)                                                                                                   aaaaaaaaaaTAAACACCCGTAAGGTTTTTTTTAAAAATGTAACAACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                       <  1:67517/142‑1 (MQ=255)                                                                                                    aagaagaaCTAAACCACCACCAGGTTATTATTATCAATGTGGCGATCATGACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAGCCGtt                                                     <  1:735976/150‑1 (MQ=255)                                                                                                     agaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                   >  2:311060/1‑151 (MQ=255)                                                                                                       aaaaaaTTAACACCCGTAAGAATTTTTTTATCAATGGGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                                  <  2:129734/147‑1 (MQ=255)                                                                                                        agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                >  1:663849/1‑151 (MQ=255)                                                                                                        agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                >  1:679989/1‑151 (MQ=255)                                                                                                          aaaTTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGTACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                              >  2:903349/1‑151 (MQ=255)                                                                                                           aaatAACACCCGTCAGGATATTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                             <  2:733956/148‑1 (MQ=255)                                                                                                               aacaccCGTAAGGATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                           <  1:870142/149‑1 (MQ=255)                                                                                                               aacaccCGTAAGGATTTTGTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         <  2:133073/151‑1 (MQ=255)                                                                                                               aacaacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTgg                                         >  1:851462/1‑151 (MQ=255)                                                                                                                        cAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTa                                 >  1:601246/1‑150 (MQ=255)                                                                                                                         aGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAAc                               >  2:324340/1‑151 (MQ=255)                                                                                                                           gATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                              >  2:11458/1‑150 (MQ=255)                                                                                                                             tttttttATAAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                                   <  1:888479/138‑1 (MQ=255)                                                                                                                             tatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  2:851235/1‑149 (MQ=255)                                                                                                                               tttttaaaaaTGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTaa                                <  1:840116/136‑1 (MQ=255)                                                                                                                               ttGTTATCAATGACACACCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAgg                          <  2:72068/150‑1 (MQ=255)                                                                                                                                   tATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     >  2:684210/1‑151 (MQ=255)                                                                                                                                     tCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCt                   >  2:494097/1‑151 (MQ=255)                                                                                                                                               gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  1:690350/1‑150 (MQ=255)                                                                                                                                                   cACTACCCCCCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:792689/1‑151 (MQ=255)                                                                                                                                                     cTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGa    >  2:398280/1‑150 (MQ=255)                                                                                                                                                       aCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:305051/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                      TCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTAC‑GTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421294‑2421592

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑