Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	422,537	T→G	6.0%	L7R (CTG→CGG)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	422,537	0	T	G	6.0%	229.8 / 7.0	67	L7R (CTG→CGG)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base G (2/2);  total (35/32)
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.26e-01

GAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAA  >  NC_000913/422389‑422672                                                                                                                                                     |                                                                                                                                        gaAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCt                                                                                                                                         >  4:274658/1‑149 (MQ=255)     tCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCg                                                                                                                                     >  3:251483/1‑149 (MQ=255)      ccTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCgg                                                                                                                                    >  3:83835/1‑149 (MQ=255)           tCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGccc                                                                                                                               >  3:130755/1‑149 (MQ=255)                 tATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATTGCACCTGCCGGTGCCCCCAttt                                                                                                                         >  3:572327/1‑149 (MQ=255)                    tCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCAATTTGtt                                                                                                                      <  2:324734/149‑1 (MQ=255)                      cGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaa                                                                                                                    >  4:207367/1‑149 (MQ=255)                      cGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTaa                                                                                                                    >  4:70334/1‑149 (MQ=255)                           gTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaa                                                                                                               >  3:530012/1‑149 (MQ=255)                               tCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaa                                                                                                           >  1:10143/1‑149 (MQ=255)                                ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaa                                                                                                          <  4:212320/149‑1 (MQ=255)                                ccGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAAcaaagcaaa                                                                                                          <  1:578090/149‑1 (MQ=255)                                  gCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGa                                                                                                        <  2:10143/149‑1 (MQ=255)                                      ccGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCaa                                                                                                    <  4:501838/149‑1 (MQ=255)                                       cGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc                                                                                                   <  4:554119/149‑1 (MQ=255)                                              ctccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAt                                                                                            <  4:540277/149‑1 (MQ=255)                                               tccCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAtt                                                                                           <  4:416657/149‑1 (MQ=255)                                                  cccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTaca                                                                                        >  1:471741/1‑149 (MQ=255)                                                  cccAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTaca                                                                                        >  2:566945/1‑149 (MQ=255)                                                       atatCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACtgtg                                                                                   <  3:70334/149‑1 (MQ=255)                                                        tatCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTgg                                                                                  <  1:110982/149‑1 (MQ=255)                                                         atCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGc                                                                                 <  1:448663/149‑1 (MQ=255)                                                          tCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCt                                                                                <  3:526790/149‑1 (MQ=255)                                                          tCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCt                                                                                <  4:351541/149‑1 (MQ=255)                                                             ttCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGAc                                                                             >  3:300021/1‑149 (MQ=255)                                                                aGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                                                                                                                         >  2:550624/1‑102 (MQ=255)                                                                aGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                                                                                                                         <  1:550624/102‑1 (MQ=255)                                                                 gatgatGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc                                                                                                   <  4:246529/123‑1 (MQ=255)                                                                 gatgatGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAAc                                                                                                   >  3:246529/1‑123 (MQ=255)                                                                    gatgaGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGaa                                                                      <  2:13672/149‑1 (MQ=255)                                                                          tGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATTGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCa                                                                >  3:436404/1‑149 (MQ=255)                                                                               cTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCa                                                           >  3:576894/1‑149 (MQ=255)                                                                               cTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCa                                                           >  3:194109/1‑149 (MQ=255)                                                                               cTGCATTAGGCTATCGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCa                                                           >  1:231078/1‑149 (MQ=255)                                                                                     tAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAt                                                     <  3:350242/149‑1 (MQ=255)                                                                                     tAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAt                                                     <  2:290537/149‑1 (MQ=255)                                                                                      aGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAtt                                                    >  3:522280/1‑149 (MQ=255)                                                                                                  ggTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACa                                        <  1:125004/149‑1 (MQ=255)                                                                                                      aTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAac                                    >  2:500742/1‑149 (MQ=255)                                                                                                      aTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAac                                    >  3:266557/1‑149 (MQ=255)                                                                                                         aGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTGCACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACg                                 <  2:487314/149‑1 (MQ=255)                                                                                                         aGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACg                                 >  2:393273/1‑149 (MQ=255)                                                                                                         aGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACg                                 >  1:487314/1‑149 (MQ=255)                                                                                                                catcaCAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACga                          >  2:34019/1‑149 (MQ=255)                                                                                                                 atcaCAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         <  1:9458/149‑1 (MQ=255)                                                                                                                    acaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaaga                       >  4:522958/1‑148 (MQ=255)                                                                                                                    acaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaaga                       <  3:522958/148‑1 (MQ=255)                                                                                                                      aGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAAt                    >  1:259289/1‑149 (MQ=255)                                                                                                                          gAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCa                <  1:302975/149‑1 (MQ=255)                                                                                                                          gAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCa                <  1:568088/149‑1 (MQ=255)                                                                                                                             ttATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTa             >  1:23802/1‑149 (MQ=255)                                                                                                                                tgatgTTAATTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGTACACGATAACGAAGAAATGTCAGTACCg          >  4:504074/1‑149 (MQ=255)                                                                                                                                tgatgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCg          >  4:21539/1‑149 (MQ=255)                                                                                                                                tgatgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCg          <  1:78031/149‑1 (MQ=255)                                                                                                                                  atgTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt        <  4:522280/149‑1 (MQ=255)                                                                                                                                    gTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAtt                                                                                           >  4:134292/1‑64 (MQ=255)                                                                                                                                    gTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCAtt                                                                                           <  3:134292/64‑1 (MQ=255)                                                                                                                                     ttAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCGTTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCa     >  1:551506/1‑149 (MQ=255)                                                                                                                                     ttAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCa     >  2:232120/1‑149 (MQ=255)                                                                                                                                      tAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCAt    <  2:23802/149‑1 (MQ=255)                                                                                                                                        aaTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCATCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATaa  <  2:198667/149‑1 (MQ=255)                                                                                                                                        aaTGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATaa  >  4:577040/1‑149 (MQ=255)                                                                                                                                         aatgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         >  2:42950/7‑125 (MQ=255)                                                                                                                                         aatgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         <  1:42950/119‑1 (MQ=255)                                                                                                                                          atgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt        >  3:529222/6‑141 (MQ=255)                                                                                                                                          atgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt        <  4:529222/136‑1 (MQ=255)                                                                                                                                                  aCCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         >  4:516501/1‑116 (MQ=255)                                                                                                                                                  aCCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         <  3:516501/116‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                        GAAATCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAA  >  NC_000913/422389‑422672

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