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

CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421295‑2421592                                                                                                                                                      |                                                                                                                                                      cATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTAcc                                                                                                                                                     <  2:397057/151‑1 (MQ=255)    ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccg                                                                                                                                                  >  2:700147/1‑151 (MQ=255)      tAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgcc                                                                                                                                                >  1:352289/1‑151 (MQ=255)         gTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGAt                                                                                                                                             >  2:713437/1‑151 (MQ=255)           aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATa                                                                                                                                            >  2:63781/1‑150 (MQ=255)              acaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgc                                                                                                                                        >  2:233584/1‑151 (MQ=255)                aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                                      >  1:711017/1‑151 (MQ=255)                         ccTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             >  2:145846/1‑151 (MQ=255)                           tGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  2:383362/1‑150 (MQ=255)                             cAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGa                                                                                                                            >  2:123286/1‑148 (MQ=255)                             aaaggccaaacccTGCCGTTTTTATTCAAAAAAAGACCCTGAGTTTGCGTCCTTTTATCGAAATAAAGAAAGAAAAAAAAAACAC‑CCGTAAGAATTTTGTTAGCAAGGGGACGGCCCACCACCCCGCCGATAAAGCGCGCACTAGCGAAg                                                                                                                          <  1:63459/139‑1 (MQ=255)                                      aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  1:623838/1‑150 (MQ=255)                                        ccggCCGTTATTATTCTTAAAAAAACCTTTCTTTTGCGTACTTTTATGGAAATAAAGAAAAAAAAAATAAACAC‑CCGTAAGGATATTTTTAGAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                                                                              <  1:38114/148‑1 (MQ=255)                                        ccTGCGGTTTTTTTTCGAAAAAAAAACCTTCTTTTCCTTCCTTTTATGGAAATTAAGGAAAAAGAAATTACCAC‑CCGTAAGGTTTTTTTTACCAAGGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             <  1:298521/136‑1 (MQ=255)                                        ccTGCCGTTTTTATTCTTAAAAAGACCTTTCTTTTGCTTCCGTTTATGGAAATAAAGAAAAAAAAAATAAACAC‑CCGTCGGAATTTTTTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  2:343552/150‑1 (MQ=255)                                         cTGCCGTTTTTATTCGTAAAAAAACCCTTCTTTTGCTTCCTTTTATCGAAATAAAGGAAGAAAAAATTAACAC‑CCGTCAGGATTTTGTTATAAATGTGACAATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                               <  2:726436/149‑1 (MQ=255)                                            ccttttttattCATAAAAAAACCCTTCTTTTGCGTCTTTTTATAGAAATAAAGAAAAAAAAAATAAACAC‑CCGAAAGGATTTTTTTAAAAATGAGACACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                          <  1:210332/145‑1 (MQ=255)                                            ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                          >  1:677412/1‑151 (MQ=255)                                                tattatTCGTAAAAAAACCCTTCTTTTGCGTACTTTTATCGGAATTAAGGAAGAAAAAATAAACAC‑CCGTCAGGATATTTTTAACAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                        <  2:135558/149‑1 (MQ=255)                                                  ttattCTAAAAAAAACCCTTCTTTTGCTTATTTTTATCGAAATAAAGGAAAAAAAAATAACCAC‑CCGTAAGGATTTTTTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      <  1:62054/149‑1 (MQ=255)                                                   tattCGTAAAAAGACCCTTCGTTTGCGTACGTTTATCGAAATTAAGGAAAAAGAAATAAACAA‑CCGTCAGGATATTGTTATAAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                                    <  2:623899/150‑1 (MQ=255)                                                        gTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              >  2:788473/1‑151 (MQ=255)                                                          aaaaaaaccttttgTTGGCTTACTTTTATCGAAATAAAGAAAAAAAAAATAACCAC‑CCGTAAGGATATTTTTATAAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                                <  1:310677/134‑1 (MQ=255)                                                          aaaaaaacccTTCTTTTGCGTATTTTTATCGAAATAAAGGAAAAAAAAATTAACAC‑CCGCCAGAATGTTGTTAAAATTGAGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCg                                                                                                 <  1:854649/138‑1 (MQ=255)                                                          aaaaaaacccTTCGTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATTAACAC‑CCGTCAGGATTTTTTTATCAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:911664/140‑1 (MQ=255)                                                          aaaaaaaaCCTTTTTTGGCGTTCTTTTATCGAAATTAAGAAAAAAAAAATAAACAC‑CCGTCAGGATTTTTTAATCAATGTGCCAATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              <  1:804446/147‑1 (MQ=255)                                                            aaaGACCCTTCGTTTGCGTACGTTTATGGAAATAAAGAAAAAAAAAATTAACAC‑CCGTCAGGATTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                          <  1:847559/151‑1 (MQ=255)                                                            aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  1:767490/1‑150 (MQ=255)                                                             aaGACCCTTCGTTTGCGTACGTTTATGGGAATAAAGAAAAAAAAAATAAACAC‑CCGGCAGGTTTTTGTTAACAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                         <  2:362855/151‑1 (MQ=255)                                                              aGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        >  1:482329/1‑151 (MQ=255)                                                                aaCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                       >  1:668502/1‑150 (MQ=255)                                                                  ccTTCGTTTGCGTTCTTTTATCGAAATAAAGGAAAAAAAAATTAACAA‑CCGTCAGGATATTTTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaa                                                                                     <  1:66506/150‑1 (MQ=255)                                                                   cTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                    >  2:919186/1‑150 (MQ=255)                                                                    ttgtttggCGTACATTTATCGAAATTAAGGAAGAATAAAAAAACAC‑CGATAGGGATATTATTAACAAGGTGACCCCCCCCGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      <  2:153163/140‑1 (MQ=255)                                                                    ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAGGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                   >  1:623252/1‑150 (MQ=255)                                                                    ttCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaaa                                                                                   >  1:461768/1‑150 (MQ=255)                                                                     tCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 >  2:61578/1‑151 (MQ=255)                                                                         ttGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              >  2:361306/1‑150 (MQ=255)                                                                         ttGCGTACGTTTAACGAAATAAAGAAAGAAAAAATAAACAC‑CCGTAAGAATTTTGTTATCAATGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACatt                                                                              <  2:796867/150‑1 (MQ=255)                                                                             gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  1:698434/1‑150 (MQ=255)                                                                             gTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  2:322514/1‑150 (MQ=255)                                                                              tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  2:528734/1‑149 (MQ=255)                                                                              tACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  2:734565/1‑149 (MQ=255)                                                                               aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                       >  2:586569/1‑151 (MQ=255)                                                                                     aTCGAAATAAAGAAAAAAAAAATAAACAC‑CCGTCAGAATTTTGTTATAAATGGGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 <  2:629153/151‑1 (MQ=255)                                                                                        ggAATAAAGGAAAAAAAAATAAACAC‑CCGTCAGAATATTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                                  <  2:308389/147‑1 (MQ=255)                                                                                        gaaaTTAAGGAAGAAGAAATTAACAA‑CCTTAAGGATTTTTTTATAAATGTGACGTTCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTaga                                                                <  2:368785/147‑1 (MQ=255)                                                                                          aaTAAAGGAAAAAAAAATAAACAC‑CCGTAAGAATTTTTTTACAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 <  2:591901/146‑1 (MQ=255)                                                                                              aaGAAAAAAGAAGTCAACAC‑CCGAAAGGATGTTATTAGCAATGTGACTTCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                        <  1:716400/151‑1 (MQ=255)                                                                                               aagaaaaaaaaAATAACAACCCATAAGGATTTTTTTATAAATGTCACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                        <  1:306589/149‑1 (MQ=255)                                                                                                  aagaaaaaATAACCAC‑CCGTCAGAATATTTTAACAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGtt                                                     <  2:310798/150‑1 (MQ=255)                                                                                                  aaaaaaaaaTAAACAC‑CCGTCAGGATTTTGTTATCAATGGGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                       <  1:456710/142‑1 (MQ=255)                                                                                                   agaagaaAGTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                   >  1:856534/1‑151 (MQ=255)                                                                                                   agaaaaaATTACCAC‑CTGTAAGTTTTTTTTTATCAATGTGCCATCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtc                                                   <  1:63778/151‑1 (MQ=255)                                                                                                     aaaaaaTAAACAC‑CCTTCAGGATTTTTTTATAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtctc                                                 <  1:630133/148‑1 (MQ=255)                                                                                                      agaaATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                >  1:409592/1‑151 (MQ=255)                                                                                                       aaaaTTAACAA‑CCGTCAGAATTTTTTTATCAATGTGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  1:471113/150‑1 (MQ=255)                                                                                                       aaaaTAAACAC‑CCGTCAGGATATTTTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:688428/150‑1 (MQ=255)                                                                                                        aaaTAAACAC‑CCGTAAGGATTTTTTTATCAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                              <  1:534494/151‑1 (MQ=255)                                                                                                           ttAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                           >  2:776839/1‑151 (MQ=255)                                                                                                           ttAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                           >  1:674436/1‑151 (MQ=255)                                                                                                                aa‑cCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaa                                       >  1:390938/1‑150 (MQ=255)                                                                                                                     gtaaGGTTATTTTTATCAATGTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAtt                                  <  2:316094/148‑1 (MQ=255)                                                                                                                         gattatTGTTAACACCGGGACGACCCCCTCTCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATGATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                        <  2:500574/138‑1 (MQ=255)                                                                                                                           atatTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             >  1:523875/1‑150 (MQ=255)                                                                                                                                     aaaTGTGCCGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                    <  1:448618/148‑1 (MQ=255)                                                                                                                                            aCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat             <  2:708703/149‑1 (MQ=255)                                                                                                                                              gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  2:39468/1‑150 (MQ=255)                                                                                                                                                tCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACg       >  1:352495/1‑151 (MQ=255)                                                                                                                                                      aCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:333349/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      CATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421295‑2421592

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