Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,596,873	Δ1 bp	10.8%	intergenic (‑136/+32)	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,873	0	C	.	10.8%	207.0 / 20.7	65	intergenic (‑136/+32)	ypfJ/purC	putative neutral zinc metallopeptidase/phosphoribosylaminoimidazole‑succinocarboxamide synthetase
Reads supporting (aligned to +/- strand):  ref base C (24/34);  new base . (4/3);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 4.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.90e-01

CTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCA  >  NC_000913/2596725‑2597018                                                                                                                                                     |                                                                                                                                                    cTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                                                                     <  2:119316/149‑1 (MQ=255) cTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGc                                                                                                                                                     <  4:450969/149‑1 (MQ=255)      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:500986/1‑149 (MQ=255)            aTAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACgg                                                                                                                                        >  3:52973/1‑149 (MQ=255)               gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCaa                                                                                                                                     >  1:336776/1‑149 (MQ=255)               gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGC‑ACGGCaa                                                                                                                                     <  2:133386/149‑1 (MQ=255)               gATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTG‑CAGCTGGC‑ACGGCaa                                                                                                                                     <  1:389076/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                                                                                                                           <  4:187318/138‑1 (MQ=255)                                    ttATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAAc                                                                                                                           >  3:187318/1‑138 (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                                                                                                    <  2:129463/149‑1 (MQ=255)                                                tCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTaca                                                                                                    >  1:291036/1‑149 (MQ=255)                                                              gagGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCggcgggc                                                                                      <  1:200020/149‑1 (MQ=255)                                                               agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAgatgat                                                                                                                               <  2:248231/107‑1 (MQ=255)                                                               agGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAgatgat                                                                                                                               >  1:248231/1‑107 (MQ=255)                                                                      ttACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                                                              <  1:355293/149‑1 (MQ=255)                                                                      ttACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCtt                                                                              >  4:328590/1‑149 (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                                                                                                          >  3:380662/1‑115 (MQ=255)                                                                            caGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGc                                                                                                          <  4:380662/115‑1 (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                                                                                                                 >  2:310450/1‑85 (MQ=255)                                                                                                   gaggaATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCa                                                                                                                 <  1:310450/85‑1 (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)                                                                                                                                                     |                                                                                                                                                    CTTGCCAACGCATAGATACCTCAACAATATATTCATTATAGAAATGATCGTAGGTGGTTCAGAGGAAGATTACCACAGGAAAGCGGGATGAAATGGCGGAGGAATGAGCGTTCAATGCATATATTCCGAGAGCGGTTGTCTTGCCGTGCCAGCT‑GC‑ACGGCAAGATGATGAACAGAAAAATCAGTCCAGCTGTACACCCAGGCGGCGGGCGACGGCTTCATAGGCTTCGATCAGGCCACCGAGGCTCTGGCGGAAACGGTCTTTGTCCATTTTCTCCAGCGTTTCTTTGTCCCA  >  NC_000913/2596725‑2597018

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