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	98.1%	199.1 / NA	54	W362* (TGG→TGA)	chiA	endochitinase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (24/28);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTT  >  NC_000913/3468522‑3469014                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        gCCCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAgg                                                                                                                                                                                                                                                     <  2:32081/251‑1 (MQ=255)        gCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCtt                                                                                                                                                                                                                                              >  1:89691/1‑251 (MQ=255)                   gACAGAAGTTTTGTAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTAc                                                                                                                                                                                                                                   >  1:392549/1‑251 (MQ=255)                     aaGAAGTTTTATAATGGGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCa                                                                                                                                                                                                                                 <  1:19580/250‑1 (MQ=255)                      aGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCAt                                                                                                                                                                                                                                >  1:67146/1‑251 (MQ=255)                      aGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCAt                                                                                                                                                                                                                                <  1:431794/251‑1 (MQ=255)                          gTTTTATAAGGGGCACCGTAATAATCTACGAGGCTGTCCGGGC‑GTCACGCTTTGCGGGACTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGan                                                                                                                                                                                                                            <  2:248188/251‑2 (MQ=255)                                taATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAAc                                                                                                                                                                                                                      >  1:147250/1‑251 (MQ=255)                                          cGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt                                                                                                                                                                                                             <  1:419268/250‑1 (MQ=255)                                          cGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAAGGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAt                                                                                                                                                                                                             <  2:169826/250‑1 (MQ=255)                                           gTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAAttt                                                                                                                                                                                                           >  2:461155/1‑251 (MQ=255)                                             aataATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCg                                                                                                                                                                                                         <  2:192656/250‑1 (MQ=255)                                              ataATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCg                                                                                                                                                                                                         >  1:192656/1‑250 (MQ=255)                                                aaTCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAACCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTa                                                                                                                                                                                                      <  1:276084/251‑1 (MQ=255)                                                       gATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCt                                                                                                                                                                                               <  1:543985/251‑1 (MQ=255)                                                          gCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGACAATAACATGGTAAATTTCGCTACCTTCCTGGc                                                                                                                                                                                            >  1:348623/1‑251 (MQ=255)                                                                      gT‑ACGCTTTGCTGGCCTGGAACCCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCtt                                                                                                                                                                                <  2:147250/251‑1 (MQ=255)                                                                         aCGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGc                                                                                                                                                                              >  2:217554/1‑251 (MQ=255)                                                                         aCGCTTTGCGGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTg                                                                                                                                                                               <  2:392549/250‑1 (MQ=255)                                                                                     ccTGGAACCCAACCTCAACACTGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCg                                                                                                                                                                  <  2:305279/251‑1 (MQ=255)                                                                                        ggAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGtt                                                                                                                                                                >  1:512304/1‑250 (MQ=255)                                                                                        ggAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGt                                                                                                                                                                 <  1:348489/249‑1 (MQ=255)                                                                                                 ccGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATcc                                                                                                                                                      <  1:217554/251‑1 (MQ=255)                                                                                                 ccGCAACAAGGACGCCGCCGTTATCAATTTCACAGGGAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATcc                                                                                                                                                      <  2:125870/251‑1 (MQ=255)                                                                                                       caaggACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGgt                                                                                                                                                <  2:40629/247‑1 (MQ=255)                                                                                                                     ttctCAATTTCACAGGTACTCGAGTTACCGTACTGACTCCGTTCCGCAGCAGTAGCTGTACGCTCGAGACGCCATGGGTAATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCAAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAATTTTCGCTACATTCCTGGCCGCTCCACGCTTGCCATGAGTAATCCTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCt                                                                                                                                 >  2:485467/4‑251 (MQ=255)                                                                                                                     ttATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCt                                                                                                                                 >  2:144612/1‑251 (MQ=255)                                                                                                                       aTCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCAGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGt                                                                                                                               >  1:7171/1‑251 (MQ=255)                                                                                                                         cAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGtt                                                                                                                              >  1:22786/1‑250 (MQ=255)                                                                                                                          aaTTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGCTGGC‑TCGCTGTTAt                                                                                                                            >  2:297660/1‑251 (MQ=255)                                                                                                                              tCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCg                                                                                                                         <  1:376430/250‑1 (MQ=255)                                                                                                                                  aGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGtt                                                                                                                    <  1:209844/251‑1 (MQ=255)                                                                                                                                      aGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTgggg                                                                                                                <  1:305347/251‑1 (MQ=255)                                                                                                                                         cGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCa                                                                                                             <  2:512304/251‑1 (MQ=255)                                                                                                                                               aCCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGgtg                                                                                                       >  1:496274/1‑251 (MQ=255)                                                                                                                                                ccGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGgtgt                                                                                                      >  1:280931/1‑251 (MQ=255)                                                                                                                                                        aCTCAATTCCGCAGCGGTAGCTGTACGCTCGTGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCAGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCgg                                                                                                >  2:299609/1‑249 (MQ=255)                                                                                                                                                              ttCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGAt                                                                                        >  2:341520/1‑251 (MQ=255)                                                                                                                                                                ccGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATc                                                                                       >  2:164050/1‑250 (MQ=255)                                                                                                                                                                  gCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAAc                                                                                    <  2:67146/251‑1 (MQ=255)                                                                                                                                                                    aGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTCCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCg                                                                                  >  2:117838/1‑251 (MQ=255)                                                                                                                                                                       ggTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCgg                                                                                 >  1:430314/1‑249 (MQ=255)                                                                                                                                                                          ngCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTc                                                                             >  2:248771/2‑250 (MQ=255)                                                                                                                                                                            cTGTACGCTCGAGACGACATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATACGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGtaac                                                                           >  2:453672/1‑247 (MQ=255)                                                                                                                                                                                aCGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACAGTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGctc                                                                      <  2:356600/251‑3 (MQ=255)                                                                                                                                                                                  gCTCGAGACGCCATGGGTTAGTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGTCCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCg                                                                    <  2:69651/251‑1 (MQ=255)                                                                                                                                                                                     cGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGc                                                                   <  1:461155/249‑1 (MQ=255)                                                                                                                                                                                              aTGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ATCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATcc                                                        >  2:556658/1‑251 (MQ=255)                                                                                                                                                                                                                 agcggtgttACCCCGGGGGAATTCTTTAGACCCAACCC‑ATCGTGCTTTTTTTTAAATATCCCAATTAAAAATAACTTGTAAAATTTCGCTCCCTCCCTGGCCTCTCCACGTTTGAAATTATTACTCGTTAGCGCTTTCCGTTGTTTTTGATGTC‑TCGTTTTTATCCGGGGTGGGGGTCCCCTTTTTTTCCGGGTTTGGATCAACGGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGc                                     <  2:458838/243‑1 (MQ=255)                                                                                                                                                                                                                   cgctgggcccccGTGGGCATTTTTTAG‑CCCAACCCCACCAGGCGTTTTTAAAACCATTACCATTAAAAATAAAAGGTTAAATTCCGCTCCCTCCCTGGCCGCTCCCCGTTTGAAATGAGTAATCGTAAACGCTACCCGCTGGTGTTGATGGC‑TGGCTGTTCTCGGGGGTGGGGGTCCCCGGTGTTTCCGGAGTTGGATCAACGGGGGCCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCa                                   <  2:347078/243‑1 (MQ=255)                                                                                                                                                                                                                           cccggTGGGCAAGCTTTAGCCCCACCCC‑AGCGGGCGTTTTTGAAAACATCACCATTAAAAATAACATGTTAATTTTCGCCACCTTCCGGGCCGCTCCACGCTTGAAATGAGTAATCGTTAACGTTCCCCGCTGGTGTGGAGGGC‑TCGCTGTTACCGGGGGTGGGGGCCCCCGGTGTTTCCGGATTTGGAGCAACGGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTAcc                           <  2:89691/247‑1 (MQ=255)                                                                                                                                                                                                                             ccttGGGCATTCTTTAGCCCCACCCC‑ATCAGGCGTTCTTGTAAACATTACCATTAAAAATAACATGTTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTAGCCGCGGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGTCCACGGGTGTTTCCGGAGTTGGATCACCCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCAt                         <  1:545428/248‑1 (MQ=255)                                                                                                                                                                                                                                  ggCAATTTTTAGCCCAAACCC‑ATCAGGCGTTCTTGTAAACATTACCATAAAAAATAACATGGAAAATTTCGCTCCCTTCCGGGCCGCCCCACGCTTGAAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTTTTACCGGGAGTGGGGGCCACGGGTGTTTCCGGAGTGGGATCAACGGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTgg                    <  1:248771/251‑1 (MQ=255)                                                                                                                                                                                                                                           tagccccacccc‑acCAGGCTTTTTTGTAAACATTACAATTAAAAAAAAAGTGGTAATTTCCCCTCCTTTCCGGGCCGCTCCACTCTTGCAATGATTAATCTTTAACGCTACCCGCTGGTGTTG‑TGGCCTCGCTTTTATCCGGTGTTGGGGCCCCCGTTGTTTCCGGATTTGGACCACCCGGAGTCCCGGCTCCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTc            <  2:496274/246‑1 (MQ=255)                                                                                                                                                                                                                                                    ccc‑aTCAGGCGTTCTTGTAAACATTACCATTAAAAATAACATGGTAATTTTCGCTACCTTCCTGGCCGCTCCCCGCTTGCAAGGAGTAATCTTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTTTTATCCGGAGTTGGGGTCCCCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTgttgccgtt  <  1:556658/251‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        GCCCCATGCGTCTTGATCGACAGAAGTTTTATAATGTGCATCGTTATAATCTACGATGCTGTCCGCGCTGT‑ACGCTTTGCTGGCCTGGAAACCATCCGCAACAATGACGCCGCCGTTATCAATTTCACAGGTAGTCGGGTTACCGTACTGACTCAATTCCGCAGCGGTAGCTGTACGCTCGAGACGCCATGGGTTATTGGAGTTTTCAGCGCTGGTACCCCGTGGGCAATCTTTAGACCCAACCC‑ACCAGGCGTTCTTGTAAACATTACCATTGAAAATAACATGGTAAATTTCGCTACCTTCCTGGCCGCTCCACGCTTGCAATGAGTAATCGTTAACGCTATCCGCTGGTGTTGATGGC‑TCGCTGTTATCCGGAGTTGGGGTCACCGGTGTTTCCGGAGTTGGATCAACCGGAGTCACAGCACCGCCATTATTATCCGGCTTAACGGAACAGGAGCCAGGGTTACCATACTGGCTGATTTCTGTTGCCGTT  >  NC_000913/3468522‑3469014

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: ‑