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	97.6%	133.6 / NA	42	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (21/19);  minor base G (0/1);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑CCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTG  >  NC_000913/3468528‑3469009                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    ggCGTCTTGATCGACAGAAGGTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGCCGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCt                                                                                                                                                                                                                                           <  2:340083/250‑1 (MQ=255)  gCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTTCGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCtt                                                                                                                                                                                                                                          <  1:314506/251‑1 (MQ=255)    gTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGt                                                                                                                                                                                                                                        <  1:201764/251‑1 (MQ=255)         gATCGACAGAAGTTTTATAATGTGCATCGTTATAATCCACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTCGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAAc                                                                                                                                                                                                                                    >  2:96659/1‑250 (MQ=255)                   aGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTAGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGa                                                                                                                                                                                                                         <  1:347054/251‑1 (MQ=255)                   aGTTTTATAATGTGCATCGTTATAATCTACGATGCTGCCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGa                                                                                                                                                                                                                         <  2:452843/251‑1 (MQ=255)                      tttATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGaaa                                                                                                                                                                                                                       >  1:66280/1‑250 (MQ=255)                              gtgCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATgg                                                                                                                                                                                                              >  1:462711/1‑251 (MQ=255)                                     gTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAttt                                                                                                                                                                                                       >  1:265308/1‑251 (MQ=255)                                          aaTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCt                                                                                                                                                                                                   >  1:11923/1‑250 (MQ=255)                                           aTCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTAc                                                                                                                                                                                                 >  1:362495/1‑251 (MQ=255)                                              cacgattgtGTCCGCGCTGTACGCTTTGCTGGCCTGGAACCCACCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTTGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTGGTAAACATTTCCATTGAAAATAACATGGTAAATTTCGCTACCt‑t                                                                                                                                                                                             <  2:350473/243‑1 (MQ=255)                                               aCGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCt                                                                                                                                                                                               <  1:13194/249‑1 (MQ=255)                                                         ccGCGCTGTACGCTTTGCTGGCCTGGAACCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTCGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTCTTGGAGTTTTCAGCGCTGGTACCCAGTGGGCACTCTTTAGACCCAACCCA‑TCAGGCGTTCTTTTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTc                                                                                                                                                                                  >  2:372915/1‑251 (MQ=255)                                                         ccGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCt                                                                                                                                                                                   <  1:389486/250‑1 (MQ=255)                                                            cgcTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGGACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTTGTCCCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCAc                                                                                                                                                                               <  1:318100/251‑1 (MQ=255)                                                                 gacggtgtGCTGACCTGTAAACAATCCGAAACAGTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGAATCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTTCCT‑TCCTGGCCGCTCCACGCtt                                                                                                                                                                           <  2:271353/243‑1 (MQ=255)                                                                  aCGCTTGGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATGCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCGACCT‑TCCTGGCCGCTCCACGCTTg                                                                                                                                                                          <  2:272543/250‑1 (MQ=255)                                                                       ttGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAAt                                                                                                                                                                      <  1:431471/249‑1 (MQ=255)                                                                                                        ccgccgTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATgg                                                                                                                                   >  2:54651/1‑251 (MQ=255)                                                                                                                  cAATTTCACAGGTAGACGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGtt                                                                                                                          >  2:150079/1‑250 (MQ=255)                                                                                                                               aGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGCAAATAACATGGTAAATTACGCTACCT‑TCCTGGCCGCTCCACGCTTGCACTGAGTAATAGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCAGGCGTTgggg                                                                                                            >  2:228458/1‑251 (MQ=255)                                                                                                                                    ggTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGGTGATGGCTCGCTGTTATCCGGAGTTGGGGTCAcc                                                                                                       <  2:1766/251‑1 (MQ=255)                                                                                                                                     gTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCg                                                                                                      >  2:91973/1‑251 (MQ=255)                                                                                                                                           gTACTGACTCAATGCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGTTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGttt                                                                                                <  2:229313/251‑1 (MQ=255)                                                                                                                                                             aGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑g                                                                             <  1:54651/251‑1 (MQ=255)                                                                                                                                                                 gTAGCTATACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGt                                                                         >  1:364702/1‑251 (MQ=255)                                                                                                                                                                   aGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTc                                                                        >  1:450398/1‑250 (MQ=255)                                                                                                                                                                     cTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGACAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGAGCAACC‑GGAGTcac                                                                      >  2:115236/1‑250 (MQ=255)                                                                                                                                                                       gTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTGTCCGGAGTTGGATCAACC‑GGAGTcaca                                                                     <  2:265308/249‑1 (MQ=255)                                                                                                                                                                            cTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTATCCCGAGGGCAATCTTTAGACCCAACTCA‑TCAAGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTAGCCGCTCCACGCTTGCAATGAGAAATCGTTACCGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGGGTTGGATCAACC‑GGAGTCACAGCACCGc                                                              >  1:215887/1‑251 (MQ=255)                                                                                                                                                                                       aTGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCAAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATcc                                                   >  2:459005/1‑251 (MQ=255)                                                                                                                                                                                                 ggattgtccggcgcTGGGCCCCCGTGGACAATCTTTAG‑CCCAACCCAAGCAGGCTTTCTTGAAAACATTACCGTTGAAAATAACAGGGTAAATTTCGCT‑CCTGTCCTGCCCCCTCCCCGCTGGCAATGAGTAATCTTAAACGCTATCCGTGGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTCCCGGAGTTG‑AACAACCCGGATTCACAGCACCGCCATTATTATCCGGCTTAACg                                          <  1:163325/240‑1 (MQ=255)                                                                                                                                                                                                  gagatttcagcacatgTACCCCAAGGGAAATCTTTAGAACCACCCCC‑TCAGGCGATCTTGTAAACCTAACCATTTAAACTAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTGGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGCGGCACAGCACCGCCATTCTTATACGGCttcaaggg                                        >  2:418745/16‑245 (MQ=255)                                                                                                                                                                                                         cagcgcgtgtcccccGGGGGCAATCTTGAGCCCCAACCCA‑TCAGGGGTTCTGGAAAACATAACCATTGAAAAAAACAGGGTAAATTTCGCTCCCT‑TCCGGGCCGCTCCACGCTTGAAAGGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGGTGTTAGCGGGTGTTGGGGTGACGGGTGTTCCGGGAGTGGAATAAACG‑GGAGTAACAGCACCGCCATTATTATCCGGCGTAACGGAACAGGAg                                 <  1:96659/240‑1 (MQ=255)                                                                                                                                                                                                                         ggagcaatcttaacccccacccca‑TCAGGCGTTTTTGTAAACACTACCAGTAAAAATAACATGGTAAATTTCGCTACCT‑CCCGGGCCGCGCCCCGCTTGCAATAAGTAATCGTTAACGCTCTCCGCGGGGGTGGAGGGCTGGTTGTTACCCGGAGTGGGGGTCCCCGGTGTTGCGGGAGTTGAACCAACC‑GGAGCCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACt                 <  2:11923/236‑1 (MQ=255)                                                                                                                                                                                                                          gggCAATCTTACGCCCCACCCCA‑CCAGGGTTTTTTGTAAACATTACCATAGAAAATAACATGGAAAATTTCGCTCCCT‑TCCTGGCCGTTCCACGCTTGCAATGAGTAATGTTTAACGCTACCCGCTGTTTTTGATGCCTCGCTGTTATCCGGATTGGGGGTCCCGGGTTTTGCGGGAGTGGGACCACCG‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACt                 <  1:91973/250‑1 (MQ=255)                                                                                                                                                                                                                                tCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGCAACTAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGCGCCAGGGTTACCATACTGGCTGAt          >  2:371862/1‑251 (MQ=255)                                                                                                                                                                                                                                tCTTTAGACCCAACCCA‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGAt          >  2:209879/1‑251 (MQ=255)                                                                                                                                                                                                                                        cccaaccca‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttg  >  1:29211/1‑251 (MQ=255)                                                                                                                                                                                                                                          cacccca‑TCAGGCTTTCTTGAAAACATTACCATAAAAAATAACATGGTAATTTTCGCTCCCT‑TCCGGGCCGCTCCCCGCTTGCAATGATTAATCGTTAACGCTATCCGCGGGTGTTGAGGGCTCGCTGTTATCCGGAGTGGGGGTCCCCGGTGTTTCGGGAGTGGGATCACCG‑GGATTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGAtt         <  1:372134/239‑1 (MQ=255)                                                                                                                                                                                                                                          caaccca‑TCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGAtt         >  2:372134/1‑242 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    TGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGTACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCCA‑CCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCT‑TCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGCTCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACC‑GGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTG  >  NC_000913/3468528‑3469009

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