Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	USA300TCH1516_ALE	1,302,967	Δ4 bp	100%	coding (325‑328/714 nt)	phnF →	putative transcriptional regulator PhnF

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	USA300TCH1516_ALE	1302967	1302970	4	73 [0]	[1] 74	phnF	putative transcriptional regulator PhnF

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	USA300TCH1516_ALE	= 1302966	0 (0.000)	69 (1.120)	49/256	NT	100%	coding (324/714 nt)	phnF	putative transcriptional regulator PhnF
	?	USA300TCH1516_ALE	1302971 =	0 (0.000)					coding (329/714 nt)	phnF	putative transcriptional regulator PhnF

TTTGAAGAGGGCACAAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAACTATGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACA  >  USA300TCH1516_ALE/1302829‑1303101                                                                                                                                          |                                                                                                                                        tttGAAGAGGGCACAAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATAATTAAAAAAACAGGATTTAAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCATCAtg                                                                                                                                       >  1:59031‑M1/1‑138 (MQ=255)      agagGGCACAAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttaga                                                                                                                                  <  1:67608‑M1/140‑8 (MQ=255)         gggCACAAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgc                                                                                                                               <  1:638443‑M1/140‑11 (MQ=255)         gggCACAAGTGTTAAGTCACAGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTCTTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACCGAATACTTAAATTTAGATGAGCAACCAGCAtgttagaaac                                                                                                                               >  1:443354‑M1/1‑130 (MQ=255)           gCACAAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCATGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaa                                                                                                                             <  2:306854‑M1/140‑13 (MQ=255)             acaAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaat                                                                                                                           >  2:385627‑M1/1‑126 (MQ=255)                aGTGTTAAGTCACTGCCCGCTTCCTTTTATCCATTGAATGAACTTATAAGTATGGGGAAATGGATTAAAAATGCAGGATGTGAATGCGGAACTGAATACTTAAATTTTGATGAGGAACCAGCAtgttagatgcaaatttg                                                                                                                        <  2:459220‑M1/140‑18 (MQ=255)                aGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttg                                                                                                                        <  1:365083‑M1/140‑18 (MQ=255)                 gtgtTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGATCTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgaaaatttgt                                                                                                                       <  2:535962‑M1/140‑19 (MQ=255)                   gtTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttg                                                                                                                     <  1:47093‑M1/140‑21 (MQ=255)                   gtTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttg                                                                                                                     <  2:129269‑M1/140‑21 (MQ=255)                       aGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcg                                                                                                                 <  1:427695‑M1/140‑25 (MQ=255)                            cTGCCCCCTTTCTTTTATCCCTTGAATGAACTTATAAGTATGGGGCAAATGAGTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaa                                                                                                            >  1:145442‑M1/1‑111 (MQ=255)                            cTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaactttgttgagcgttgaa                                                                                                            <  1:363366‑M1/140‑30 (MQ=255)                               cccccTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTNTGATGAGCAACCAGCAtggtagatgcaaatttgttgagcattgaagaa                                                                                                         >  1:614391‑M1/1‑108 (MQ=255)                                ccccTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttaaagaag                                                                                                        >  1:118959‑M1/1‑107 (MQ=255)                                    tttctttTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggata                                                                                                    <  1:156977‑M1/140‑38 (MQ=255)                                    tttctttGATCCATTGAATGGCCTGAGAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggata                                                                                                    <  1:163187‑M1/140‑38 (MQ=255)                                     ttttgttatccaTTGAATGAACTTAGAAGTATTGTGAAAATGATTAAAAATGCAGGATTGGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatat                                                                                                   <  1:359837‑M1/135‑39 (MQ=255)                                     ttctttTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatat                                                                                                   <  1:94671‑M1/140‑39 (MQ=255)                                     ttctttTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatat                                                                                                   <  1:641096‑M1/140‑39 (MQ=255)                                        ttttATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatcca                                                                                                >  2:252553‑M1/1‑99 (MQ=255)                                        ttttATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACAGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatatgttgagcgttgaagaaggatagcca                                                                                                >  2:433211‑M1/1‑99 (MQ=255)                                         tttATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAATTTTGTTTTACCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccag                                                                                               >  1:345778‑M1/1‑98 (MQ=255)                                           tATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccggta                                                                                             >  1:91321‑M1/1‑96 (MQ=255)                                               cATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAAATTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaacca                                                                                         >  1:265162‑M1/1‑92 (MQ=255)                                               cATTGAATGAACTTACAAGCATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaacca                                                                                         >  1:541617‑M1/1‑92 (MQ=255)                                               cATTGAATGAAATTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAAAGCGGAAATGAATACTTAAATTTTGATGAGCAACAAGCAtgttagatgcaaatttgttgcgcgttgaagaaggatatacagtaaaca                                                                                         >  2:176615‑M1/1‑92 (MQ=255)                                               cATTGAATCAACTTAGAATTATTAGGCAAATGATTAAAAATGCAGCATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgtggagcgttgaagaaggatatccagtaacca                                                                                         <  2:149228‑M1/140‑49 (MQ=255)                                                 ttGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGGTTTGAATGCGGAACTGAATACTTAAATTTTGATGAACAACCAGCAtgttaaatgcaaatttgttgagcgttgaagaaggatatccattaaccatt                                                                                       >  1:694921‑M1/1‑90 (MQ=255)                                                 ttGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggaaatccagtaaccatt                                                                                       >  1:25410‑M1/1‑90 (MQ=255)                                                     atgaaCTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatag                                                                                   <  1:131648‑M1/140‑55 (MQ=255)                                                        aaCTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaac                                                                                >  1:490771‑M1/1‑83 (MQ=255)                                                        aaCTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaac                                                                                >  1:273809‑M1/1‑83 (MQ=255)                                                        aaCTTATAAGTATTGGGAAAATGATTAAAAATGCAGGATTTGAATGCGGAACGGAATACTTAAATTTTGATGAGCAACCATCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaac                                                                                <  2:684927‑M1/140‑58 (MQ=255)                                                         aCTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacg                                                                               >  1:273202‑M1/1‑82 (MQ=255)                                                                    ttGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattctagaacgattacgaactg                                                                    >  2:47637‑M1/1‑71 (MQ=255)                                                                    ttGGGCAAATGATCAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactg                                                                    <  2:167623‑M1/140‑70 (MQ=255)                                                                        gCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccga                                                                >  2:726123‑M1/1‑67 (MQ=255)                                                                        gCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatataatgagcgttgaagaaggatatccagtaacaattatagaacggttacgaactgccga                                                                >  1:137953‑M1/1‑67 (MQ=255)                                                                            aTGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattccgaactgccaaggga                                                            >  2:46214‑M1/1‑63 (MQ=255)                                                                             tgaatcaaaattcaggATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggag                                                           <  1:163148‑M1/134‑79 (MQ=255)                                                                             tGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgatagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggag                                                           >  1:665115‑M1/1‑62 (MQ=255)                                                                              gATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgtaagatgcaaatttgttgagcgttgaagaaggatatccattaaccattatagaacgattacaaactagcgatggaaa                                                          >  2:414398‑M1/1‑61 (MQ=255)                                                                               atcaaaaatgcAGGATTTGAATGCGGAACTGAATACTTAAATTGTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaataaggatatccagtaaccatcatagaccgcttactaactgccgatggagaa                                                         <  2:156410‑M1/137‑81 (MQ=255)                                                                               aataaAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaacgagatccagtaaccattatagaacgattacgaactgccgatggagaa                                                         <  2:298181‑M1/138‑81 (MQ=255)                                                                               aTTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAGATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaa                                                         <  1:513801‑M1/140‑81 (MQ=255)                                                                               aTTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatgtggaa                                                         >  2:138099‑M1/1‑60 (MQ=255)                                                                                ttAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaac                                                        >  1:611905‑M1/1‑59 (MQ=255)                                                                                 tAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatagccagtaaccattatagaacgattacgaactgccgatggagaacc                                                       <  1:383969‑M1/140‑83 (MQ=255)                                                                                      aTGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcg                                                  <  2:677921‑M1/140‑88 (MQ=255)                                                                                         cAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagcaccggtcgtct                                               >  2:184138‑M1/1‑50 (MQ=255)                                                                                         cAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagacccggtcgtct                                               >  1:726130‑M1/1‑50 (MQ=255)                                                                                                 gAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAtggtagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgttta                                       <  1:376834‑M1/140‑99 (MQ=255)                                                                                                   aTGCGGAACTGAATACTTNTATTTTGATGAGCAACCTGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtagtcattcgaaaactattacgaactgccgatggagaaccggtcgtctattgtttaga                                     <  1:646078‑M1/140‑101 (MQ=255)                                                                                                      cGGAATTGAATAGTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccatcatagaacgattacgaactgccgagggagaaccggtcgtctattgtttatataa                                  <  2:443354‑M1/140‑104 (MQ=255)                                                                                                         aaCTGAATACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaat                               <  2:577508‑M1/140‑107 (MQ=255)                                                                                                             gAATTCTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgct                           <  1:101218‑M1/140‑111 (MQ=255)                                                                                                              aaTACTTAAATTTTGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgcta                          >  2:393485‑M1/1‑29 (MQ=255)                                                                                                                  cTTAAATTTTTATGAGCAACCAGCAtgttagatgcaaattggttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaa                      <  2:232961‑M1/140‑116 (MQ=255)                                                                                                                         tttGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgacgaaggctatccagtcaccattatagaacgattacgaactgatgatgtacaaccgttcgtcaattgtttagataaaattgcttaaatagaatta               >  1:275559‑M1/1‑18 (MQ=255)                                                                                                                         tttGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaatta               >  2:493041‑M1/1‑18 (MQ=255)                                                                                                                         tttGATGAGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattaagaattgtcttagaagaaccggtcgtctattgtttagataaaattgctaaaaaagaatta               >  1:680577‑M1/1‑18 (MQ=255)                                                                                                                          ctgatggcCAACCCGCGtgttagatgtatatttgtttagcgttgaagaatgatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaa              <  2:498291‑M1/139‑124 (MQ=255)                                                                                                                            gatgaGTAACGAGCAtgttagatggaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattccgacctgtcgaaggagaaccggtcgtctattgtttagataaaataaataaaaaagaattaaca            <  2:545689‑M1/140‑126 (MQ=255)                                                                                                                            gatgaGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaacagaattaaca            >  2:428784‑M1/1‑15 (MQ=255)                                                                                                                            gatgaGCAACCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaaca            <  2:345778‑M1/140‑126 (MQ=255)                                                                                                                            gatgaGCAAACAGCAtgtcagatgcgaagtgtttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaaca            <  2:74862‑M1/140‑126 (MQ=255)                                                                                                                            aagcagcaaCCAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctgagctatatttaaaattgatataaaagaattaaca            <  2:666148‑M1/136‑126 (MQ=255)                                                                                                                                    aCCAGCAtgttcgatgcacatttgttgcgcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaacatgtacagg    >  1:661033‑M1/1‑7 (MQ=255)                                                                                                                                     ccAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaacatgtacagag   >  1:451850‑M1/1‑6 (MQ=255)                                                                                                                                     ccAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaaaatgtacagag   >  1:451950‑M1/1‑6 (MQ=255)                                                                                                                                      cAGCAtgttagatgcaaatttgttgagcgttgaagaaggatatccagtaaccattatagaacgattacgaactgccgatggagaaccggtcgtctattgtttagataaaattgctaaaaaagaattaacatgtacagagt  <  1:726123‑M1/140‑136 (MQ=255)                                                                                                                                          |                                                                                                                                        TTTGAAGAGGGCACAAGTGTTAAGTCACTGCCCCCTTTCTTTTATCCATTGAATGAACTTATAAGTATTGGGCAAATGATTAAAAATGCAGGATTTGAATGCGGAACTGAATACTTAAATTTTGATGAGCAACCAGCAACTATGTTAGATGCAAATTTGTTGAGCGTTGAAGAAGGATATCCAGTAACCATTATAGAACGATTACGAACTGCCGATGGAGAACCGGTCGTCTATTGTTTAGATAAAATTGCTAAAAAAGAATTAACATGTACA  >  USA300TCH1516_ALE/1302829‑1303101

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