Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,596,878:1	+G	9.4%	intergenic (‑141/+27)	ypfJ ← / ← purC	putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,596,878	1	.	G	9.4%	247.3 / 22.2	64	intergenic (‑141/+27)	ypfJ/purC	putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base . (23/35);  new base G (3/3);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGC  >  NC_000913/2596730‑2597023                                                                                                                                                      |                                                                                                                                                   cAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCt                                                                                                                                                     >  3:435404/1‑149 (MQ=255)  aaCGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑g                                                                                                                                                   >  3:19458/1‑149 (MQ=255)      cATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACg                                                                                                                                              <  1:48281/149‑1 (MQ=255)       aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACgg                                                                                                                                             >  3:52973/1‑149 (MQ=255)       aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACgg                                                                                                                                             >  3:500986/1‑149 (MQ=255)          gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCaa                                                                                                                                          >  1:336776/1‑149 (MQ=255)          gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGC‑ACGGCaa                                                                                                                                          <  1:389076/149‑1 (MQ=255)          gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGC‑ACGGCaa                                                                                                                                          <  2:133386/149‑1 (MQ=255)            tACCTCAACACTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGAGTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAga                                                                                                                                        <  4:500986/149‑1 (MQ=255)                tCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGTCGTGCCAGCT‑GC‑ACGGCAAgatgat                                                                                                                                    <  3:405189/149‑1 (MQ=255)                     aaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGCTTACCACAGGAAAGCGGGATGAAATGGCGGAGCAATGAGCGTTCAGTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCaa                                                                                                                                          <  4:430066/138‑1 (MQ=255)                     aaTATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCaa                                                                                                                                          >  3:430066/1‑138 (MQ=255)                        atatTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGaa                                                                                                                            >  2:200020/1‑149 (MQ=255)                          atTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGaaaa                                                                                                                          <  4:435404/149‑1 (MQ=255)                           ttCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGaaaaa                                                                                                                         >  2:355293/1‑149 (MQ=255)                             cATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATc                                                                                                                       >  3:176366/1‑149 (MQ=255)                              aTTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      <  1:266673/149‑1 (MQ=255)                               ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAAc                                                                                                                                >  3:187318/1‑138 (MQ=255)                               ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAAc                                                                                                                                <  4:187318/138‑1 (MQ=255)                               ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAg                                                                                                                     <  2:265373/149‑1 (MQ=255)                                 ataGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGAGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTc                                                                                                                   <  3:163024/149‑1 (MQ=255)                                  taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTcc                                                                                                                  >  1:50437/1‑149 (MQ=255)                                  taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGT‑ACGGCaa                                                                                                                                          <  1:473191/125‑1 (MQ=255)                                  taGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCT‑GGTACGGCaa                                                                                                                                          >  2:473191/1‑125 (MQ=255)                                      aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              <  4:542191/149‑1 (MQ=255)                                      aaTGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              <  1:156009/149‑1 (MQ=255)                                           tCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                                                                         >  1:291036/1‑149 (MQ=255)                                           tCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                                                                         <  2:129463/149‑1 (MQ=255)                                                         gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                                                                           <  1:200020/149‑1 (MQ=255)                                                          agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAgatgat                                                                                                                                    >  1:248231/1‑107 (MQ=255)                                                          agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAgatgat                                                                                                                                    <  2:248231/107‑1 (MQ=255)                                                                 ttACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                                                                   >  4:328590/1‑149 (MQ=255)                                                                 ttACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                                                                   <  1:355293/149‑1 (MQ=255)                                                                  tACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTc                                                                                  <  1:471178/149‑1 (MQ=255)                                                                     cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                               >  1:248710/1‑149 (MQ=255)                                                                     cacaGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATa                                                                               <  3:328590/149‑1 (MQ=255)                                                                       caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                               <  4:380662/115‑1 (MQ=255)                                                                       caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                               >  3:380662/1‑115 (MQ=255)                                                                           aaaGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTc                                                                         <  2:322477/149‑1 (MQ=255)                                                                               cGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATc                                                                     <  4:19458/149‑1 (MQ=255)                                                                                 ggATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGGTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCCTAGGCTTCGATCAg                                                                   <  2:219417/149‑1 (MQ=255)                                                                                  gATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAgg                                                                  >  4:534884/1‑149 (MQ=255)                                                                                       aaTGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAcc                                                             >  3:510200/1‑149 (MQ=255)                                                                                             ggaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGct                                                       >  1:330803/1‑149 (MQ=255)                                                                                              gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      <  1:310450/85‑1 (MQ=255)                                                                                              gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCa                                                                                                                      >  2:310450/1‑85 (MQ=255)                                                                                              gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGctc                                                      <  2:25149/149‑1 (MQ=255)                                                                                                         gTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCAGCGAGGCTCTGGCGGAAACg                                           <  4:52973/149‑1 (MQ=255)                                                                                                          ttCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACgg                                          <  4:233049/149‑1 (MQ=255)                                                                                                             aaTGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCt                                       <  3:172039/149‑1 (MQ=255)                                                                                                              attcaaatatTCCGAGGGCGGTTTTCTTGCCGGTCCAGCT‑GC‑ACGGCAAGATTATGAACAGAAAAATCAGTCCAGCTTTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCtt                                      <  3:451514/143‑1 (MQ=255)                                                                                                                   tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              >  3:157933/1‑72 (MQ=255)                                                                                                                   tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCt                                                                                                              <  4:157933/72‑1 (MQ=255)                                                                                                                   tatatTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTcc                                 <  3:463669/149‑1 (MQ=255)                                                                                                                       ttCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCAttt                             >  2:346982/1‑149 (MQ=255)                                                                                                                               ttcatcatCTTGCCGTG‑CAGCTGGC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGc                     >  4:528768/8‑149 (MQ=255)                                                                                                                                     gTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCt               <  2:291036/149‑1 (MQ=255)                                                                                                                                      tCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCtt              >  2:220188/1‑149 (MQ=255)                                                                                                                                       cTTGTCGCGCCAGCT‑GC‑GCGGCAAGATGCTGAACAGAAAAATCAGTCCACCTGTACACCCAGGCGGCGGGCGACGGCTCCCTAGGCTACGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCCTTTTCTCCAGCGTTTCttt             <  2:279324/149‑1 (MQ=255)                                                                                                                                       cTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCttt             <  4:61945/149‑1 (MQ=255)                                                                                                                                        ttGCCGTG‑CAGCTGGC‑ACGGCAAGAGGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTg            >  2:288370/1‑149 (MQ=255)                                                                                                                                         tGCCGTG‑CAGCTGGC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGt           >  3:9930/1‑149 (MQ=255)                                                                                                                                          gCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTc          <  4:176366/149‑1 (MQ=255)                                                                                                                                           ccGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTcc         <  2:330803/149‑1 (MQ=255)                                                                                                                                             gTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCca       >  1:328351/1‑149 (MQ=255)                                                                                                                                                  aGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAggc  <  3:534884/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   CAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCACAGGC  >  NC_000913/2596730‑2597023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑