Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	78,543	C→T	A1914T (GCC→ACC)	unknown ←	Autotransporter adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	78,543	0	C	T	91.9%	176.4 / ‑4.2	74	A1914T (GCC→ACC)	unknown	Autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (50/18);  minor base G (0/6);  total (50/24)
Fisher's exact test for biased strand distribution p-value = 7.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAAAAGTGCGGTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGCGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCGTGGCTTGTCATGCTTGGACCGT  >  1/78312‑78788                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      tAAAAAAGTGCGGTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    <  1:842451/251‑1 (MQ=255) aaaaaaaGGGGGGTGAAGAAACACATTATTTTTCCCCCGCATTTTTGTTTCAGGTTAGCTGAATAAGCGAATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    <  1:779922/250‑1 (MQ=255) aaaaaaaGGGCGGTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTAAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCCCCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    <  1:870257/250‑1 (MQ=255)     aaaGTGCGGTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGcc                                                                                                                                                                                                                                               >  1:574751/1‑236 (MQ=255)             gTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAGGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    >  1:645804/1‑239 (MQ=255)             gTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCaaa                                                                                                                                                                                                                        >  1:456763/1‑251 (MQ=255)              tGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAAc                                                                                                                                                                                                                       <  1:832906/251‑1 (MQ=255)                  gaaACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    >  1:641598/1‑234 (MQ=255)                  gaaACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACccg                                                                                                                                                                                                                    >  1:771411/1‑250 (MQ=255)                     acacATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGccc                                                                                                                                                                                                                                              >  1:205275/1‑221 (MQ=255)                        caTTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGcc                                                                                                                                                                                                                                               <  1:154570/217‑1 (MQ=255)                          ttattTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGtctc                                                                                                                                                                                                                                         >  1:729471/1‑221 (MQ=255)                            attTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    <  1:829119/224‑1 (MQ=255)                            attTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAgg                                                                                                                                                                                                         >  1:819910/1‑251 (MQ=255)                                       gCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGtctc                                                                                                                                                                                                                                         >  1:97023/1‑208 (MQ=255)                                        cACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTg                                                                                                                                                                                             >  1:849966/1‑251 (MQ=255)                                          cTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGc                                                                                                                                                                                            >  1:790400/1‑250 (MQ=255)                                               gTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCtt                                                                                                                                                                                      >  1:706211/1‑251 (MQ=255)                                               gTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCtt                                                                                                                                                                                      >  1:18936/1‑251 (MQ=255)                                                     gtttagcaaaaaaagcttttttccccAGAAAACGCCCGCGCTTGCCGTACCGCCACCTTCCCTTT‑GTGGGTGGGTGCTAAAAGAACCTTT‑CG‑CGCCGTTTTTTCCGTTAGGGGAAAGGCGGGAGAACCCCACAGCTA‑CGGCGTTTTCGTTTTGATATGGCCCGACAGGGGCGACTAGGGTGGATTCCCCGTCTCTAGGGGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGa                                                                                                                                                                                 <  1:888135/226‑1 (MQ=255)                                                     ggTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    >  1:838224/1‑199 (MQ=255)                                                      tttagtaaaagaaGCGCTTTTTACAAGAATTAGCCCGCCCCGGCCGCCCCGCCAACTTCACTTG‑GTGGGTTGGGGC‑GAAAAAAAATTTCCG‑CGCCGTATTTTCCGTTATCGGAAAGGCGTGAGTACCCCACAGCTA‑CAGCGTTTTCGCTTTGATATTGACCGACAGCGGCGACTAGGGTGGATTGCCCGTCTTTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGAc                                                                                                                                                                                <  1:685216/239‑1 (MQ=255)                                                             gAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATc                                                                                                                                                                         >  1:243834/1‑250 (MQ=255)                                                             gAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATc                                                                                                                                                                         >  1:46862/1‑250 (MQ=255)                                                                tAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCAtt                                                                                                                                                                     >  1:34414/1‑251 (MQ=255)                                                                      cATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCAttttttt                                                                                                                                                                >  1:864342/1‑250 (MQ=255)                                                                            aCCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCAttt                                                                                                                                                                    >  1:526843/1‑240 (MQ=255)                                                                            aCCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCAcgcg                                                                                                                                                          >  1:479006/1‑250 (MQ=255)                                                                                      cccgcccccTGCCGTCCCGCCAACTCCACTTT‑GTGAGTGGGTGATAAAAAAAACTTC‑CC‑CGCCGTTTTTTCCTTTATCGGAAATGCGTGAGTACCCCACAGCTA‑CAGCGGTTTCGTTTTGATATTGACCGACAGGGGCGACTAGGGTGGATTGCCCGTCTCTATGCGGCTGCACCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCttatcaatgt                                                                                                                                               <  1:548816/244‑1 (MQ=255)                                                                                      ccccGCGCCTGCCGT‑CCCCCCACTTCCCCTTGGTGGGTGGGTGCTAAAAAAACCTTC‑CA‑CGCCGTTTTTTCCTTTATCGGAAAGGCGTGAGTACCCCACAGCTA‑CAGCGTTTTCGTTTTGATATTGCCCGACAGCGGCGACTAGGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCttatcaatgt                                                                                                                                               <  1:534534/250‑1 (MQ=255)                                                                                             ccTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCttatcaatgttatcaat                                                                                                                                        >  1:372094/1‑251 (MQ=255)                                                                                                           aaCTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAggg                                                                                                                                                                                                        >  1:165173/1‑173 (MQ=255)                                                                                                           aaCTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAAt                                                                                                                            >  1:337492/1‑249 (MQ=255)                                                                                                            atttccttgt‑gtgAGTTGGTGGTAAAAAAAATTTT‑CG‑CGCCTTTTTTTCCGTTATCGAAAATGCGTGAGTAACCCAGAGCTA‑CAGCGGTTTGCCTTTAATATTGCCCGACGGCGGCGACAAGGGTGGTTTCCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGc                                                                                                                         <  1:20735/242‑1 (MQ=255)                                                                                                            aCTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATc                                                                                                                                                                         >  1:590891/1‑203 (MQ=255)                                                                                                             tttcccctttttgggtgggTGCTAAAAAAATCTTT‑CC‑CGCCTTTTTTTCCGTTATCGAAAAGGCGTGAGTACCCCACAGCTA‑CAGCGTTTTCGTTTTGATATTGACCGACAGCGGCGACTAGGGTGGTTTGCCCGTCTCTATGCGGCTGCATAAACCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGc                                                                                                                         <  1:59288/234‑1 (MQ=255)                                                                                                                    tt‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGa                                                                                                                 >  1:267513/1‑251 (MQ=255)                                                                                                                    tt‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGa                                                                                                                 >  1:375356/1‑251 (MQ=255)                                                                                                                       gtgAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACg                                                                                                               >  1:651752/1‑251 (MQ=255)                                                                                                                        ggtgTTGGTGC‑GAAAAAAACCTT‑CCCCGCCTTATTTCCCTTTATCGAAAATGCGTGAGTACCCCAC‑GCCACCAGCGTTTTCGTTTTGATATTGACCGACGGCGGCGACAAAGGGGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACg                                                                                                               <  1:825543/247‑1 (MQ=255)                                                                                                                                            cTTT‑CA‑CGCCGTTTTTTCCGTTATCGAAAAGGCGTGAGTACCCCACAGCTA‑CAGCGTTTTCGTTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAACCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          <  1:365443/251‑1 (MQ=255)                                                                                                                                            cTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    >  1:81676/1‑113 (MQ=255)                                                                                                                                            cTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          >  1:800386/1‑251 (MQ=255)                                                                                                                                             tttccA‑CGCCGTTTTTTCCGTTATCGGAAATGCGTGAGAAACCCACAGCTA‑CAGGGGTTTCGCTTTGATATTGACCGACGGCGGCGACTAGGGGGGATGGCCCGTCTTTATGGGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          <  1:894000/247‑1 (MQ=255)                                                                                                                                              tttccccGCCGTTTTTTCCGTTATCGGAAATGCGTGAGTACCCCCCAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACGGCGGCGACTAAGGGGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          <  1:502976/245‑1 (MQ=255)                                                                                                                                              tttccccGCCGTTTTTTCCGTTATCGGAAAGGCGTGAGTACCCCACAGCTA‑CAGCGTTTTCGTTTTGATATTGCCCGCCAGCGGCGACTAAGGGGGATTGCCCGTCTCTATGCGGCTGCATCAACCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          <  1:542979/245‑1 (MQ=255)                                                                                                                                                  cccGCCGTATTTTCCTTTATCGGAAATGCGTGAGTACCCCACAGCTA‑CGGCGGTTTCGCTTTGATATTGCCCGACGGCGGCGACTAAGGTGGATTGCCCGTCTCTAGGCGGCTGCATCAACCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTg                                                                                      <  1:507807/249‑1 (MQ=255)                                                                                                                                                     cccGTATTTCCCGTTATGGGAAATGGGTGAGTACCCCACGGCTA‑CGGCGTTTTCGCTTTGATATTGCCCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTAGGGGGCTGCAACAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCAt                                                                                   <  1:395/250‑1 (MQ=255)                                                                                                                                                        gTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGa                                                                                                                                                                                                           >  1:614415/1‑128 (MQ=255)                                                                                                                                                        gTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCttatcaat                                                                                                                                                 >  1:54089/1‑186 (MQ=255)                                                                                                                                                        gTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGAc                                                                                                  >  1:283892/1‑233 (MQ=255)                                                                                                                                                        gTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCg                                                                                >  1:363833/1‑251 (MQ=255)                                                                                                                                                        gTATTTTCCGTTATCGGAAAGGCGTGAGTACCCCACAGCTA‑CAGCGTTTTCGTTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          <  1:464419/241‑1 (MQ=255)                                                                                                                                                          tttttCCTTTATCGAAAATGCGTGAGTACCCCCCAGCTA‑CAGCGTTTTCGTTTTGATATTGCCCGACGGCGGCGACTAGGGGGGATTGCCCGTCTCTATGCGGCTGCATCAACCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTc                                                                              <  1:393284/250‑1 (MQ=255)                                                                                                                                                          aTTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCttatcaatgt                                                                                                                                               >  1:262504/1‑186 (MQ=255)                                                                                                                                                          aTTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGt                                                                               >  1:681813/1‑250 (MQ=255)                                                                                                                                                           ttttCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCg                                                                             >  1:288101/1‑251 (MQ=255)                                                                                                                                                               ccGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTg                                                                          >  1:870733/1‑250 (MQ=255)                                                                                                                                                                     tCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTcc                                                                    >  1:496779/1‑250 (MQ=255)                                                                                                                                                                       ggAATTGCGTGAGTAACCCACGGCTA‑CGGCGGTGTCGCTTTGATATTGCCCGACGGCGGCGACTGAGGGGGAGTGCCCGTCTCTGGGCGGCTGCACCAACCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCt                                                                  <  1:4880/250‑1 (MQ=255)                                                                                                                                                                       ggAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCtt                                                                 >  1:549823/1‑251 (MQ=255)                                                                                                                                                                        gAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGc                                                                         >  1:67660/1‑242 (MQ=255)                                                                                                                                                                            tGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCa                                                            >  1:124524/1‑251 (MQ=255)                                                                                                                                                                                 gAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGc                                                                                                                                                                                                                                    >  1:292373/1‑78 (MQ=255)                                                                                                                                                                                   gTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGAt                                                                                                                                                                                                                                                   >  1:441297/1‑61 (MQ=255)                                                                                                                                                                                      aCCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCaa                                                  >  1:734783/1‑251 (MQ=255)                                                                                                                                                                                                               ttGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATAttt                                                                                          >  1:797503/1‑187 (MQ=255)                                                                                                                                                                                                               ttGATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATg                           >  1:685744/1‑250 (MQ=255)                                                                                                                                                                                                                 gATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCttatcaatg                                                                                                                                                >  1:598511/1‑131 (MQ=255)                                                                                                                                                                                                                 gATATTGACCGACAGCGGCGACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGGAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCg                        >  1:423091/1‑251 (MQ=255)                                                                                                                                                                                                                                cggcgACTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCGTGGCTTGTCa              <  1:357808/246‑1 (MQ=255)                                                                                                                                                                                                                                     agaaaGGTGGATGGCCGGTCTTTATGGGGCTGATTCAATCCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCGTGGCTTGTCATGCTTGGAc     <  1:354775/247‑1 (MQ=255)                                                                                                                                                                                                                                      cTAAGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCGTGGCTTGTCATGCTTGGACCg   >  1:209658/1‑251 (MQ=255)                                                                                                                                                                                                                                        aaGGTGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCGTGGCTTGTCATGCTTGGACCGt  >  1:255745/1‑250 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      TAAAAAAGTGCGGTGAAGAAACACATTATTTTTCCACCGCACTTTTGTTTCAGGTTAGCTGAATAAGCGCATTTTACCAGAAATAGCCCGCGCCTGCCGTACCGCCAACTTCACCTT‑GTGAGTTGGTGCTGAAAGAAACTTT‑CA‑CGCCGTATTTTCCGTTATCGGAAATGCGTGAGTAACCCACAGCTA‑CAGCGGTTTCGCTTTGATATTGACCGACAGCGGCGACTAAGGCGGATTGCCCGTCTCTATGCGGCTGCATCAAACCCGCCGTAGCGAGGGCGGAAGCGGTGCCGGCTTTGCGACGTTTATCCATTTTTTTCACGCGCTTATCAATGTTATCAATGCGATTATTAATGGCACTGCCGACGTCTTGAAGCTGACTCATATTTACTGCATCCGTCGGTGCTTTCCCTTTCGCAACGTTAGCAATTCGTTTGCCGCCGGCGTTGATGCCGTGGCTTGTCATGCTTGGACCGT  >  1/78312‑78788

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