Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,468,768	C→T	W362* (TGG→TGA)	chiA ←	endochitinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,468,768	0	C	T	96.9%	232.8 / NA	65	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (34/28);  minor base G (0/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGT  >  NC_000913/3468524‑3469013                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cccATGCGTCTGGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCg                                                                                                                                                                                                                                                  <  2:230599/251‑1 (MQ=255)    ataggtgttGAGCGACATAATTTGTATAATGTGCATTGGTATAATATACGCGTCTGGCCGCGCTGTGCGTTTTGCTGGCTGGGTAACCGTCCCCACTGATTACGCCGCCGGTATCCATGTCACAGGTAGTCGGGTTACCGTACTGACTCAATACCGCAGCGGTAGCTGTGAGCTCGGGACGCGATGGGTTATTGGAGTTTTCAGCGCTGGTCCCCCGTGGGCAATCTTTAGTCCCAACCGAGCAGGCGttg                                                                                                                                                                                                                                               <  2:420977/244‑2 (MQ=255)         tCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGt                                                                                                                                                                                                                                           <  2:615489/250‑1 (MQ=255)            tGATTGACAGAAGTTTTAAAATGTTCATCGTTATAATCTACGTTGTTGTCCGCGCTGTCCGCTTTGCTGGCCTGGAAACCATCCGTAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGTGTTCTTGTa                                                                                                                                                                                                                                          <  2:771736/248‑1 (MQ=255)                cGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACAtt                                                                                                                                                                                                                                    >  1:269408/1‑250 (MQ=255)                       aGTTTTATAATGTGCTTCTTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAGCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGa                                                                                                                                                                                                                            <  2:427785/251‑1 (MQ=255)                            tataATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAAt                                                                                                                                                                                                                        >  1:346622/1‑250 (MQ=255)                            tataATGTGAATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACCTTACCATTGAAAAt                                                                                                                                                                                                                        >  2:430873/1‑250 (MQ=255)                            tataAGGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGTCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTTAACATTACCATTGAAAAt                                                                                                                                                                                                                        <  2:789491/250‑1 (MQ=255)                             atagtgtgCTTCGTTATAATCGACGAGGCTGTCCGCGCGGTACCCTTTGCTGGCCTGGAAACCGTCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATa                                                                                                                                                                                                                       <  2:580832/246‑1 (MQ=255)                             ataATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATaa                                                                                                                                                                                                                      <  1:772667/251‑1 (MQ=255)                             ataATGTGCATCGTTATAATCTACCATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATa                                                                                                                                                                                                                       <  1:204005/250‑1 (MQ=255)                               aaTGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAAc                                                                                                                                                                                                                     >  1:219508/1‑250 (MQ=255)                                 tgtgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATg                                                                                                                                                                                                                  >  2:527920/1‑251 (MQ=255)                                  gtgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATgg                                                                                                                                                                                                                 >  1:720236/1‑251 (MQ=255)                                     cATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTa                                                                                                                                                                                                               >  1:403878/1‑250 (MQ=255)                                     cATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAAGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTaa                                                                                                                                                                                                              <  1:527920/251‑1 (MQ=255)                                             taATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACTCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCt                                                                                                                                                                                                      >  1:711500/1‑251 (MQ=255)                                              aaTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCt                                                                                                                                                                                                      >  2:93063/1‑250 (MQ=255)                                               accgactaTGCTGTCCGCGCGCTACGCTTTTCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTAc                                                                                                                                                                                                    <  2:720236/244‑1 (MQ=255)                                                         cTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGcc                                                                                                                                                                                          >  2:36480/1‑251 (MQ=255)                                                            cccGCTCTGTTCGCTTTGCGGGCATGGAACCCATCGACAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTAGTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTGAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCt                                                                                                                                                                                       <  2:511365/250‑1 (MQ=255)                                                              cgcgCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTcc                                                                                                                                                                                     >  2:58843/1‑251 (MQ=255)                                                                    gTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAttt                                                                                                                                                                                                          <  1:704782/224‑1 (MQ=255)                                                                    gTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAttt                                                                                                                                                                                                          >  2:704782/1‑224 (MQ=255)                                                                    gTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCtt                                                                                                                                                                               >  1:83736/1‑251 (MQ=255)                                                                                 gCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATc                                                                                                                                                                  <  1:548720/251‑1 (MQ=255)                                                                                     ggAAGCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCAGCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGt                                                                                                                                                                <  1:435746/249‑1 (MQ=255)                                                                                      gAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTaa                                                                                                                                                             <  1:517257/251‑1 (MQ=255)                                                                                           cATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCt                                                                                                                                                         >  2:233641/1‑250 (MQ=255)                                                                                              ccGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATcc                                                                                                                                                     >  2:624668/1‑251 (MQ=255)                                                                                                          cgccgcCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGAt                                                                                                                                         <  2:403878/251‑1 (MQ=255)                                                                                                              gccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCt                                                                                                                                     >  2:53538/1‑251 (MQ=255)                                                                                                               ccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTc                                                                                                                                    <  1:758644/251‑1 (MQ=255)                                                                                                                 gTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGc                                                                                                                                  >  2:539646/1‑251 (MQ=255)                                                                                                                       aaTTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTAt                                                                                                                            >  2:602993/1‑251 (MQ=255)                                                                                                                                 gTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTgg                                                                                                                  >  1:732941/1‑251 (MQ=255)                                                                                                                                       gggTTACCGTACTGACTCAATTCCTCCGCGGTAGCTGTACGCTCGAGACGCCATGGGTTACTGGAGTTTTCAGCGCTGGTACCCCGTGGTCAATCTTTAGACCCAACCCATAAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTATTCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAc                                                                                                            <  2:479848/251‑1 (MQ=255)                                                                                                                                       gggTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAc                                                                                                            >  2:299651/1‑251 (MQ=255)                                                                                                                                                 aCTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTcc                                                                                                  <  2:748041/251‑1 (MQ=255)                                                                                                                                                      cTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTc                                                                                                                                    <  2:100347/212‑1 (MQ=255)                                                                                                                                                      cTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTc                                                                                                                                    >  1:100347/1‑212 (MQ=255)                                                                                                                                                        cAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTg                                                                                           >  1:80752/1‑251 (MQ=255)                                                                                                                                                          aTTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTCTCCGGAGTTGGa                                                                                         <  1:539646/251‑1 (MQ=255)                                                                                                                                                               gCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGt                                                                                                                                                                <  2:499573/175‑1 (MQ=255)                                                                                                                                                               gCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGt                                                                                                                                                                >  1:499573/1‑175 (MQ=255)                                                                                                                                                                           gTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGc                                                                       <  2:269408/251‑1 (MQ=255)                                                                                                                                                                               gCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGTCGTTCTTGTAAACATTGCCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGcgccn                                                                   <  2:346622/251‑5 (MQ=255)                                                                                                                                                                                  cGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCAc                                                                     <  1:77258/246‑1 (MQ=255)                                                                                                                                                                                  cGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCAc                                                                     >  2:77258/1‑246 (MQ=255)                                                                                                                                                                                  cGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGcc                                                                 >  1:533945/1‑250 (MQ=255)                                                                                                                                                                                      aCGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCattat                                                            >  2:205217/1‑251 (MQ=255)                                                                                                                                                                                           aTGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATcc                                                       >  1:717780/1‑251 (MQ=255)                                                                                                                                                                                                  aTTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTa                                                 >  2:255453/1‑250 (MQ=255)                                                                                                                                                                                                                ccttgttacccgttggacAATCTTAAGACCCACCCCATCAGCCGTTCTTGAAAACATTCCAATTAAAAAAAACATGTTAAATTCCGCTCCCTCCCTGCCCGCTCCCCGCTTGCAATGATTAATCGTAAACGCTACCCGCGGGTGTGGAGGGCTCGGTGTTATCCGTAGTGGGGTTCCCGGGTTTTTCCGGAGTTGGA‑CAAGCGGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCa                                  <  1:104526/234‑1 (MQ=255)                                                                                                                                                                                                                       aCCCCGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAc                           >  1:322453/1‑251 (MQ=255)                                                                                                                                                                                                                          ncGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATcc                                                                                                                          >  2:352023/2‑154 (MQ=255)                                                                                                                                                                                                                          ccGTGGGCAATCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATcc                                                                                                                          <  1:352023/154‑1 (MQ=255)                                                                                                                                                                                                                                  aaTCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTg                >  2:720532/1‑251 (MQ=255)                                                                                                                                                                                                                                   aTCTTTAGACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTACTCGTTAACGCTATCCGCTGGCGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACTGGCGCCAGGGTTACCATCCTGGCTgc               >  2:778171/1‑250 (MQ=255)                                                                                                                                                                                                                                        tagcccaaacccaTCAGGCGTTCTTGAAAACATTACCATTAAAAATAACATGGTAAATTTCGCTCCCTCCCTGGCCGCTCCACGTTGGCAATGAGTAATCTTTAACGTTATCCGCTGGTGTTGATGGCTCGCTTTTATCCGGAGTTGGGGTCCCCGGTGTTTCGGGGGTGGGATCAC‑CCGGGGTCACAGCTCCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTc           <  1:779046/243‑1 (MQ=255)                                                                                                                                                                                                                                          gACCCAACCCATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgt        >  2:516137/1‑251 (MQ=255)                                                                                                                                                                                                                                          gACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAACTTTCGCTACCTTCCTGGCCGCTCCACGCTTGCACTGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CTGGAGCCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgt        >  2:289826/1‑251 (MQ=255)                                                                                                                                                                                                                                           ccccaacccaTCAGGTGTTTTTTAAAACATTCCCATAAAAAATAACATGATAATTTTCCCTCCCTCCCTGGCCGCTCCACGTTTGAAATGATTAATCGTTAACGCTCTCCGCGGGTTTTGATGGCTCGCTGTTATCCGGAGTGGGGGTCCCCGGTGTTTCCGGAGTTGAATCAC‑CCGGAGCCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgtt       <  2:322453/250‑1 (MQ=255)                                                                                                                                                                                                                                                acccaTCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccgt  >  1:435083/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       CCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAA‑CCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGT  >  NC_000913/3468524‑3469013

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