Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,479,717	G→C	S176S (TCG→TCC)	unknown →	Flp pilus assembly protein RcpC/CpaB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,479,717	0	G	C	100.0%	134.3 / NA	49	S176S (TCG→TCC)	unknown	Flp pilus assembly protein RcpC/CpaB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (29/20);  total (29/20)

ACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAACGATCCTAGATTCTTTACATCCAGCCTAGATCCGAATCAAGAAGTTGCATATCGCATCTATGTTAAATCTGGTGATGATTACATTCTAAATACTATTTTAAGTGGAGACTATGTTGCTATCTATAGCCAACAAATTGCTTCAGATTCGCGCAACACTTATGATCGTAAAGATTTAGTTAAAATTTTAGGTAAAGTATTGGTTTTACAAGCAAGAAACTTTACTCAATCAACAGAACAAGGTGCTCCTGAAGCTGAAAAGGCAGCA‑GA‑CGGGATTCCAAACAATCAAGAATATGCTGGTTATGTAGCACTAAAAATAACTGCTGAACAAGCAAAAAAATTCTATTCACTGGATAGAGAATCAAAGTTAGTTATATTACCTGCCGAAAAT  >  1/1479475‑1479957                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   aCTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCaatac                                                                                                                                                                                                                                            <  1:769804/251‑4 (MQ=255)      aaaTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATg                                                                                                                                                                                                                                       <  1:98335/251‑1 (MQ=255)        aTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCaatac                                                                                                                                                                                                                                            >  1:187440/1‑241 (MQ=255)        aTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGAt                                                                                                                                                                                                                                     <  1:10668/251‑1 (MQ=255)        aTGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGAt                                                                                                                                                                                                                                     <  1:626647/251‑1 (MQ=255)           gTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGt                                                                                                                                                                                                                                  >  1:218883/1‑251 (MQ=255)                  gggATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATg                                                                                                                                                                                                                                       <  1:209340/239‑1 (MQ=255)                   ggATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGa                                                                                                                                                                                                                             <  1:667185/248‑1 (MQ=255)                     aTTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAg                                                                                                                                                                                                                        >  1:553043/1‑251 (MQ=255)                           aGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGtt                                                                                                                                                                                                                      <  1:8770/247‑1 (MQ=255)                             gCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTaaagtt                                                                                                                                                                                                                >  1:283197/1‑248 (MQ=255)                              cTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGttt                                                                                                                                                                                                               >  1:263498/1‑251 (MQ=255)                                aTTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTa                                                                                                                                                                                                             >  1:105407/1‑251 (MQ=255)                                 tttAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAg                                                                                                                                                                                                            <  1:35820/251‑1 (MQ=255)                                         cAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggc                                                                                                                                                                                                     <  1:366099/250‑4 (MQ=255)                                               aTGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggcagtggtt                                                                                                                                                                                              >  1:605186/1‑241 (MQ=255)                                                          ggTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacg                                                                                                                                                                                   >  1:412097/1‑249 (MQ=255)                                                          ggTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaag                                                                                                                                                                                       >  1:141345/1‑247 (MQ=255)                                                           gTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacgg                                                                                                                                                                                  >  1:682175/1‑248 (MQ=255)                                                              cTTTTCTATCCAACAAATCTATTATCTTTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGAt                                                                                                                                                                                                                                     <  1:432130/197‑1 (MQ=255)                                                                tttCTATCCAACAAATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattt                                                                                                                                                                             >  1:564018/1‑243 (MQ=255)                                                                        caacaaATCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCTTATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggc                                                                                                                                                                                                     <  1:841199/219‑4 (MQ=255)                                                                             aaTCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTaaggcagtggtt                                                                                                                                                                                              >  1:716179/1‑211 (MQ=255)                                                                              aTCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCCAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattttgctcaagtagcag                                                                                                                                                               >  1:87582/1‑229 (MQ=255)                                                                               tCTATTATCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattttgctcaagtagcag                                                                                                                                                               >  1:154701/1‑228 (MQ=255)                                                                                 tattatCTCTCCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTAcaagcacggaattttgctcaagtagcagaac                                                                                                                                                            <  1:47873/251‑26 (MQ=255)                                                                                      tctctcCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggc                                                                                                                                                       <  1:740715/251‑2 (MQ=255)                                                                                         ctcCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggca                                                                                                                                                      >  1:91141/1‑247 (MQ=255)                                                                                         ctcCGAATGATCCTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCCATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggca                                                                                                                                                      >  1:54289/1‑247 (MQ=255)                                                                                                    ccTAGATTCTTTACTTCCAGCTTGAATCCAAGTCAAGAAGTCGCATATCGCATCTATGTTAGATCTGGCGATGATTACATTCTAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAggcactgaaactgagaa                                                                                                                                         >  1:167761/1‑236 (MQ=255)                                                                                                                                                                                      tAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTc                                                  <  1:543537/251‑1 (MQ=255)                                                                                                                                                                                      tAAATACTATTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTc                                                  <  1:2960/251‑1 (MQ=255)                                                                                                                                                                                              aTTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTA‑GATCTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCa                                            >  1:454188/1‑249 (MQ=255)                                                                                                                                                                                              aTTTTAAGCGGAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTA‑GATCTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCa                                            >  1:448427/1‑249 (MQ=255)                                                                                                                                                                                                      cggAGACTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATagag                                  <  1:885580/250‑1 (MQ=255)                                                                                                                                                                                                        gagaCTATGTATCTGTTTATAGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTGTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAgaga                                 <  1:453155/250‑1 (MQ=255)                                                                                                                                                                                                          gactatgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTc                                                  >  1:279359/16‑231 (MQ=255)                                                                                                                                                                                                            ctatgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGAt                                      <  1:459579/228‑1 (MQ=255)                                                                                                                                                                                                            ctatgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATcga                            >  1:62278/14‑249 (MQ=255)                                                                                                                                                                                                            ctatgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATcga                            >  1:773401/14‑249 (MQ=255)                                                                                                                                                                                                            ctatgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATcga                            >  1:479787/14‑249 (MQ=255)                                                                                                                                                                                                            ctatgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATcga                            >  1:894938/14‑249 (MQ=255)                                                                                                                                                                                                              atgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATCGAag                          >  1:780671/12‑251 (MQ=255)                                                                                                                                                                                                              atgtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATCGAag                          <  1:551970/240‑1 (MQ=255)                                                                                                                                                                                                                gtatctgtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATCGAagtt                        >  1:606079/10‑251 (MQ=255)                                                                                                                                                                                                                      gtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATCGAagttagtcat                  >  1:877339/4‑248 (MQ=255)                                                                                                                                                                                                                      gtttataGCCAGCAAACCGCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATCGAagttagtcat                  >  1:464845/4‑248 (MQ=255)                                                                                                                                                                                                                                agcaaaccgCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAgaga                                 >  1:445055/9‑226 (MQ=255)                                                                                                                                                                                                                                      ccgCTTCAGATTCCCGCAATACTTATGATCGTAAAGATTTAGTTAAAGTTTTAGGTAAGGCAGTGGTTCTACAAGCACGGAATTTTGCTCAAGTAGCAGAACAAG‑‑GC‑ACTGAAACTGAGAAAGCGTTAGGT‑CTGG‑CTCCAAATAGTGAAGTATATGCCGGTTATGTAGCATTAAAAATAACTGCTGAACAAGCTAAAAAATTCTATTCACTGGATAGAGAATCGAAGTTAGTCATATTACCTGTCGAAaat  <  1:725284/249‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                   ACTAAAAATGGTATTCAGGGATTACAAGGCTATTTAATGGCAGAAAATGTCAGAGCTGGTTCTTTTCTATCCAACAAATCTATTATCTCTCCGAACGATCCTAGATTCTTTACATCCAGCCTAGATCCGAATCAAGAAGTTGCATATCGCATCTATGTTAAATCTGGTGATGATTACATTCTAAATACTATTTTAAGTGGAGACTATGTTGCTATCTATAGCCAACAAATTGCTTCAGATTCGCGCAACACTTATGATCGTAAAGATTTAGTTAAAATTTTAGGTAAAGTATTGGTTTTACAAGCAAGAAACTTTACTCAATCAACAGAACAAGGTGCTCCTGAAGCTGAAAAGGCAGCA‑GA‑CGGGATTCCAAACAATCAAGAATATGCTGGTTATGTAGCACTAAAAATAACTGCTGAACAAGCAAAAAAATTCTATTCACTGGATAGAGAATCAAAGTTAGTTATATTACCTGCCGAAAAT  >  1/1479475‑1479957

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