Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,020,821	A→G	6.1%	S94S (TCA→TCG)	ygfM →	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,020,821	0	A	G	6.1%	166.4 / 3.0	67	S94S (TCA→TCG)	ygfM	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (34/27);  new base G (2/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.58e-01

GCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGCCTGCC  >  NC_000913/3020675‑3020969                                                                                                                                                   |                                                                                                                                                     gCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACg                                                                                                                                                    >  1:422241/1‑149 (MQ=255)        ccTAACAGGATCTGGAACAGGACTGTGATGACTGGGATAACGGTGCACTGCGGATTGGCGCAAAGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGGAGCGCAGCATGATGCAGTCGGATTTGTTTAATCACACCATGtt                                                                                                                                             >  2:95361/1‑149 (MQ=255)         cTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTc                                                                                                                                            >  1:798338/1‑149 (MQ=255)         cTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTc                                                                                                                                            <  1:121047/149‑1 (MQ=255)          ttACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGGGAAGCCCGCGGTTTTGTTTACGCCCGCCATGTTCg                                                                                                                                           >  2:495978/1‑149 (MQ=255)              cggATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAAt                                                                                                                                       <  1:168225/148‑1 (MQ=255)              aGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCTCTAGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGAGAAGCCATCGGTTTTGTTTACTCACGACAAGTTCGTAAt                                                                                                                                       >  1:308360/1‑149 (MQ=255)               ggATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATc                                                                                                                                      <  2:489602/149‑1 (MQ=255)                gtacTGGAACTGGACTGGGGTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCa                                                                                                                                     <  2:564920/146‑1 (MQ=255)                gATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCa                                                                                                                                     >  2:929677/1‑149 (MQ=255)                  tCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGt                                                                                                                                   <  2:388953/149‑1 (MQ=255)                       aaCTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCATTCGAcc                                                                                                                              >  1:929588/1‑149 (MQ=255)                       aaCTGGACTGGGTTGACTGGGATAACGATGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTTAAGCCCTCGGTTTTGTTTACTCACGCCATGATCGTAATCAGTCGAcc                                                                                                                              >  1:147779/1‑149 (MQ=255)                        actggactggGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCGCGCCATGTTCGTAATCAGTCGACCa                                                                                                                             >  1:90068/1‑149 (MQ=255)                        actggactggGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCGCGCCATGTTCGTAATCAGTCGACCa                                                                                                                             >  1:93646/1‑149 (MQ=255)                              ctggGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATtggtg                                                                                                                       <  1:70114/149‑1 (MQ=255)                                gggTTGACTGGGATAACGGTGCACTGCGGCTTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGAACATTggtggt                                                                                                                     <  2:182348/149‑1 (MQ=255)                                 ggTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGgtggtg                                                                                                                    <  2:288859/149‑1 (MQ=255)                                  gTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGa                                                                                                                   >  2:79896/1‑149 (MQ=255)                                      aCTGGGACACCGGTGCACTGCGGATTGGCGCAATGGCTCGCTTGCAGCCACTGCGGGAGGCGCGATTTATTCCTGCCGCGCTGCGTGAAGCCCTCGGTTGTGTTTACTCACGCCATGTTCGTAATCAGGCGACCATTGGTGGTGAAATc                                                                                                               <  2:475236/149‑1 (MQ=255)                                        tGGGATAACGGTGCACTGCGGATGGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTGTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATcgc                                                                                                             <  2:759130/149‑1 (MQ=255)                                         gggATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATcgcc                                                                                                            <  1:84602/149‑1 (MQ=255)                                         gggATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATcgcc                                                                                                            <  1:460835/149‑1 (MQ=255)                                                   tGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAgg                                                                                                  <  2:207261/149‑1 (MQ=255)                                                        tGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAagag                                                                                             >  1:603679/1‑149 (MQ=255)                                                            gATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCTTGTTCGTAATCAGTCGCCCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCgg                                                                                         >  1:347529/1‑149 (MQ=255)                                                            gATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGCCCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCgg                                                                                         >  1:298912/1‑149 (MQ=255)                                                             aTTGGCACAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGt                                                                                        >  2:342244/1‑149 (MQ=255)                                                              ttGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGtg                                                                                       >  1:664369/1‑149 (MQ=255)                                                               tGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGtgc                                                                                      <  1:738801/149‑1 (MQ=255)                                                                 gcgcAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCCGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGtgctg                                                                                    >  1:341571/1‑149 (MQ=255)                                                                 gcgcAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGtgctg                                                                                    >  2:80289/1‑149 (MQ=255)                                                                 gcgcAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGAATGCCGTCTCCCATAGGTGGAGAAATAGCCGCACGCCAGGAAGAGTCGGttcgg                                                                                    >  2:946806/1‑145 (MQ=255)                                                                     aaTGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGAGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTc                                                                                <  1:641006/149‑1 (MQ=255)                                                                         tctcGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGt                                                                            >  2:591623/1‑149 (MQ=255)                                                                          gtctcTTGAAGGCACTGCGTGATGCGCGATTTATTCCTGCAGCGGTGCGTGAAGCCGTCGGTTTTGTTTACTCACGCCTTGTTCGTAAGCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTc                                                                           <  2:634069/145‑1 (MQ=255)                                                                          ctcGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTc                                                                           <  1:316340/149‑1 (MQ=255)                                                                          ctcGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGGGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGGCGGTGCTGCTTCCCGTc                                                                           <  1:224543/149‑1 (MQ=255)                                                                             gCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCAGTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTg                                                                        >  2:170186/1‑149 (MQ=255)                                                                              cTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGt                                                                       <  1:929677/149‑1 (MQ=255)                                                                                 gaaGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCGCGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTgg                                                                    <  2:90068/147‑1 (MQ=255)                                                                                 gCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCGCGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTgg                                                                    <  2:93646/149‑1 (MQ=255)                                                                                   aGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGATATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCa                                                                  >  2:1002396/1‑149 (MQ=255)                                                                                     ccACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGGTCGTAATCAGTCGGCCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCagg                                                                >  2:453576/1‑147 (MQ=255)                                                                                           cGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCCCGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGc                                                          <  2:877846/149‑1 (MQ=255)                                                                                                 gcgcGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACt                                                    >  2:103372/1‑149 (MQ=255)                                                                                                  cgcgATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCTGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAAcag                                                   >  1:108636/1‑147 (MQ=255)                                                                                                  cgcgATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTg                                                   >  1:472037/1‑149 (MQ=255)                                                                                                  cgcgATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGAAGCTGAACTg                                                   >  1:465608/1‑149 (MQ=255)                                                                                                  cgcgAATTATTCCAGCAGCGCCGCGTGAAGCCCGCGGGGTTGTTTACTCACGCCATGTTCGTAAACAGTCGACCATTGGGGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTg                                                   <  1:357070/149‑1 (MQ=255)                                                                                                        ttATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGGCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGttttt                                             <  1:597588/149‑1 (MQ=255)                                                                                                         tATTGCTGGTGCGCGTCGTGAAGCCCTCGGTTTTGTTTACTCACGCCTTGTTATGAATCAGTCGACCAGTGGTTGTGAAATCGCCGCCCGCCAGGTAGAGTCGGTGCTGCATCCCGTACTGTTGGCACTGGATTCTGAACTGGTTTTTg                                            <  1:843770/149‑1 (MQ=255)                                                                                                         tATTCCTGCAGCGCTGCGTGAAGCCCTCGGGTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTg                                            <  1:341758/149‑1 (MQ=255)                                                                                                         tATTACTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTAGTAATCAGTCGACCATTGGTGGTGAAATCGCAGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTg                                            <  2:971633/149‑1 (MQ=255)                                                                                                              cTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGCCCATTGGTGGTGAAATCGCCGCCCGCGAGGAAGAGTCGGTGCTGCGTCCCGTCCTGTTGGCACTGGATGGTGAACTGGTTTTTGGCAAc                                       >  2:169221/1‑149 (MQ=255)                                                                                                              cTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTCGCAAc                                       >  2:332090/1‑149 (MQ=255)                                                                                                              cTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCAAGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCAAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACGGGATGCTGAACGGGTTTTTGGCAAc                                       >  2:931943/1‑149 (MQ=255)                                                                                                                      cTGCGTGAAGCCCTCGGGTTTGTTTCCTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAAc                               <  2:959903/149‑1 (MQ=255)                                                                                                                           tGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGCAACGCTGt                          >  1:731645/1‑149 (MQ=255)                                                                                                                              aGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCaa                       <  2:603679/149‑1 (MQ=255)                                                                                                                                  cTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGa                   >  2:386140/1‑149 (MQ=255)                                                                                                                                         ttGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCGGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTAc            >  1:622113/1‑149 (MQ=255)                                                                                                                                           gTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCt          >  1:654219/1‑149 (MQ=255)                                                                                                                                              tACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTACTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGc       <  2:832701/149‑1 (MQ=255)                                                                                                                                               aCTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGcc      >  2:728779/1‑149 (MQ=255)                                                                                                                                               aCTCACGCCATGTTCGTAATCAGTCGACAATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTCTTGGCACTGGATCCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATAGAGGACTACCTGGcc      >  1:291587/1‑149 (MQ=255)                                                                                                                                                cTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGCCt     >  1:391258/1‑149 (MQ=255)                                                                                                                                                   aCGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGCCTGcc  >  2:623007/1‑149 (MQ=255)                                                                                                                                                   aCGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACGGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGCCTACCTGGCCTGcc  >  2:771176/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCATTTCCTTACAGGATCTGGAACTGGACTGGGTTGACTGGGATAACGGTGCACTGCGGATTGGCGCAATGTCTCGCTTGCAGCCACTGCGTGATGCGCGATTTATTCCTGCAGCGCTGCGTGAAGCCCTCGGTTTTGTTTACTCACGCCATGTTCGTAATCAGTCGACCATTGGTGGTGAAATCGCCGCCCGCCAGGAAGAGTCGGTGCTGCTTCCCGTCCTGTTGGCACTGGATGCTGAACTGGTTTTTGGCAACGGCGAAACGCTGTCAATCGAGGACTACCTGGCCTGCC  >  NC_000913/3020675‑3020969

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