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

CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421298‑2421592                                                                                                                                                    |                                                                                                                                                      ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAcccc                                                                                                                                                   <  1:1044591/150‑1 (MQ=255) ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAcccc                                                                                                                                                   >  2:713993/1‑150 (MQ=255) ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAcccc                                                                                                                                                   >  1:147776/1‑150 (MQ=255) ccTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTAcccc                                                                                                                                                   <  2:268979/150‑1 (MQ=255)   tAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCccgcc                                                                                                                                                >  2:623810/1‑151 (MQ=255)    aaGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACAATCCACTACCccgcc                                                                                                                                                <  1:663340/150‑1 (MQ=255)         aaaCAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATaa                                                                                                                                           >  2:694294/1‑150 (MQ=255)            caacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAgcg                                                                                                                                       >  2:1003662/1‑151 (MQ=255)             aacGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAgcgc                                                                                                                                      >  2:374745/1‑151 (MQ=255)                       cTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCg                                                                                                                             >  1:870261/1‑150 (MQ=255)                                cacacccccccGCTTTTTTTCGAAAAAAAACCCTTCTTTTGCGTATTTTTACAGAAATAAGGGAAAAAAAATTT‑AACACCCG‑CAGGAATTTTTTAACAATGGTGCCACCCCCCAACCCCGCCGATAAAGCGCGCACTAGCGaa                                                                                                                           <  2:689202/134‑1 (MQ=255)                                 acaACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATCTTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  >  1:515211/1‑151 (MQ=255)                                 acaACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAAc                                                                                                                  >  2:1085779/1‑151 (MQ=255)                                     cctgcgtttattatTCGTAAAAAAACCCGTTTTTGGCGTACTTTTATGGAAATAAAGAAAGAAAAAATA‑AACACCCGGAGGGATATTGTTAGCAAG‑GTTACGATCAACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATg                                                                                                              <  1:442934/143‑1 (MQ=255)                                      cTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATgc                                                                                                             >  2:750199/1‑151 (MQ=255)                                         ccGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCt                                                                                                          >  2:429527/1‑151 (MQ=255)                                               ttattaGTAAAAAGAACCTTCTTTTGCTTACGTTTATAGAAATAAAGAAAGAAAAAATA‑AACACCCGACAGAATATTTTAACCATG‑GTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGa                                                                                                        <  2:818088/147‑1 (MQ=255)                                                tattCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTAt                                                                                                   >  2:687472/1‑151 (MQ=255)                                                  tacgaaaaaaaacccTTCGTTTCCTTACTTTTAT‑GGAAATAAAGAAAAAAAAATTAAACACCGGTAAGTATATTGTTATAAAT‑GAGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTa                                                                                                    <  1:561222/135‑1 (MQ=255)                                                     gTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCAc                                                                                              >  1:538561/1‑151 (MQ=255)                                                       ataAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            >  1:712471/1‑151 (MQ=255)                                                       aaaaaaacccTTCTTTGGCGTCTGTTTATCGAAATTAAGAAAAAAAAAATA‑AACACCCGTAAGGATATTGTAATAAAT‑GTGCCCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGc                                                                                                <  1:179793/139‑1 (MQ=255)                                                       aaaaaGACCCTTCGTTTGCGTACTTTTATCGAAATAAAGAAAAAAAAAATA‑AACACCCGTCAGTATATTTTTATCATT‑GTGACGCCCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCa                                                                                               <  2:244205/146‑1 (MQ=255)                                                        aaaaaaccttttgTTTGCGTTTTTTTATAGAAATAAAGGAAAAAAAAATA‑ACCACCCGGCGGGTTTTTGTTATCATG‑GTGACGCCCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           <  1:362232/139‑1 (MQ=255)                                                         aaaGACCCTTCTTTTGTGTCTTTTTACCGAAATAAAGAAAAAAAAAATA‑AACCCCCGTGAGAATTTTTTTACAAAT‑GTGACAACCCCTTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCg                                                                                            <  2:748031/149‑1 (MQ=255)                                                         aaaGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GCGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGt                                                                                           >  2:930092/1‑150 (MQ=255)                                                           aGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAAc                                                                                        >  2:763737/1‑151 (MQ=255)                                                           aGAACCTTCGTTTGCGGACGTATATCGGAATTAAGGCAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCACTCCGATAAAGCGCGCAGTAGCTAAGATATCAACGATGAGCTGCATTATCGCACCGTAAc                                                                                        >  1:291667/1‑151 (MQ=255)                                                              aCCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGa                                                                                      >  1:131124/1‑150 (MQ=255)                                                                   cGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAc                                                                                 >  1:493108/1‑150 (MQ=255)                                                                        gcttACTTTTATCGAAATAAAGAAAAAAAAAATA‑AACACCCGCAGGGATATTTTTATAAAT‑GTGCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACattat                                                                            <  2:308031/147‑1 (MQ=255)                                                                         cGTATGTTTATCGAAATAAGGAAAAAAAAAATA‑AACACCCGTCAGGATTTTGTTATCAAT‑GTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          <  2:203906/151‑1 (MQ=255)                                                                           tACGTTTATCGGAATTAAGGAAGACGAAATT‑AACAAACGCAAGGATATTGTTATCTAT‑GGGACGATCCACTACCCAGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAAcatcat                                                                            >  1:696552/1‑147 (MQ=255)                                                                           tACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTc                                                                          >  1:443533/1‑149 (MQ=255)                                                                             cGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGc                                                                       >  2:882617/1‑150 (MQ=255)                                                                               tttATCGAAATAAAGAAAAAAAAAATA‑ACCACCCGTCAGTATATTTTTATCATT‑GTAACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCa                                                                      <  2:580305/149‑1 (MQ=255)                                                                                tttgcgaaaTAAAGGAAGAAAAAAAT‑ACCACCCTAAAGAATTTTTTTAAAAAT‑GAGACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATg                                                                    <  1:237343/143‑1 (MQ=255)                                                                                     ggAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 >  1:666345/1‑148 (MQ=255)                                                                                      gAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAgaga                                                              >  1:24778/1‑150 (MQ=255)                                                                                       aaTTAAAAAAAAAAAAATT‑AACACCCGTCAGGATATTTTTATCAGT‑GTGACTACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTag                                                                 <  1:749862/146‑1 (MQ=255)                                                                                       aaTAAAGAAAGAAAAAATT‑AACACCCGTCAGGTTTTTTTTATAAAT‑GTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGaa                                                             <  2:292756/150‑1 (MQ=255)                                                                                       aaTAAAGAAAAAATAAAAT‑ACCACCGGTAAGGATATTGTTATCGAG‑GTCACGACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCGCCGTAACGAAAAACATTATTCAGCATGTAGAGAAt                                                            <  1:317943/151‑1 (MQ=255)                                                                                         ttAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCACCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                           >  1:344366/1‑150 (MQ=255)                                                                                         taaagaaagaAAAAATA‑ACCACCCGGCAGGATTTTTTTATCAAG‑GTCACCACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTagag                                                               <  1:193486/144‑1 (MQ=255)                                                                                          tAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAc                                                         >  2:243477/1‑151 (MQ=255)                                                                                           aagaaagaAAAAATT‑AACACCCGACAGGATATTGTTATCAAT‑GTGACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                           <  2:505486/148‑1 (MQ=255)                                                                                           aaGAAAAAAAAAATT‑AACACCCGTCAGAATTTTGTTACCAAT‑GTGACAACCCACACCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATa                                                           <  2:68069/148‑1 (MQ=255)                                                                                           aaGAAAAAAAAAATA‑AACACCCGCCAGGATATTTTTATAAAT‑GTGACGACCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAAcc                                                        <  2:759887/151‑1 (MQ=255)                                                                                               aagaagaaATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    >  2:332449/1‑151 (MQ=255)                                                                                               aagaaaaaATC‑AACACCCGTCAGGATATTTTTAACAAG‑GTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGttt                                                    <  2:13078/151‑1 (MQ=255)                                                                                                 gaagaaATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTtct                                                  >  1:610497/1‑151 (MQ=255)                                                                                                   agaaATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTctct                                                >  2:873380/1‑151 (MQ=255)                                                                                                    gaaATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCtt                                               >  2:624576/1‑151 (MQ=255)                                                                                                       aaacaacaccCGTCAGTTTATTTTTATCAAT‑GTAACAACCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTg                                             <  1:332258/147‑1 (MQ=255)                                                                                                        tt‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCt                                           >  2:728026/1‑151 (MQ=255)                                                                                                           aacaccCGTAAGGATGTGGATAAAAAT‑GTGACGCCCCCCCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTg                                          <  1:379836/150‑1 (MQ=255)                                                                                                               acCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaa                                      >  1:1084393/1‑150 (MQ=255)                                                                                                               acCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGACAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGACTAACCGTTTCTCTTTGCTGGaca                                      >  1:472419/1‑148 (MQ=255)                                                                                                                ccttCAGAATTTTTTTACCAGT‑GAAACGCTCCTCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaaa                                    <  2:105389/148‑1 (MQ=255)                                                                                                                ccGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGaaaa                                     >  2:1084394/1‑150 (MQ=255)                                                                                                                                acaaaT‑GTAACGCCCCACCCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgtg                     <  1:256954/147‑1 (MQ=255)                                                                                                                                   att‑gtgACCACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGAtgt                      <  2:740959/144‑1 (MQ=255)                                                                                                                                   aaT‑GTGACAACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTc                  <  1:143007/150‑1 (MQ=255)                                                                                                                                    aaggtgACGATCCACTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGc                    <  1:777068/145‑1 (MQ=255)                                                                                                                                       gtgCCGACCCCCTCCCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGata              <  2:226126/151‑1 (MQ=255)                                                                                                                                          aCGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCg           >  2:126369/1‑151 (MQ=255)                                                                                                                                            gATCCACTACCCCTCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGACTTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGc          >  1:270495/1‑150 (MQ=255)                                                                                                                                                 aCTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTg     >  2:378001/1‑150 (MQ=255)                                                                                                                                                   tACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAc  >  2:842604/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCTAAGTAAAACAACGGAGAACCTGCAATGGCACAACCTGCCGCTATTATTCGTATAAAGAACCTTCGTTTGCGTACGTTTATCGGAATTAAGGAAGAAGAAATT‑AACAACCGTCAGGATATTGTTATCAAT‑GTGACGATCCACTACCCCGCCGATAAAGCGCGCACTAGCGAAGATATCAACGATGCGCTGAATTATCGCACCGTAACGAAAAACATTATTCAGCATGTAGAGAATAACCGTTTCTCTTTGCTGGAAAAATTAACTCAGGATGTGCTCGATATCGCACGTGAAC  >  NC_000913/2421298‑2421592

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