Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,017,074	A→G	5.7%	L409P (CTG→CCG)	ycbZ ←	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,017,074	0	A	G	5.7%	202.0 / 2.9	70	L409P (CTG→CCG)	ycbZ	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (40/26);  new base G (2/2);  total (42/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

ACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGATGACCCGGAAA  >  NC_000913/1016933‑1017222                                                                                                                                              |                                                                                                                                                     aCCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGc                                                                                                                                               >  1:682235/1‑149 (MQ=255)   cGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTGCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCtt                                                                                                                                             <  2:1114362/149‑1 (MQ=255)     gCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGc                                                                                                                                           <  1:683916/149‑1 (MQ=255)          tCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGat                                                                                                                                      <  2:44194/149‑1 (MQ=255)          tCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGat                                                                                                                                      <  2:455362/149‑1 (MQ=255)               gAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCTCCTTCACTGTATGACTGCTCAAATGTCAGCGAAGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATGATGAACGCTTGCATGatcatca                                                                                                                                 >  2:301585/1‑149 (MQ=255)                   gcgcTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTTCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGCGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGCACGCTTGCATGATCATCATccc                                                                                                                             >  1:810388/1‑149 (MQ=255)                       ttATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTc                                                                                                                         <  1:2172/149‑1 (MQ=255)                         aTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCGGGAACGCTTGCATGATCATCATCCCTTTCGc                                                                                                                       >  1:161163/1‑149 (MQ=255)                         aTCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCGGGAACGCTTGCATGATCATCATCCCTTTCGc                                                                                                                       >  1:156531/1‑149 (MQ=255)                                   agagTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTg                                                                                                             >  1:1031412/1‑149 (MQ=255)                                       ttCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTgccgc                                                                                                         <  2:632893/149‑1 (MQ=255)                                            ccATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGc                                                                                                    >  1:1171253/1‑149 (MQ=255)                                                    gCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCttt                                                                                            <  1:529597/149‑1 (MQ=255)                                                      aCTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCAGACATCAGGCACGCTTGCATGCTCTTCATCCCTTTCGCATGGATATTGCCGCCAAGCTACGCGTTGc                                                                                          >  2:666155/1‑149 (MQ=255)                                                       cTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCg                                                                                         >  1:359839/1‑149 (MQ=255)                                                          tctcCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTc                                                                                      <  1:1075875/149‑1 (MQ=255)                                                          tctcCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTc                                                                                      >  1:438548/1‑149 (MQ=255)                                                          tatcCTTCAACTGCACTGTATGACTGATCAAATGTCAGCGACGCTGAGAAGGGGATCTGTTGCTCAAGGTGTAGTTCCGACATCAGGAACGCATGCATGATCATCCTCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTc                                                                                      <  2:785588/147‑1 (MQ=255)                                                           ctcCATCAACTTGACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCg                                                                                     <  1:551127/149‑1 (MQ=255)                                                              cATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCCTCCCTTTCGCATGGATATTGCCTCCAAGCTACGCTATGCGTTCGATg                                                                                  >  1:26378/1‑149 (MQ=255)                                                              cATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATg                                                                                  >  1:429044/1‑149 (MQ=255)                                                                tCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTTCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCTGCTTTGCGTTCGATGTc                                                                                <  2:235498/149‑1 (MQ=255)                                                                   aCTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGt                                                                             <  2:40107/149‑1 (MQ=255)                                                                    cTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTg                                                                            >  1:1175177/1‑149 (MQ=255)                                                                    cTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAAAGCTTGCATGATCATCATCCCTTTCGCATGGATATTTCCGCCAAGCTCCGCTTTGCGTTCTATGTCGGTg                                                                            >  1:862690/1‑149 (MQ=255)                                                                         cTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTc                                                                       <  1:333169/149‑1 (MQ=255)                                                                           gTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACTCTTGCATTATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCAc                                                                     <  2:860680/149‑1 (MQ=255)                                                                            tATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGGGTTCGATGTCGGGGAATTCAcc                                                                    >  1:1186794/1‑149 (MQ=255)                                                                              tGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCAt                                                                  >  2:409193/1‑149 (MQ=255)                                                                               gACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATc                                                                 <  1:481758/149‑1 (MQ=255)                                                                                  tGCTCAAATGTCAGCGATGCTGAGCAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCCTGATCATCAGCCCGTTCGCATGGATATTGCCGCCACGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCCTCGcc                                                              >  2:889856/1‑149 (MQ=255)                                                                                    cTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCaa                                                            <  1:828304/149‑1 (MQ=255)                                                                                     tCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAat                                                           <  1:170710/149‑1 (MQ=255)                                                                                      cAAATGTCAGCGATGCTGAGTAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCTTCATCCCTTTCGCAAGGATATTGCCGCCAAGCTACGGTGTGCGTGCGATGTCTGTGAATTCCCCTTCGCCAata                                                          >  2:1017563/1‑149 (MQ=255)                                                                                       aaaTGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAatat                                                         >  1:717722/1‑149 (MQ=255)                                                                                             cAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGTGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACaa                                                   <  1:1003322/149‑1 (MQ=255)                                                                                              aGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAAc                                                  >  2:388814/1‑149 (MQ=255)                                                                                              aGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCACCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATTCACaag                                                  >  2:687493/1‑148 (MQ=255)                                                                                              aGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATCGATATTACCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCCTCGCCAATATGCACAAc                                                  >  2:808674/1‑149 (MQ=255)                                                                                                      tGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTaa                                          >  1:181906/1‑149 (MQ=255)                                                                                                       tagaAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAAt                                         <  2:826881/148‑1 (MQ=255)                                                                                                       gagaAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGAGTTCGATGTCTGTGAATTCACCATCGCCAATATGCAAAACGCAGCTAAt                                         >  2:888126/1‑149 (MQ=255)                                                                                                        agaAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGACCGCTTGCATGATCATCATCCCTTTCGCATGGAGATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGGGAATTCCCCAGCGCCAATATGCACCCCGCAGCTAATg                                        >  2:1107848/1‑149 (MQ=255)                                                                                                        agaAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATg                                        >  1:638555/1‑149 (MQ=255)                                                                                                        agaAGGGGATCTGGTGCTCAAGCTGAAGATCCGACATCAGGAACGCGTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCCAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATg                                        <  2:163538/149‑1 (MQ=255)                                                                                                         gaAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCGGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGc                                       <  2:161163/149‑1 (MQ=255)                                                                                                         gaAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCGGGAACGCTTGCATGATCATCATCCCTTTCGCAGGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGc                                       <  2:156531/149‑1 (MQ=255)                                                                                                         gaAGGGGATCTGTTGCTCAAGCTGAAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGc                                       <  1:311479/149‑1 (MQ=255)                                                                                                           aGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCga                                     <  2:269096/149‑1 (MQ=255)                                                                                                           aGGGGAGCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGTGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCga                                     <  1:1017605/149‑1 (MQ=255)                                                                                                            ggggATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCgag                                    <  1:28579/149‑1 (MQ=255)                                                                                                             gggATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCgaga                                   >  1:613609/1‑149 (MQ=255)                                                                                                             gggATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCgaga                                   >  1:1058673/1‑149 (MQ=255)                                                                                                              ggATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGaa                                  >  1:1112133/1‑149 (MQ=255)                                                                                                               gATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGAAGCTAATGCGAGAAg                                 >  2:574477/1‑149 (MQ=255)                                                                                                               gATCTGTTGCTCAAGCAGTAGTTCCGACATCAGGAACGCTTGCATGATCATAATCCCATACGCATGGATATTGCCGCCAACATCCGCTTTGCGATCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAg                                 >  2:1183986/1‑149 (MQ=255)                                                                                                                 tCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGt                               >  1:39741/1‑149 (MQ=255)                                                                                                                   tGTTGCTCAAGCTGTAGTTCCGACAGCAGGAACGCTTGCATGATCATCCTCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTc                             <  1:867410/149‑1 (MQ=255)                                                                                                                     ttGCTCAAGCTGTAGTTCCGACAGCAGGATCGCTTGCATGATCATCATCCGTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGGTCGACGTCGGGGAATTCACCAGCGCCTCGAGGCACAACGCAGCTCATGCGAGAAGGTTCGc                           >  1:766651/1‑149 (MQ=255)                                                                                                                       gCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCa                         >  1:488092/1‑149 (MQ=255)                                                                                                                         tCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATCCGAGAAGGTTCGCCaaa                       >  1:238818/1‑149 (MQ=255)                                                                                                                             gCTGTAGTTCCGACATCAGGAACGCTGGCATGATCATCATCCCTTGCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGCTGTCGGTGAATTCACCATCGCAAATGTGCCCACCGCCTCGAATGCGAGAAGGTGCGCCAAAcgcg                   >  2:573524/1‑149 (MQ=255)                                                                                                                              cTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATAGCCGCCAAGCTCCGCTGTGCGTTCGATGTCGGTGAATTCACAATCGCCAATATGCCCATCGCAGCTAATGCGAGAAGGTTCGCCCAAAgccc                  >  2:713370/1‑147 (MQ=255)                                                                                                                                gTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGt                >  1:298176/1‑149 (MQ=255)                                                                                                                                    ttCCGAGAACAGGAACGCTTGCTTGATCATCATCCCTTTAGCATTGATAATTCCTCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGAt            <  1:799795/149‑1 (MQ=255)                                                                                                                                       cGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGTGCGTGGATGAc         >  2:824892/1‑149 (MQ=255)                                                                                                                                       cGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTGCGATGTCGGTGAATTCCCCATCGCCAATATTCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGATGAc         >  2:1048912/1‑149 (MQ=255)                                                                                                                                        gACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGATGAcc        >  1:184209/1‑149 (MQ=255)                                                                                                                                          cATCAGGCACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCGTTGCGATCGATGTCGGTGAATTCACCATCGCCAATATGCACACCGCAGCGAAGGCGAGCAGGTTCGCCAAAAGCGCGGGGATgtcgcg      >  2:1141078/1‑144 (MQ=255)                                                                                                                                              aGGAACGCTTGTATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGATGACCCGGaaa  >  1:843018/1‑149 (MQ=255)                                                                                                                                              aGGAACGCTTGAATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGCTGATTTCACCATCGCCAATAGGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGATGACCCGGaaa  >  2:887111/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                     ACCGGCACATCGGCGAGGGCGCTTATCAGGGCGCAGAGTTCAGCCATCGAGGCACTATCTCCATCAACTTCACTGTATGACTGCTCAAATGTCAGCGATGCTGAGAAGGGGATCTGTTGCTCAAGCTGTAGTTCCGACATCAGGAACGCTTGCATGATCATCATCCCTTTCGCATGGATATTGCCGCCAAGCTCCGCTTTGCGTTCGATGTCGGTGAATTCACCATCGCCAATATGCACAACGCAGCTAATGCGAGAAGGTTCGCCAAAAGCGCGTGGATGACCCGGAAA  >  NC_000913/1016933‑1017222

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