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

CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421298‑2421589                                                                                                                                                     |                                                                                                                                                   ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccg                                                                                                                                               >  1:217349/1‑151 (MQ=255)  cTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCccgc                                                                                                                                              >  2:708244/1‑151 (MQ=255)     aGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGa                                                                                                                                           <  2:946647/151‑1 (MQ=255)        aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATa                                                                                                                                         >  2:78127/1‑150 (MQ=255)        aaaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATa                                                                                                                                         >  1:237408/1‑150 (MQ=255)           acaacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAg                                                                                                                                      >  1:693232/1‑150 (MQ=255)               cGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCTGCCGATAAAgcgcgc                                                                                                                                 >  2:286050/1‑151 (MQ=255)                ggAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAgcgcgc                                                                                                                                 >  1:40876/1‑150 (MQ=255)                  agaCCCTGAAAGGCCACACCCTGCCATTATTTTTAGTAAAAAGACCCTTCGTTTGCTTTCTTTTAGCGA‑AATTAAGA‑AAAAAAAAATAAACAC‑CGGTCAGGATATTGTTATAATTGTGCCGATCCCCTCCCCCGCCGATAAAGCGCGCACt                                                                                                                              <  2:759262/151‑1 (MQ=255)                    acccTGAAATGACACACCCTGCCTTTTTTATTCAAAAAAAGACCCTTCGTTTGCGTATTTTTATCGA‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTAAGGATTTTTTTATAAAGGTAACGACCCCCTCCCCCGCCGATAAAGCGCGCACt                                                                                                                              <  2:166177/147‑1 (MQ=255)                     aCCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGc                                                                                                                           >  2:118124/1‑151 (MQ=255)                      ccTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                          >  1:336160/1‑151 (MQ=255)                           aatgcaacacCCTGCCTTTTTTTTTCATAAAAAGACCCTTGTTTGGCGTTTTTTTATCAA‑AATAAAGA‑AAAAAAAAATAACCAC‑CCGTCAGAATTTTGTTAACAATGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                     <  2:749165/145‑1 (MQ=255)                           aaggccacaCCCTGCCGTTTTTTTTCGTAAAAAAAACTTTCTTTTGCTTATTTTTATCGA‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTCAGGATATTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGa                                                                                                                      <  1:207821/145‑1 (MQ=255)                            agggCACATCCTGCCTCTTTTATTCGAAAAAAGACCCTTCTTTTGCTTACTTTTATCGA‑AATAAAAA‑AAAAAAAAATAAACAC‑CCTCCAGAATTTTTTAAAAAATGTGACGCCCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGat                                                                                                                     <  1:905138/148‑1 (MQ=255)                                     cccgcctttattatTCGTAAAAAGACCCTTCTTTTGCGTCCTTTTATGGA‑AATAAAAG‑AAAAAAAAATAAACAC‑CCGTGAGGATATTTTTAAAAATGTGCCGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGatat                                                                                                                   <  2:611246/135‑1 (MQ=255)                                     ccTGCCTTTATTATTCGAAAAAAAACCTTTCGTTTGCGTACTTTTACCGA‑AATAAAGG‑AAGAAAAAATAAACAC‑CCGTCAGGATATTTTTAAAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACg                                                                                                              <  1:423364/148‑1 (MQ=255)                                     ccTGCCGTTTTTATTCGAAAAAAAACCTTTCGTTGGCGTACTTTTATCGA‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTAAGGATTTTTTTAACAAGGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGAt                                                                                                            <  2:540335/150‑1 (MQ=255)                                         ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                       >  2:633875/1‑151 (MQ=255)                                          cgtttttattCTAAAAAAGACCCTTCGTTTGCTTGCTTTTAGGGA‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTAAGAATTTTTTTAAAATTGTGCCGATCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                       <  2:799927/145‑1 (MQ=255)                                           gCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      >  2:901885/1‑150 (MQ=255)                                             ttttattCGAAAAAAGACCCTTCTTTTGCGTTGTTTTATCGA‑AATAAAGG‑AAGAAAAAATTAACAC‑CCGGAAGGATTTTGTTATAAATGGGACTACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTg                                                                                                      <  1:886218/146‑1 (MQ=255)                                             ttttattCATAAAAAAACCCTTCGTTTGCGTACGTTTATCGG‑AATTAAGA‑AAAAAAAAATAAACAC‑CCGTCAGTATATTTTTATAAATGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                     <  1:464108/147‑1 (MQ=255)                                               ttattaGAAAAAAGACCCTTCGTTTGCGTACTTTTATCGA‑AATAAAGA‑AAAAAAAAATTAACAC‑CCGTCAGGATTTTTTTAAAAATGTGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgcg                                                                                                         <  1:206982/143‑1 (MQ=255)                                               ttattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAAtt                                                                                                  >  1:669970/1‑150 (MQ=255)                                                 attCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               >  1:449527/1‑151 (MQ=255)                                                  ttCGAATAAAGAACTTTCTTTTGCGT‑CTTTTTTGGGTAATTAAGG‑AAGAAAAAATAACCAG‑CCGTCAGGATTTTTTTATAAAGGTGACATCCCACTACACCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATc                                                                                               <  1:258188/150‑1 (MQ=255)                                                     gTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                           >  2:904306/1‑151 (MQ=255)                                                       ataAAGACCTTTTTTTGGCGTACTTTTATCAG‑AATGAAGG‑AAGAAAAAATAACCAC‑CCTTCAGTATTTTTTTATCAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                           <  1:246839/149‑1 (MQ=255)                                                       aaaaaGACCCTTCGTTTGCGTATTTTTATCGA‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTCAGGTTTTTTTTATCAATGGGCCACCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                         <  2:728206/149‑1 (MQ=255)                                                       aaaaaGACCCTTCGTTTGCGTACTTTTATCGA‑AATAAAGA‑AAGAAAAAATTAACAC‑CCGTAAGAATTTTTTTATAATTGGGACGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                             <  1:702423/145‑1 (MQ=255)                                                        taAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                        >  1:406856/1‑151 (MQ=255)                                                         aaaaacccTTCTTTTGCGTACGTTTATCGA‑AATAAAGA‑AAAAAAAAATAAACAA‑CCGTCAGGATATTGTTACCATTGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTa                                                                                       <  1:228729/145‑1 (MQ=255)                                                         aaaaacccTTCGTTTGCGTGCTTTTATCGG‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTCAGAATTTTTTTATAAATGTGCCGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                        <  2:506122/144‑1 (MQ=255)                                                           aGACCCGTCTTTGGCTTATGTTTATGAA‑AATAAAGA‑AAAAAAAAAAAAACAC‑CCGTGAGGATATTGTTATAAATGTGACCCTCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                     <  2:687313/151‑1 (MQ=255)                                                           aGAACCTTCGTTTGCTTATTTTTATGAA‑AATAAAGA‑AAAAAAAAATTACAAC‑CCGTAAGAATTTTTTTATAAATGTGCCGACCCCCACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                     <  2:137261/151‑1 (MQ=255)                                                              aCCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                   >  2:420061/1‑150 (MQ=255)                                                                cTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                 >  1:389722/1‑150 (MQ=255)                                                                  tcttttGCTTATTTTTATAGA‑AATAAAGA‑AAAAAAAAATAAACAC‑CCGTCAGGTTTTTTTTATCATTGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGaaa                                                                                 <  2:966882/145‑1 (MQ=255)                                                                    gTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                            >  1:894900/1‑151 (MQ=255)                                                                    gTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACat                                                                            >  2:634289/1‑151 (MQ=255)                                                                     tttGCGTCTTTTTATGGG‑ATATAA‑G‑AAAAAAAAATTAAAAACCCGTCAGGATTTTGTTATCAAGGTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACa                                                                             <  2:633654/149‑1 (MQ=255)                                                                     tttGCGTCTGTTTATAGG‑AATTAAGG‑AAAAAAAAATAAACAC‑CCGTCAGGATATTGTTATAAATGTGACAACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                              <  1:201099/148‑1 (MQ=255)                                                                     ttggCGTCCTTTTAT‑GG‑AAATAAGGAAAGAAAAAATAACCAC‑CCGTTAGGATATTGTTAACAATGTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACg                                                                                    <  2:693995/139‑1 (MQ=255)                                                                       tGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                         >  1:777608/1‑151 (MQ=255)                                                                         cGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                       >  1:199673/1‑151 (MQ=255)                                                                         cGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                       >  1:359923/1‑151 (MQ=255)                                                                          gTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                       >  2:924170/1‑150 (MQ=255)                                                                             cGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                    >  1:36731/1‑150 (MQ=255)                                                                               tttATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                 >  2:865449/1‑151 (MQ=255)                                                                                  aTCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                >  2:894746/1‑149 (MQ=255)                                                                                   tCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGt                                                                >  1:269211/1‑148 (MQ=255)                                                                                     gg‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                              >  1:188043/1‑148 (MQ=255)                                                                                      g‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                           >  1:907381/1‑150 (MQ=255)                                                                                      g‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                           >  1:970492/1‑150 (MQ=255)                                                                                        aaTAAAGA‑AAAAAAAAATTAACAC‑CCTTCAGAATATTTTTATAAAGGTAACCACCCACGCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                              <  2:55214/146‑1 (MQ=255)                                                                                          taaaGA‑AAAAAAAAATAAACAC‑CCGACAGGATTTTGTTAAAAAGGTAACGTTCCATTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                        <  1:814127/148‑1 (MQ=255)                                                                                            aaGG‑AAAAAAAAATAAACAC‑CCGTCAGAATTTTTTTAAAAATGTGCCGCTCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                        <  1:845656/148‑1 (MQ=255)                                                                                            aaGA‑AAAAAAAAATAAACAA‑CGTTAGGGTTTTTTTTATCAATGTGCCCACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                     <  2:582115/151‑1 (MQ=255)                                                                                             aGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCg                                                    >  2:713213/1‑151 (MQ=255)                                                                                                   gaagaaATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                               >  1:758999/1‑151 (MQ=255)                                                                                                     agaaATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                             >  1:1907/1‑151 (MQ=255)                                                                                                        aaatAACAC‑CCGTCAGAATTTTTTTAAAAATGACACCACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                          <  2:621723/148‑1 (MQ=255)                                                                                                           tAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       >  2:123781/1‑151 (MQ=255)                                                                                                            aacac‑cCGTAAGATTTTTTTTATAATTGGGCCGACCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                       <  1:544711/150‑1 (MQ=255)                                                                                                                   cgtagggATTTTTTTATAAATGTGACAACCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAAt                                <  2:777504/146‑1 (MQ=255)                                                                                                                           tatTGTAATAAATTTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTc                          <  1:303269/149‑1 (MQ=255)                                                                                                                               gTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtg                    >  1:601562/1‑151 (MQ=255)                                                                                                                                  aTCAATGTGACGATCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                   <  2:217345/149‑1 (MQ=255)                                                                                                                                    cAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc               >  2:85944/1‑151 (MQ=255)                                                                                                                                     attgtgACAACCCCCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg              <  1:692517/149‑1 (MQ=255)                                                                                                                                     aaTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg              >  1:733329/1‑151 (MQ=255)                                                                                                                                      aTGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCg              >  1:670666/1‑150 (MQ=255)                                                                                                                                        gtgCCGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGa             <  1:500105/149‑1 (MQ=255)                                                                                                                                          gACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATc         >  2:527139/1‑151 (MQ=255)                                                                                                                                            cGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg        <  1:570667/150‑1 (MQ=255)                                                                                                                                            cGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGatat          <  1:279697/148‑1 (MQ=255)                                                                                                                                             gATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc       >  2:831273/1‑150 (MQ=255)                                                                                                                                                 cACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg  >  2:927609/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                   CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGG‑AATTAAGG‑AAGAAGAAATTAACAA‑CCGTCAGGATATTGTTATCAATGTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTG  >  NC_000913/2421298‑2421589

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