Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,421,426	0	C	A	49.2%	59.0 / 32.2	62	I34I (ATC→ATA)	folX	D‑erythro‑7,8‑dihydroneopterin triphosphate 2'‑epimerase and dihydroneopterin aldolase
Reads supporting (aligned to +/- strand):  ref base C (22/9);  new base A (0/30);  total (22/40)
Fisher's exact test for biased strand distribution p-value = 1.17e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTC  >  NC_000913/2421276‑2421565                                                                                                                                                       |                                                                                                                                            aGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATc                                                                                                                                             >  2:422968/1‑151 (MQ=255)  ggTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCa                                                                                                                                            <  1:498867/151‑1 (MQ=255)            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)                                                                                                                                                       |                                                                                                                                            AGGTTCTCGTGTATTATTTCATCCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATTAACAACCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTC  >  NC_000913/2421276‑2421565

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: ‑