Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	307,931	0	T	.	64.4%	32.7 / 25.6	50	coding (1402/2526 nt)	ecpC	putative fimbrial usher protein EcpC
Reads supporting (aligned to +/- strand):  ref base T (9/7);  new base . (19/10);  total (28/22)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CATCAAACTGTTTGCGTGCTGACAGGTTTGCCATGGTGTAGCCGTTTTGATGAGTCATCCCTGCGCTAAACCAGTTGCCCAGTGGTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTTCCCGTTTTCCGACCAGCGATCCATAT  >  NC_000913/307571‑308218                                                                                                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                         catcaAACTGTTTGCGTGCTGACAGGTTTGCCATGGTGTAGCCGTTTTGATGAGTCATCCCTGCGCTAAACCAGTTGCCCAGTGGTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑‑ATGGCTGTTGTAACGGCGGTCATC‑TTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑‑TGA‑TTGCCAATGCGGG‑‑‑‑‑‑‑TTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGa                                                                                                                                                                                                                                          >  1:32862/1‑418 (MQ=255)                               ccATGGTGTAGCCGTTTTGATGAGTCATCCCTGCGCTAAACCAGTTGCCCAGTGGTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCg                                                                                                                                                                                                                                                                                  <  1:32863/354‑1 (MQ=255)                                                       tCATCCCTGCGC‑AAACCAGTTGCCCAGTGGTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCca                                                                                                                                                                                                                                                                                       <  1:32864/325‑1 (MQ=255)                                                                                    ggTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑‑ATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCA‑‑‑TGA‑TTGCC‑ATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TT‑ACCCACAGCGactaagcctc                                                                                                                                                                                                                                                                        >  1:32868/1‑296 (MQ=255)                                                                                    ggTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACC‑GCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAGTTT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCC‑ATGCGGG‑‑‑‑‑‑‑TTTTCCTGA‑TTaa                                                                                                                                                                                                                                                                                           >  1:750054/1‑289 (MQ=255)                                                                                     gTAGCGAGAGATCGAGGAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑t                                                                                                                                                                                                                                                                                              >  1:750059/1‑288 (MQ=255)                                                                                     gTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACC‑GCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTG‑CCAGCTTCGACCACAG‑TGAG‑TTCCAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCC‑ATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑tt                                                                                                                                                                                                                                                                                             >  1:750057/1‑286 (MQ=255)                                                                                          gagagaTCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTTGATAGTAATCTGCCGTGTAAT‑‑ATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑‑TGA‑TTGCCA‑‑GCGGG‑‑‑‑‑‑TTTTTCCTGA‑TT‑ACCCACAGCG‑ACTAAAGCCT‑CCCGGTAGAGTgg                                                                                                                                                                                                                                                         >  1:32869/1‑313 (MQ=255)                                                                                                               ttttCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑‑ATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑‑TTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGAtg                                                                                                                                                                                                                                                 >  1:32871/1‑304 (MQ=255)                                                                                                                tttCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACATTTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAATTTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGaaa                                                                                                                                                                                                                                                                               >  1:750060/1‑278 (MQ=255)                                                                                                                  tCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑‑ATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCG‑ACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑C‑TATTTTGCAGGTT‑ACGTTGATCG‑CCCCCCCa                                                                                                                                                                         >  1:32872/1‑368 (MQ=255)                                                                                                                            gCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACC‑GCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTa                                                                                                                                                                                                                                                                             >  1:750062/1‑266 (MQ=255)                                                                                                                                 gCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGgc                                                                                                                                                                                                                                                        >  1:750064/1‑281 (MQ=255)                                                                                                                                                          aTACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGG‑‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGc                                                                                                                                                                                                  <  1:32873/307‑1 (MQ=255)                                                                                                                                                                     cAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAAc                                                                                                                                                                                                                                                                               >  1:750069/1‑224 (MQ=255)                                                                                                                                                                      aGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCTCCCCGGTAGAGTGGCGCTGATGCTGCCGAtgct                                                                                                                                                                                                                                      >  1:750071/1‑264 (MQ=255)                                                                                                                                                                             gCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAAccc                                                                                                                                                                                                                                                                                        >  1:750073/1‑207 (MQ=255)                                                                                                                                                                                   cAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑Ct                                                                                                                                                                                                                         <  1:750068/262‑1 (MQ=255)                                                                                                                                                                                             gCTGTAGACATTTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCAC‑G‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGc                                                                                                                                                                                                                >  1:32877/1‑260 (MQ=255)                                                                                                                                                                                              cTGTAGACA‑TTTTGATAGTAATCTGCCGTG‑AAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGGCGG‑‑‑‑‑‑‑TTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATA‑TTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTg                                                                                                                                     <  1:32876/328‑1 (MQ=255)                                                                                                                                                                                                 tAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGG‑‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGtt                                                                                                                                                                                             <  1:750072/273‑1 (MQ=255)                                                                                                                                                                                                   gACA‑TTTTGATAGTAATCTGCCGTGTAATAAATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTgg                                        >  1:32878/1‑417 (MQ=255)                                                                                                                                                                                                                             taat‑aaTGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGt                                                                                                                                                                                                                                                           >  1:750074/1‑188 (MQ=255)                                                                                                                                                                                                                                    tGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGGTGG‑‑‑‑‑‑‑TTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTg                                                                                                                                                                                      <  1:32879/247‑1 (MQ=255)                                                                                                                                                                                                                                      gCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAg                                                                                                                                                           >  1:750075/1‑272 (MQ=255)                                                                                                                                                                                                                                                       catcatTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAAct                                                    >  1:32882/1‑355 (MQ=255)                                                                                                                                                                                                                                                          catTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAggatcacctgag                                                                                                                                                                                     >  1:750078/1‑218 (MQ=255)                                                                                                                                                                                                                                                          catTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACG‑TGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGc                                                                                                                                        >  1:750079/1‑270 (MQ=255)                                                                                                                                                                                                                                                            ttGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGT‑‑‑‑‑GTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑t                                                                                                                          <  1:32880/282‑1 (MQ=255)                                                                                                                                                                                                                                                                             tgtgCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATG‑GCG‑‑‑‑‑GTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGcc                                                                                                                                                                                                                                            <  1:32881/155‑1 (MQ=255)                                                                                                                                                                                                                                                                               tgCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGGTTTGAGTG‑TGCCGCCGATTTGCACGGTTTGTCGGCATCGCTACGTCGCAA‑TTGATTTGCCAATGCGGGTTTTTTTTTTTCCTGATTTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTg                                                                                                                                                                                                                  >  1:32889/1‑189 (MQ=255)                                                                                                                                                                                                                                                                                 cccAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTa                                                                                                                                                                                            <  1:750077/198‑1 (MQ=255)                                                                                                                                                                                                                                                                                    aGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATG‑GCG‑‑‑‑‑GTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATA‑TTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TATCCCCGTTGCCGCCCAGCTAAGCGTACTc                                                                                       <  1:32885/290‑1 (MQ=255)                                                                                                                                                                                                                                                                                     gCTTCGACCACAGTTGAG‑TT‑CAGG‑TTGAGTGTTGCCGCCGA‑TTGCACGG‑TGTTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATTC‑GG‑‑‑‑‑CGTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAGCCATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACg                                                                                                                                                       <  1:32883/231‑1 (MQ=255)                                                                                                                                                                                                                                                                                        tCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCaa                                                                                                                                                                        >  1:750083/1‑209 (MQ=255)                                                                                                                                                                                                                                                                                          gACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTgccg                                                                                                        >  1:750084/1‑268 (MQ=255)                                                                                                                                                                                                                                                                                             cacaG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGt                                                                                                                                                              <  1:750081/215‑1 (MQ=255)                                                                                                                                                                                                                                                                                               caG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTc                                                                                                                                                                                                                                  >  1:750085/1‑147 (MQ=255)                                                                                                                                                                                                                                                                                                aG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATG‑GCG‑‑‑‑‑GTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATA‑TTTGCAGG‑TAACGTTGATCG‑CCCCCC‑‑AGCGGCAGCGTCAG‑ACGGG‑TTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTgccg                                                                                                        <  1:32886/256‑1 (MQ=255)                                                                                                                                                                                                                                                                                                   tGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAGCCTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAGAACGGGTTTCACCCACCGCCTGATTATCGTACTCCACTAT‑CCCGTTgccg                                                                                                        >  1:750086/1‑263 (MQ=255)                                                                                                                                                                                                                                                                                                   tGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGG‑‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATA‑TTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAg                                                          <  1:32888/303‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     aG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCg                                                                                            <  1:750082/269‑1 (MQ=255)                                                                                                                                                                                                                                                                                                     aG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGc                                                                                >  1:750087/1‑281 (MQ=255)                                                                                                                                                                                                                                                                                                           cAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCa                                                       >  1:32891/1‑302 (MQ=255)                                                                                                                                                                                                                                                                                                                                   tGCACGG‑‑TGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGG‑‑‑‑‑‑‑‑TTTTCCTGA‑‑TAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCGCCCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTA‑CGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTT‑‑CCGTTTTCCGACCAGCGATCCatat  <  1:32890/328‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                                         g‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGGCGG‑‑‑‑‑‑‑TTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTC‑TTTTCCCGTTTTCCGACCAGCGATCCa     <  1:32892/323‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                                              tCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGa                                                                                  <  1:750088/245‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                                                ggCATCGCTACGTCGCAA‑TTGA‑TTGCCAAT‑C‑GG‑‑‑‑‑TGTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCa                                                                                                                                                                                                 <  1:32893/135‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                                                      gCTACGTCGCAA‑TTGA‑TTGCCAATGCGG‑‑‑‑‑‑‑TTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCgg                                                                                                                                                                    <  1:32894/157‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                                                                                 gcgg‑‑‑‑‑TGTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑‑CCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTg                                                              <  1:32895/229‑1 (MQ=255)                                                                                                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                         CATCAAACTGTTTGCGTGCTGACAGGTTTGCCATGGTGTAGCCGTTTTGATGAGTCATCCCTGCGCTAAACCAGTTGCCCAGTGGTAGCGAGAGATCGA‑GAGCGATATATTTTCCTGTATTGGCGTTGCTGTCGCCGTTGTTATAGCGCTGAATACCGGCCCGCAGGCCAAGCGAACCAAAGGTACCGCTGTAGACA‑TTTTGATAGTAATCTGCCGTGTAAT‑AATGGCTGTTGTAACGGCGGTCATCATTGTAGCTGATGCTGAATGTGCCCAGCTTCGACCACAG‑TGAG‑TT‑CAGG‑TTGAGTG‑TGCCGCCGA‑TTGCACGG‑TTGTCGGCATCGCTACGTCGCAA‑TTGA‑TTGCCAATGCGGG‑‑‑‑‑TTTTTTCCTGA‑TTAACCCACAGCGAACTAAAGCCT‑CCCGGTAGAGTGGCGCTGATGCTGCCGATGCTGCTCCATGAG‑CTGTCACTGGCCAG‑CATATTTTGCAGGTTAACGTTGATCG‑CCCCCCCAAGCGGCAGCGTCAG‑ACGGGTTTCACCCACCGCCTGATTATCGTA‑TCCA‑TAT‑CCCGTTGCCGCCCAGCTAAGCGTACTCAATGAGCCGGAGGTCGAGGCACCTGCTAGCCAACTCTCTTTAGCTGGTCGCGTCTTTTTCCCGTTTTCCGACCAGCGATCCATAT  >  NC_000913/307571‑308218

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