Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	423,295	G→T	5.3%	D260Y (GAC→TAC)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	423,295	0	G	T	5.3%	267.9 / 2.8	75	D260Y (GAC→TAC)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base G (32/39);  new base T (2/2);  total (34/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

CGATACCGAGGATTATCGCCATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTG  >  NC_000913/423150‑423441                                                                                                                                                  |                                                                                                                                                   cGATACCGAGGATTATCGCCATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGATGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGa                                                                                                                                                   <  4:95637/147‑1 (MQ=255)  gATACCGAGGATTATCGCCATGTCGATCCGCAGTTTGGCGGTGATGGGGCGGTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACa                                                                                                                                                 <  3:25335/148‑1 (MQ=255)             ttATCGCCATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGt                                                                                                                                      <  3:107662/148‑1 (MQ=255)                 cGCCATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGCCACAATACGCAGCAGCTGGGAATGCGGCTGGTGCCGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCa                                                                                                                                  >  6:214275/1‑148 (MQ=255)                   ccATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACg                                                                                                                                >  7:261958/1‑148 (MQ=255)                   ccATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAAACACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACg                                                                                                                                >  7:257403/1‑148 (MQ=255)                     aTGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACggg                                                                                                                              <  2:139506/148‑1 (MQ=255)                       gTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTTACAGGCATAATCGTGGCACGggtg                                                                                                                            >  7:187247/1‑148 (MQ=255)                       gTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTTACAGGCATAATCGTGGCACGggtg                                                                                                                            >  7:182501/1‑148 (MQ=255)                       gTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGggtg                                                                                                                            >  5:160862/1‑148 (MQ=255)                       gTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGggtg                                                                                                                            >  4:147761/1‑148 (MQ=255)                        tCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGggtg                                                                                                                            <  7:164610/147‑1 (MQ=255)                          gATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGgtggtg                                                                                                                         <  5:301916/148‑1 (MQ=255)                          gATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGgtggtg                                                                                                                         <  4:177679/148‑1 (MQ=255)                           cgttctaGTTTGGGGGGGTTGGGGCGGTGGTGCGGTTGGGCCACAATTCGGAGGAGGGGGGTATGGGGCTGGGGCTGGCCGGCGTGGTTTACCACAGGGGCGAATCCCATGCCTTGTTTGACAGGCATAATCGTGGCACGGgtggtg                                                                                                                         <  5:7454/141‑1 (MQ=255)                           aTCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGc                                                                                                                        >  8:304400/1‑148 (MQ=255)                             ccGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTTGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCt                                                                                                                       <  6:160862/147‑1 (MQ=255)                               gCAGTTTGGCGGTGAGGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCTCATGCCTGGTTTGACAGGCATAATCGTGGCACTGGTGGTGCTTGt                                                                                                                    <  2:56722/148‑1 (MQ=255)                                cAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCTCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTc                                                                                                                   <  8:261958/148‑1 (MQ=255)                                cAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTc                                                                                                                   <  8:257403/148‑1 (MQ=255)                                   tttGGCGGTGATGGGGCGTTGCTGCGTTTGCGCCACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACgg                                                                                                                               >  1:86760/1‑133 (MQ=255)                                   tttGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACgg                                                                                                                               <  2:86760/133‑1 (MQ=255)                                   tttGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTcac                                                                                                                 >  1:311491/1‑147 (MQ=255)                                    ttGCCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCGGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACaa                                                                                                               >  5:95443/1‑148 (MQ=255)                                       gCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAAcc                                                                                                             <  1:128453/147‑1 (MQ=255)                                           tGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGaa                                                                                                        >  6:111299/1‑148 (MQ=255)                                             aTGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAAt                                                                                                       <  5:304877/147‑1 (MQ=255)                                               ggggCGTTGCTGCGTTTGCGACACAATACGCAGCCGCTGGGAATGCGGATGGTGCTTGACGGCGTGTTTAACCACAGTTGCGATGCCCATGCCTGGTGTGACAGGCATAATCGTGGCACGGGTGGTGCTGGTCACCACCCCGAATCGc                                                                                                    >  7:232251/1‑148 (MQ=255)                                                        ctgcGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGcgcg                                                                                           <  7:304400/148‑1 (MQ=255)                                                         tgcGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGa                                                                                          >  4:139184/1‑148 (MQ=255)                                                                   acacAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCg                                                                                >  2:151617/1‑148 (MQ=255)                                                                    cacaATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGt                                                                               >  8:176042/1‑148 (MQ=255)                                                                    cacaATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGt                                                                               <  5:245955/148‑1 (MQ=255)                                                                     acaATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGtt                                                                              <  7:176042/148‑1 (MQ=255)                                                                      caATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGttt                                                                             <  8:164503/148‑1 (MQ=255)                                                                      caATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGttt                                                                             <  3:270591/148‑1 (MQ=255)                                                                        aTACGCAGCAGCGGGGAATGCGGCTTGCGCTGGCCTGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGGGGCACGGGTGGTGCGTGTCACAACCCCGAAGCGCCCTGGCGCGACTGGTACTCGTTTAg                                                                           <  8:232251/148‑1 (MQ=255)                                                                             cagcagCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGtgatg                                                                      <  1:242653/148‑1 (MQ=255)                                                                              agcagcTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGtgatga                                                                     <  1:151617/148‑1 (MQ=255)                                                                                cagcTGGGAATGCGGCTGGTGCTGGACGGCGTCTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGatgatg                                                                   >  1:317101/1‑148 (MQ=255)                                                                                     gggAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACg                                                              >  3:33919/1‑148 (MQ=255)                                                                                      ggAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCCCgg                                                             >  1:247892/1‑148 (MQ=255)                                                                                        aaTTCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCGATGCCTGGTTTGACAGGCATAATCGGGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGgcg                                                           <  4:218856/148‑1 (MQ=255)                                                                                            cGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTc                                                        >  4:20364/1‑147 (MQ=255)                                                                                                 ggTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGgc                                                                                              <  8:108579/104‑1 (MQ=255)                                                                                                 ggTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGgc                                                                                              >  7:108579/1‑104 (MQ=255)                                                                                                   tGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTa                                                                            >  8:92921/1‑120 (MQ=255)                                                                                                   tGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTa                                                                            <  7:92921/120‑1 (MQ=255)                                                                                                     cTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGTTACTCGTTTAGTGATGATGGCACGGCGCTCGACtggcttg                                              <  6:95443/148‑1 (MQ=255)                                                                                                       ggACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGAc                                                                                         >  4:25191/1‑103 (MQ=255)                                                                                                       ggACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGAc                                                                                         <  3:25191/103‑1 (MQ=255)                                                                                                       ggACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACtggcttggc                                            >  8:80638/1‑148 (MQ=255)                                                                                                        gACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTTACAGGCATAATCGTGGCACGGGTGTTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACtggcttggct                                           <  8:187247/148‑1 (MQ=255)                                                                                                        gACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTTACAGGCATAATCGTGGCACGGGTGTTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACtggcttggct                                           <  8:182501/148‑1 (MQ=255)                                                                                                         aCGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTa                                          <  5:214275/148‑1 (MQ=255)                                                                                                              gtgtTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGc                                       >  7:294699/1‑146 (MQ=255)                                                                                                                  ttAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCtt                                 >  7:208881/1‑148 (MQ=255)                                                                                                                    aaCCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTg                                <  3:20364/147‑1 (MQ=255)                                                                                                                     aCCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGc                               <  5:66178/147‑1 (MQ=255)                                                                                                                        acaGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCg                             >  8:274550/1‑146 (MQ=255)                                                                                                                             ggCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTg                       >  2:201169/1‑147 (MQ=255)                                                                                                                              gCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGa                     <  3:147761/148‑1 (MQ=255)                                                                                                                                gATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGcgcg                                                                                           >  1:76193/1‑76 (MQ=255)                                                                                                                                gATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGcgcg                                                                                           <  2:76193/76‑1 (MQ=255)                                                                                                                                         gCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGaa                           <  2:166472/131‑1 (MQ=255)                                                                                                                                         gCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGaa                           >  1:166472/1‑131 (MQ=255)                                                                                                                                         gCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTAt                 >  5:72942/1‑141 (MQ=255)                                                                                                                                         gCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTAt                 <  6:72942/141‑1 (MQ=255)                                                                                                                                         gCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCg           >  7:274021/1‑147 (MQ=255)                                                                                                                                          ccTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCgg          <  2:311491/147‑1 (MQ=255)                                                                                                                                             ggTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAACTTGCCGAAGCTGGATTATCAGTCGGAAAg      >  7:286459/1‑148 (MQ=255)                                                                                                                                              gTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAg      <  8:294699/147‑1 (MQ=255)                                                                                                                                                ttGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACtggc                                                 <  5:82190/102‑1 (MQ=255)                                                                                                                                                ttGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACtggc                                                 >  6:82190/1‑102 (MQ=255)                                                                                                                                                 tGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCt   <  3:139184/147‑1 (MQ=255)                                                                                                                                                  gACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTg  <  5:111299/147‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CGATACCGAGGATTATCGCCATGTCGATCCGCAGTTTGGCGGTGATGGGGCGTTGCTGCGTTTGCGACACAATACGCAGCAGCTGGGAATGCGGCTGGTGCTGGACGGCGTGTTTAACCACAGTGGCGATTCCCATGCCTGGTTTGACAGGCATAATCGTGGCACGGGTGGTGCTTGTCACAACCCCGAATCGCCCTGGCGCGACTGGTACTCGTTTAGTGATGATGGCACGGCGCTCGACTGGCTTGGCTATGCCAGCTTGCCGAAGCTGGATTATCAGTCGGAAAGTCTG  >  NC_000913/423150‑423441

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