Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_007779	4,575,860	A→C	23.7%	V132G (GTG→GGG)	Y75_RS22690 ←	PLP‑dependent aminotransferase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779	4,575,860	0	A	C	23.7%	129.0 / 2.4	60	V132G (GTG→GGG)	Y75_RS22690	PLP‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (27/16);  new base C (2/15);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 4.37e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCT  >  NC_007779/4575627‑4576094                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                gcagcaTCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGcaccac                                                                                                                                                                                                                                             <  1:484729/238‑1 (MQ=255) gcagcaTCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCaccacc                                                                                                                                                                                                                                            <  1:2221263/239‑1 (MQ=255)      acctccaccGAACCTGAATAACAGGGGGATTCGCCGCCGACAATATCCCCGGTTTTACCCACTGCCATTACCGCCAGCGACATCGAGTTATGGCGGCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCt                                                                                                                                                                                                                                       <  2:860184/231‑1 (MQ=255)                               gggATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGG‑GACCCCCCAGCCCCCGTCCGGCATCCTGGGGGGACCCCGCTCTCGCCATACCCGGCGATCCGGCAAAC‑AGGCCTAACCGACAAAGGG‑TTTGTGCATT‑TTCTTCCCCACCCCGCTCTGCTCCCCCCACCGCGGGTTccg                                                                                                                                                                                                               >  2:1704671/1‑236 (MQ=255)                                                                     gCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGC‑GAACCCCCCGCCCCCGTGGGGGAACCGGGGGGGACACTGCCCGCGCAATAACCGCCGC‑CCCGCCCAAAACCCATAAACCA‑AAACGGCTTTGC‑GATTATGGCTCACCCCCCCGCGCCGCCCCCCCCCCAGCGGGTTC‑ACGCCGGGCC‑CGTCGACATCCGGGGCGCTTTTGCTTAAccg                                                                                                                                                                      >  2:118322/1‑237 (MQ=255)                                                                                   cGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCCCCCAGCCGCCGTCGCGGCTCAAGCGGGCGATCCGCTCTCTCCAAACTCCCCGA‑CCGGCCAACAAATCACAAACCA‑GAACAG‑GTTGCGGAATTTTTTTCACCCCCCCGCCTCGCCCCCCCCCCACCCGGTTC‑CGGGTTGGCGCCTTCCCCTCCCGGTGGTCTTTT‑CTTTACGGGAAAATGGCACTGt                                                                                                                                                        >  2:2212652/1‑239 (MQ=255)                                                                                          gttttggggcggcCGCTGGTGTGGGTGTTGTCTTGGGCGGT‑GACCCCCGAGCGCCTGTCGAGCATCAGGGGGGCAATCTGCTCCCGCAATCCTCGCGGC‑CCGGCTAACAAGTCATACCCGA‑GACCGG‑TTGGCAGATTATGCGGCCCCCCCCGGTTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGc                                                                                                                                                     <  2:928381/228‑1 (MQ=255)                                                                                                tGGCAGCCGCTGGTGATGATGATGTCATCGGCGGTGGACCA‑CGAGCCCCCGTCGGGCCACAGAGGGGGAAACGGCCCTTTCCATAAACGCCGC‑CCGGCGAACAAATCCTAACCGG‑AGACGG‑TTTGCCAATTATGGTGCACCACCCCGCATCGCGCCCGCCACCGCGGTTTC‑AGGCTGGGGG‑CGTCGTCACCCGGTGTGCTTTTGCTTTACAGAAAACAGGAACGGTCGcccgtcgcca                                                                                                                                            >  2:930954/1‑229 (MQ=255)                                                                                                     gCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCACCCAGCCGCCGGCGGGGATCAAGCGGGGAAACTGCTCTCGCCATACTCGCTGA‑CCCGCCAAAAAAACATAAACCA‑GAAAGG‑GGTGCAGATTATTCTGCCCCCCCCCGCCCTTCTCCCCCCCCCCCGGGTTCTGGGCTGGG‑G‑CGTCCCCATCCCGGGGGGGTTTGGTTTACGGGACAAAGGAAACGTCGCTATTCGCCCCCCGc                                                                                                                                       >  2:1706038/1‑238 (MQ=255)                                                                                                      ccGCTGGTGAGGATGATGTCATCGGCGGT‑GCCCCCCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCa                                                                                                                                      <  1:1253800/238‑1 (MQ=255)                                                                                                             ggatgatgatGTCATGGGCGGT‑GCCCCCCGAGCCGCTGCCGAGCACCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAg                                                                                                                               <  1:225000/237‑1 (MQ=255)                                                                                                                                            aGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGATAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCgg                                                                                                 >  1:1797193/1‑238 (MQ=255)                                                                                                                                              gcGCTGTCGAGCAACAGGCGGGAAATCTGCTCCCCCAAACCCCGCTGC‑CGGGCTAACAAGTCATACCCGA‑TAACGG‑TTGGCAGATTAGCCTGCCCCGCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGgcg                                                                                              <  2:311362/238‑1 (MQ=255)                                                                                                                                               cGCTGTCGGGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGgcg                                                                                              <  1:1490353/238‑1 (MQ=255)                                                                                                                                                 cTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCt                                                                                            >  2:1457571/1‑238 (MQ=255)                                                                                                                                                  tGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTg                                                                                           >  2:1789724/1‑238 (MQ=255)                                                                                                                                                         cATCAGGGGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGCCAAAACCAA‑GAACGG‑TTTGCAGATTATGCTGCCCCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGa                                                                                    <  1:770495/238‑1 (MQ=255)                                                                                                                                                          aTCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GACCGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGa                                                                                    <  1:469739/237‑1 (MQ=255)                                                                                                                                                              ggCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTc                                                                              >  2:643316/1‑239 (MQ=255)                                                                                                                                                                 gggAAATCTGCCCGCGCAATACTCGCTGA‑CCGGCTACCAAGTCATAACCGA‑GACCGG‑TTTGCAGATTATGCTGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCAt                                                                            <  1:1576214/238‑1 (MQ=255)                                                                                                                                                                            tctcGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑Tg‑g                                                              >  2:1189989/1‑239 (MQ=255)                                                                                                                                                                             ctcGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcggcgggaactgc                                                               >  2:2176013/1‑232 (MQ=255)                                                                                                                                                                               cGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGtgactgcgg                                                             >  1:1582238/1‑231 (MQ=255)                                                                                                                                                                                aaaaaccccGCTGC‑CAGGCAAACAGGTCAAAACCGA‑GAGCGG‑TTGGCAGTATTTGCGGCACCCCCCGGTTTCTCCCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑Tg‑gcggc                                                          <  2:1022068/231‑1 (MQ=255)                                                                                                                                                                                 aaaTACCCGCTGC‑CCGCCTAACAGGCCATACCCGA‑GAAGGT‑TTTGCAGATTATGCTGCCCCCCCCGCTTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑t                                                        <  1:400668/238‑1 (MQ=255)                                                                                                                                                                                  aaTACTCGCTGC‑CGGGCTAACAAGTCATACCCGA‑GAACGG‑TTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑t                                                        <  1:165132/238‑1 (MQ=255)                                                                                                                                                                                  aaTACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCGGCACCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑t                                                        <  1:1704671/238‑1 (MQ=255)                                                                                                                                                                                     aCTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCg                                                                                                                                                                                         >  2:403016/1‑109 (MQ=255)                                                                                                                                                                                     aCTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCg                                                                                                                                                                                         <  1:403016/109‑1 (MQ=255)                                                                                                                                                                                     aCTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGCAA‑Ta‑gcggcctggg                                                     >  1:597711/1‑234 (MQ=255)                                                                                                                                                                                        cGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCa                                                                                                                             >  1:104139/1‑166 (MQ=255)                                                                                                                                                                                        cGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCttt                                                 >  2:1999336/1‑239 (MQ=255)                                                                                                                                                                                        cGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGCAC‑Tg‑gcgggcggggcgt                                                  >  1:462748/1‑230 (MQ=255)                                                                                                                                                                                         gCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGAC‑TG‑GCGGC‑Cgggggttt                                                >  2:845081/1‑234 (MQ=255)                                                                                                                                                                                         gCTGA‑CCGGCTAACAAGTCATAACCGA‑GAAAGG‑TTTGCAGATTATGATGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACACGTCATCGGCGGCAC‑Tc‑gccggccggggtg                                                  >  1:1491501/1‑227 (MQ=255)                                                                                                                                                                                          cTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGaa                                                                                                                                                            <  1:1644716/133‑1 (MQ=255)                                                                                                                                                                                          cTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGaa                                                                                                                                                            >  2:1644716/1‑133 (MQ=255)                                                                                                                                                                                          cTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGCTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTCA‑Cg‑gcggcgtgggtgtac                                                >  1:1568061/1‑229 (MQ=255)                                                                                                                                                                                            tgcccGGCTAACAAGTCATACCCGA‑GAACGG‑TTTGCATATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACg                                              <  1:1270959/236‑1 (MQ=255)                                                                                                                                                                                               ccGGCTACAACGTCAGACCGGA‑AAACGG‑TTTGCAGTTTGTGCTGCCCCCCCGGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTtg                                           <  1:1520994/239‑1 (MQ=255)                                                                                                                                                                                               ccGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGGAC‑TG‑GCGGC‑gggggggttcgatg                                           >  2:1138514/1‑229 (MQ=255)                                                                                                                                                                                                    tAACAAGTCATACCCGA‑GACCGT‑TTTGCAGATTATGCTGCCCCCCCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGa                                      <  1:1706038/239‑1 (MQ=255)                                                                                                                                                                                                          gTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTaa                                                                                                                                                                         >  1:1053558/1‑105 (MQ=255)                                                                                                                                                                                                          gTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTaa                                                                                                                                                                         <  2:1053558/105‑1 (MQ=255)                                                                                                                                                                                                           tCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCCCCCCCCGGCTTAGCCCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCCCCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑Tg‑g                                                              <  1:1435181/209‑1 (MQ=255)                                                                                                                                                                                                             aTAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTg                                                                                                                                                                              <  1:352342/97‑1 (MQ=255)                                                                                                                                                                                                             aTAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTg                                                                                                                                                                              >  2:352342/1‑97 (MQ=255)                                                                                                                                                                                                             aTAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCa                                                                                                                                >  1:81301/1‑143 (MQ=255)                                                                                                                                                                                                                aCCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTACTTG‑CCGCC‑CGGGCTTTTCTTTTGGCGCCAAATAAAcc                            >  1:1890199/1‑238 (MQ=255)                                                                                                                                                                                                                       aaCGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGCCATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTCC‑TGCGGGGCGGGGGCTCGACCTTGTGTGCCAAAAAAAccaaagcg                      >  2:750032/1‑233 (MQ=255)                                                                                                                                                                                                                           g‑TTTGCGGATATTGCGGCCCCACCGGGCTTGGTTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCt                 <  1:1138514/238‑1 (MQ=255)                                                                                                                                                                                                                           g‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcgg                                                                         <  1:368020/185‑1 (MQ=255)                                                                                                                                                                                                                           g‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATcgg                                                                         >  2:368020/1‑185 (MQ=255)                                                                                                                                                                                                                              tggCATATTACGCTGCACCCCCCGGTTTGGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCg              <  2:1890199/237‑1 (MQ=255)                                                                                                                                                                                                                                gCAGATTATGCGGCCCCCCCGGGCTTAGCCCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGa            <  1:1838700/239‑1 (MQ=255)                                                                                                                                                                                                                                 cAGATTCTGCTGCCCCCCCCGGTTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGa            <  1:845081/238‑1 (MQ=255)                                                                                                                                                                                                                                     ttATGCTGCCCCCCCCGTCTTGGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtgat       <  2:1568061/239‑1 (MQ=255)                                                                                                                                                                                                                                     ttATGCTGCACCCCGCGGCTTAGCCTCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtga        <  1:919702/238‑1 (MQ=255)                                                                                                                                                                                                                                     ttATCCTGCCCCGCCGGGTTTGGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtg         <  1:29115/237‑1 (MQ=255)                                                                                                                                                                                                                                      aaggcgccaccacCCGGCTGAGCCCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CGTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGtgatg      <  1:866172/232‑1 (MQ=255)                                                                                                                                                                                                                                         gctgcACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGt                                                                                                             >  2:2127601/1‑136 (MQ=255)                                                                                                                                                                                                                                         gctgcACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGt                                                                                                             <  1:2127601/136‑1 (MQ=255)                                                                                                                                                                                                                                           tgcACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCt  <  1:1189989/238‑1 (MQ=255)                                                                                                                                                                                                                                           cccaccccCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCt  <  1:452891/236‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                GCAGCATCTGCATCGAACCGTAATAACAGGGGGATTCGACCGCGACAATATCGCCCGGTTTACACACTGCCATTAACGCCAGCGACATCGAGTTATGGCAGCCGCTGGTGATGATGATGTCATCGGCGGT‑GACCACCGAGCCGCTGTCGAGCATCAGGCGGGCAATCTGCTCTCGCAATACTCGCTGA‑CCGGCTAACAAGTCATAACCGA‑GAACGG‑TTTGCAGATTATGCTGCACCACCCGGCTTAGCTCCCGCCACAGCGGTTTC‑AGGCTGGGCG‑CTTCGACATCCGGTGTGCTTTTGCTTAACGGAACAATGGAACTGTCGCTATGCGCCTCCAGCATATCCAGCACCTGATCCCACTGGGTAATTTCCACCGGGCGCTGCACCGGACGCGTCATCGGCGGTAC‑TG‑GCGGC‑TGGGCTTTACGTTGTGCGACAAAATAACCCGAACGCGGCTGCGGAGTGATGAGCT  >  NC_007779/4575627‑4576094

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