Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,437	0	T	C	52.4%	50.3 / 84.9	63	I38T (ATC→ACC)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base T (22/8);  new base C (0/33);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 1.11e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-09

TATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATC  >  NC_000913/2421287‑2421582                                                                                                                                                       |                                                                                                                                                  tattatTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGAt                                                                                                                                                   >  1:640554/1‑151 (MQ=255)       ttCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACt                                                                                                                                              >  2:806294/1‑150 (MQ=255)            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)                                                                                                                                                       |                                                                                                                                                  TATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATC  >  NC_000913/2421287‑2421582

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