Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	422,540	C→A	6.0%	P8Q (CCG→CAG)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	422,540	0	C	A	6.0%	234.9 / 5.0	67	P8Q (CCG→CAG)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base A (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 = 5.69e-01

TCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAA  >  NC_000913/422393‑422672                                                                                                                                                    |                                                                                                                                     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                                                                                                                         <  1:550624/102‑1 (MQ=255)                                                            aGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCAttt                                                                                                                         >  2:550624/1‑102 (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                                                     <  2:290537/149‑1 (MQ=255)                                                                                 tAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAt                                                     <  3:350242/149‑1 (MQ=255)                                                                                  aGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCAtt                                                    >  3:522280/1‑149 (MQ=255)                                                                                              ggTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACa                                        <  1:125004/149‑1 (MQ=255)                                                                                                  aTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAac                                    >  3:266557/1‑149 (MQ=255)                                                                                                  aTCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAac                                    >  2:500742/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                       <  3:522958/148‑1 (MQ=255)                                                                                                                acaGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaaga                       >  4:522958/1‑148 (MQ=255)                                                                                                                  aGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAAt                    >  1:259289/1‑149 (MQ=255)                                                                                                                      gAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCa                <  1:568088/149‑1 (MQ=255)                                                                                                                      gAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCa                <  1:302975/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                         <  1:42950/119‑1 (MQ=255)                                                                                                                                     aatgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         >  2:42950/7‑125 (MQ=255)                                                                                                                                      atgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt        <  4:529222/136‑1 (MQ=255)                                                                                                                                      atgggggCACCGGCAGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGAt        >  3:529222/6‑141 (MQ=255)                                                                                                                                              aCCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         >  4:516501/1‑116 (MQ=255)                                                                                                                                              aCCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACgaa                         <  3:516501/116‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                     TCCTGATCCGGCTATTCCGAAAGTTATCCGCCCCCGTCCTCCTCCCCCAAATATCCTTCAGATGATGAGTGATCCTGCATTAGGCTATGGCAAGGTGATCAGATTTTCATCACAGGGGAATTATGATGTTAAATGCATGGCACCTGCCGGTGCCCCCATTTGTTAAACAAAGCAAAGATCAACTGCTCATTACACTGTGGCTGACGGGCGAAGACCCACCGCAGCGCATTATGCTGCGTACAGAACACGATAACGAAGAAATGTCAGTACCGATGCATAA  >  NC_000913/422393‑422672

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