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

CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421298‑2421592                                                                                                                                                  |                                                                                                                                                      ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTAcccc                                                                                                                                                   <  1:407372/150‑1 (MQ=255)  cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                                 >  1:702086/1‑151 (MQ=255)       tAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGAt                                                                                                                                             >  2:387796/1‑150 (MQ=255)                    aaCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACt                                                                                                                                 >  1:820344/1‑149 (MQ=255)                           atgggCACACCCTCCCTTTTTTATTCTAAAAAAAACCCTTCGTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTCAGGATATTGTTATAAAGGTGCCGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGaa                                                                                                                           <  2:271050/145‑1 (MQ=255)                           aaggccccacccTGCCTTTATTTTTCTTAAAAAAACCCTTCGTTGGCTTCCTTTTATCGAATTTAGGGAAAAAAAAATAAACACCCGTCAGGATTTTTTTAAAATTGTGACGCCCCACCCCCCCGCCGATAAAGCGCGCACTa                                                                                                                                <  2:48069/133‑1 (MQ=255)                            aTGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                        >  1:375405/1‑150 (MQ=255)                                  caACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  2:817834/1‑151 (MQ=255)                                   accctcccttttttattCGTAAAAAAAACCTTCTTTTGCTTACTTTTATCGGAATTAAGAAAAAAAAAATTAACACCCGTAAGGATATTTTTATAAATGTGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 <  1:258479/138‑1 (MQ=255)                                   acccTGCCGTTATTATTCGAAAAAAGAACCTTCGTTTGCTTATGTTTATCGGAATTAAGAAAGAAAAAATTAACCCCCGGCAGGATATGTTTATAAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  <  2:822779/147‑1 (MQ=255)                                   acccTGCCGTTATTATTCGAAAAAAAACCCTTCGTTTGCGTACTTTTATCGAAATAAAGGAAGAAGAAATTAACAACCGTCAGGATATTTTTATCAATGGGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                                   <  1:532809/146‑1 (MQ=255)                                   acccTCCCGTTTTTATTCGTAAAAAGACCCTTCGTTTGCGTACTTTTATCGTAATTAAGGAAGAAGAAATTAACACCCGTCAGAATATTTTTATCAAGGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                                   <  1:361117/146‑1 (MQ=255)                                   aaCCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                                 >  1:608735/1‑150 (MQ=255)                                     cttccctttattatTCGAAAAAAGAACTTTCGTTTGCTTAGTTTTATGGAAATAAAGGAAAAAAAAATAAAACCCCGACAGGATATTTTTATAAAGGGGACATCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCaa                                                                                                                   <  2:841966/138‑1 (MQ=255)                                      cTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             >  1:419912/1‑151 (MQ=255)                                          cgttttttttCTTAAAAAGACCCTTCGTTCACGTACTTTTATGGAAATTAAAAAAGAAAAAATTAACACCCGTCAGGTTTTTTTTCACAAGGCGACGATCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             <  1:139780/139‑1 (MQ=255)                                          cGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                         >  2:874926/1‑151 (MQ=255)                                             tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      >  1:464140/1‑151 (MQ=255)                                             tattatTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCACGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAt                                                                                                      >  2:428432/1‑151 (MQ=255)                                                tattCGTAAAAAGACCCTTCGTTTGCGTACTTTTATCGGAATAAAGGAAAAAGAAATAAACACCCGTCAGAATATTGTTATAAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                                    <  1:414441/150‑1 (MQ=255)                                                 attCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                                  >  1:45223/1‑151 (MQ=255)                                                      tataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAcc                                                                                             >  2:468943/1‑151 (MQ=255)                                                       ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            >  2:713629/1‑151 (MQ=255)                                                       aaaaaGACCTTTCTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAACCACCCGTAAGAATTTTTTAATAAATGTGACTACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              <  2:592158/147‑1 (MQ=255)                                                       aaaaaGACCCTTGGTTTGCGTACTTTTATGGAAATTAAGAAAAAAAAAATAAACACCCGTCAGGATATTGTTACCAATGGGACGACCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              <  1:193490/147‑1 (MQ=255)                                                          aaGAACTTTTGTTGGCGATCGTTGACCGAAATAAAGGAAAAAAAAATCAACAACCCAGGGATTATGTTTAGTAATTCGACATTCCACCGCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTaa                                                                                         <  1:415120/151‑1 (MQ=255)                                                           aGACCCTTGTTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGCAAGAATTTTTTTATAAATGTGACACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        <  2:170199/151‑1 (MQ=255)                                                               cctttttttGCGTCCTTTTATCGAAATAAAGAAAAAAAAAATAAACACCCGTAAGAATTTTTTAATAAATGTGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                    <  2:900649/145‑1 (MQ=255)                                                               ccTTCTTTTGCGTACTTTTATGGAAATAAAGGAAAAAAAAATAAACACCCGTAAGGATATTGTTATAAATGGGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                    <  1:558373/151‑1 (MQ=255)                                                                     ttcgcgtgtgTTTGTCGGAATAAAGGAAAAAAAAACCAACAACCGGGAGAATATTATTATAAATGTAACCACCCCTTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 <  1:497059/139‑1 (MQ=255)                                                                          gtcgttttATGAGAATAACGGAAGAAGAAATAAACACCCTCAAGAATATGTTTAAAACTGGGAGAACCCACTGCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  1:22254/145‑1 (MQ=255)                                                                           tACTTTTATCGAAATTAAGAAAAAAAAAATAAACACCCGTCAGAATATTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  1:860700/149‑1 (MQ=255)                                                                            acttttATCGAAATAAAGAAAAAAAAAATTAACACCCGTCAGAATATTGTTATAAAGGAGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAg                                                                        <  1:430528/147‑1 (MQ=255)                                                                            aCGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                       >  2:178864/1‑151 (MQ=255)                                                                                ttATGGAAATAAAGGAAAAAAAAATAAACACCCGTCGGAATTTTTTTACAAATGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                    <  1:571277/150‑1 (MQ=255)                                                                                ttATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTGTCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   >  2:87588/1‑151 (MQ=255)                                                                                taatggaaaTTAAGAAAGAAGAAATTAACACCCGTCAGAATATTTTTATAAGGGGGACGGTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                    <  1:323249/143‑1 (MQ=255)                                                                                 tATCGGAATAAAGAAAAAAAAAATAAACAACCGTCAGAATTTTGTTAAAAATGTGACTACCCCCAACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                   <  2:777417/150‑1 (MQ=255)                                                                                   tggaaaTAAAGAAAAAAAAAATAAACACCCTTCAGAATTTTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTa                                                                  <  1:387603/145‑1 (MQ=255)                                                                                       aaTTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             >  2:114699/1‑150 (MQ=255)                                                                                       aaTAAAGAAAAAAAAAATAAACACCCGTAAGTATTTTTTAATAAATGGGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                            <  2:375367/151‑1 (MQ=255)                                                                                               aaaaaaaaaTAAACACCCGTCAGAATTTTGTTATCAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                       <  1:809192/142‑1 (MQ=255)                                                                                               aaaaaaaaaTAAACACCCGTAAGAATTTTTTTATAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    <  2:419870/145‑1 (MQ=255)                                                                                                 gaagaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                                  >  1:522574/1‑151 (MQ=255)                                                                                                 aaaaaaaTTAACCCCGGTCAGGATTTTTTTATAAAGGTGACTCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                                  <  2:16445/147‑1 (MQ=255)                                                                                                   agaaATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                >  2:841252/1‑151 (MQ=255)                                                                                                    aaaaTAAACCCCCGTAAGAATTTTGTTATAAATGTGACGTCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:940661/150‑1 (MQ=255)                                                                                                     aaaTAAACACCCGTCAGAATATTGTTATCAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               <  2:342051/150‑1 (MQ=255)                                                                                                      aaTTAACACCCGTAAGTATTTTGTTATAAAGGAGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                             <  2:441656/151‑1 (MQ=255)                                                                                                        ttAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                           >  2:642761/1‑151 (MQ=255)                                                                                                          aacaccCGTCAGAATATTTTTATAAATGTGACGACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGc                                            <  2:916681/148‑1 (MQ=255)                                                                                                            caacCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                        >  1:396832/1‑150 (MQ=255)                                                                                                               ccGTCAGGATTTTGTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     <  1:437913/150‑1 (MQ=255)                                                                                                               ccGTCAGGATTTTGTTATAAATGGGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                    <  1:973107/151‑1 (MQ=255)                                                                                                               ccGTCAGAATTTTTTTATAAATGTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGa                                        <  2:883751/147‑1 (MQ=255)                                                                                                                cGTAAGGATTTGTTTATAAAGGTGACACCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      <  2:554942/148‑1 (MQ=255)                                                                                                                  acAGGATTGTGTTACAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCttt                                              <  1:125888/137‑1 (MQ=255)                                                                                                                        tatTGTTATAAAGGGGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                             <  2:576474/149‑1 (MQ=255)                                                                                                                          ttGTTATAAATGTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACt                              <  1:431962/146‑1 (MQ=255)                                                                                                                                caaaTGTGCCCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAt                        <  1:322523/144‑1 (MQ=255)                                                                                                                                  aaTGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                    <  1:431514/148‑1 (MQ=255)                                                                                                                                       gACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc            >  2:634632/1‑151 (MQ=255)                                                                                                                                              cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:654880/1‑151 (MQ=255)                                                                                                                                               aCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  1:86479/1‑150 (MQ=255)                                                                                                                                                 tACCCCGCCGCTAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  2:757959/1‑151 (MQ=255)                                                                                                                                                 tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  1:846061/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421298‑2421592

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