Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,596,220	(TCA)3→2	coding (182‑184/927 nt)	livH ←	branched‑chain amino acid ABC transporter permease

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3596212	3596214	3	97 [0]	[0] 98	livH	branched‑chain amino acid ABC transporter permease

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3596215 =	0 (0.000)	87 (0.700)	63/182	1.4	100%	coding (189/927 nt)	livH	branched‑chain amino acid ABC transporter permease
	?	NC_000913	= 3596211	0 (0.000)					coding (193/927 nt)	livH	branched‑chain amino acid ABC transporter permease

TAAGCCACCCGTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATG  >  NC_000913/3596117‑3596314                                                                                                   |                                                                                                    cggtaagccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccca                                                                                                    >  2:2242303‑M1/99‑100 (MQ=255) cggtaagccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccca                                                                                                    <  1:2134465‑M1/2‑1 (MQ=255)  ggtaagccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgcccat                                                                                                   >  1:1046422‑M1/98‑100 (MQ=255)   gtaagccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgcccatc                                                                                                  <  1:1348327‑M1/4‑1 (MQ=255)    taagccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgcccatca                                                                                                 >  1:2553114‑M1/96‑100 (MQ=255)     aagccacccgttcgatactccagccttaggcgctggcaatgacgattgcgccgacgaagcccgcagctaccagcagccagccggtatcaatgcccatcat                                                                                                >  2:2985989‑M1/95‑100 (MQ=255)      agccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCatcatc                                                                                               >  2:250092‑M1/94‑100 (MQ=255)       gccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCAtcatca                                                                                              <  1:2982179‑M1/8‑1 (MQ=255)       gccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCAtcatca                                                                                              <  2:2170710‑M1/8‑1 (MQ=255)        ccacccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAg                                                                                             <  1:1990631‑M1/9‑1 (MQ=255)          acccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAgcg                                                                                           <  1:3020450‑M1/11‑1 (MQ=255)           cccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAgcgc                                                                                          <  2:269597‑M1/12‑1 (MQ=255)            ccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGcgcg                                                                                         <  1:1304503‑M1/13‑1 (MQ=255)            ccgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGcgcg                                                                                         <  2:2779479‑M1/13‑1 (MQ=255)             cgttcgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCgcgg                                                                                        >  1:1344337‑M1/87‑100 (MQ=255)                 cgatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCgat                                                                                    <  1:2515392‑M1/18‑1 (MQ=255)                  gatactccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCgatg                                                                                   >  2:2915917‑M1/82‑100 (MQ=255)                      ctccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCa                                                                               <  2:2662724‑M1/23‑1 (MQ=255)                      ctccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCa                                                                               <  2:1442606‑M1/23‑1 (MQ=255)                      caccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagctgccagccggtatcaatgccCATCATCAGCGCGGCgatgat                                                                                 >  2:906090‑M1/78‑98 (MQ=255)                        ccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagtagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATa                                                                             <  1:1681342‑M1/25‑1 (MQ=255)                        ccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATa                                                                             >  1:2651510‑M1/76‑100 (MQ=255)                        ccagccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATa                                                                             <  1:3115745‑M1/25‑1 (MQ=255)                          agccgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATaaa                                                                           >  1:2575739‑M1/74‑100 (MQ=255)                             cgtaggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATga                                                                        >  2:177312‑M1/71‑100 (MQ=255)                               taggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATgaga                                                                      <  2:2370749‑M1/32‑1 (MQ=255)                                ttccgatctcaatgacgattgcgccgacgaatcccgctgctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGAc                                                                     >  2:1297541‑M1/68‑100 (MQ=255)                                 ggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACg                                                                    <  2:3035998‑M1/34‑1 (MQ=255)                                 ggcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACg                                                                    <  2:1434348‑M1/34‑1 (MQ=255)                                  gcgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGt                                                                   >  2:774644‑M1/66‑100 (MQ=255)                                   cgctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGt                                                                   >  2:2652056‑M1/65‑99 (MQ=255)                                    gctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAg                                                                 <  2:1969803‑M1/37‑1 (MQ=255)                                     ctggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGc                                                                <  1:3317349‑M1/38‑1 (MQ=255)                                       ggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTg                                                              <  2:2247826‑M1/40‑1 (MQ=255)                                       ggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTg                                                              <  2:1217627‑M1/40‑1 (MQ=255)                                       ggcaatgacgattgcgccgacgaatcccgcagctaccagcagccagccgggatcaatgccCATCATCAGCGCGGCGACGATCATAAATGAGACGTAGCt                                                               <  2:1626327‑M1/39‑1 (MQ=255)                                         caatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCAaaaatgagacagatcggaa                                                            <  1:1297545‑M1/42‑10 (MQ=255)                                         caatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGcc                                                            <  1:80090‑M1/42‑1 (MQ=255)                                         caatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGcc                                                            >  2:1508815‑M1/59‑100 (MQ=255)                                         caatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGcc                                                            >  2:1086558‑M1/59‑100 (MQ=255)                                          aatgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCa                                                           <  2:2243077‑M1/43‑1 (MQ=255)                                           atgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCaa                                                          <  1:1340182‑M1/44‑1 (MQ=255)                                           atgacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCaa                                                          <  2:2942897‑M1/44‑1 (MQ=255)                                             gacgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATc                                                        >  1:2928700‑M1/55‑100 (MQ=255)                                              acgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCa                                                       >  1:1067717‑M1/54‑100 (MQ=255)                                              acgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCa                                                       >  1:1109474‑M1/54‑100 (MQ=255)                                              acgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCa                                                       >  2:1510573‑M1/54‑100 (MQ=255)                                               cgattgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCat                                                      >  1:228675‑M1/53‑100 (MQ=255)                                                 attgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCatat                                                    >  1:3124246‑M1/51‑100 (MQ=255)                                                  ttgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCatat                                                    >  2:3200685‑M1/50‑99 (MQ=255)                                                  ttgcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCAtata                                                   >  1:3191415‑M1/50‑100 (MQ=255)                                                    gcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATaaa                                                 >  1:2170721‑M1/48‑100 (MQ=255)                                                    gcgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATaaa                                                 <  2:755712‑M1/53‑1 (MQ=255)                                                     cgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAAc                                                >  2:2238727‑M1/47‑100 (MQ=255)                                                     cgccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAAAGAGACGTAGATGCCAATCATATaaa                                                 <  2:1937700‑M1/53‑1 (MQ=255)                                                       ccgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCt                                              <  2:3147615‑M1/56‑1 (MQ=255)                                                        cgacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTc                                             <  2:1030411‑M1/57‑1 (MQ=255)                                                         gacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCg                                            <  2:2642554‑M1/58‑1 (MQ=255)                                                         gacgaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCg                                            <  1:3257825‑M1/58‑1 (MQ=255)                                                            gaatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCg                                         <  2:1598320‑M1/61‑1 (MQ=255)                                                             aatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGt                                        <  2:55915‑M1/62‑1 (MQ=255)                                                             aatcccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGt                                        <  2:2268917‑M1/62‑1 (MQ=255)                                                                cccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTggg                                     >  2:24126‑M1/36‑100 (MQ=255)                                                                 ccgcagctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGc                                    >  1:1001109‑M1/35‑100 (MQ=255)                                                                     agctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAg                                >  2:1278377‑M1/31‑100 (MQ=255)                                                                      gctaccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGaa                                 >  2:2510176‑M1/30‑98 (MQ=255)                                                                      gctaccagcagccagccggtatcaatgccCACCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGt                               <  1:2254034‑M1/71‑1 (MQ=255)                                                                        taccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTg                             >  1:1149393‑M1/28‑100 (MQ=255)                                                                        taccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTg                             <  2:1448628‑M1/73‑1 (MQ=255)                                                                        taccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAA‑TGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGa                            <  2:310197‑M1/73‑1 (MQ=255)                                                                          ccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGAt                           >  1:1040064‑M1/26‑100 (MQ=255)                                                                          ccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGAt                           >  1:830927‑M1/26‑100 (MQ=255)                                                                          ccagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGAt                           <  1:829997‑M1/75‑1 (MQ=255)                                                                           cagcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAACGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATc                          >  2:838346‑M1/25‑100 (MQ=255)                                                                            agcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCa                         >  2:2571728‑M1/24‑100 (MQ=255)                                                                             gcagccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCAt                        <  2:380989‑M1/78‑1 (MQ=255)                                                                                gccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGcc                     <  1:131076‑M1/81‑1 (MQ=255)                                                                                 ccagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCg                    <  2:359694‑M1/82‑1 (MQ=255)                                                                                  cagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGa                   <  2:1482884‑M1/83‑1 (MQ=255)                                                                                  cagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGa                   >  2:1304499‑M1/18‑100 (MQ=255)                                                                                  cagccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGa                   <  1:1450223‑M1/83‑1 (MQ=255)                                                                                   agccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGAt                  >  2:3317331‑M1/17‑100 (MQ=255)                                                                                     ccggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATaa                <  2:3016004‑M1/86‑1 (MQ=255)                                                                                       ggtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATg              >  2:475222‑M1/13‑100 (MQ=255)                                                                                        gtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGc             >  2:3020433‑M1/12‑100 (MQ=255)                                                                                        gtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGc             >  1:1217631‑M1/12‑100 (MQ=255)                                                                                        gtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGc             <  1:1294655‑M1/89‑1 (MQ=255)                                                                                        gtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGc             >  1:2662737‑M1/12‑100 (MQ=255)                                                                                        gtatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGc             >  2:1340178‑M1/12‑100 (MQ=255)                                                                                         tatcaatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGcc            <  1:1886706‑M1/90‑1 (MQ=255)                                                                                            caatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTa         <  1:1545460‑M1/93‑1 (MQ=255)                                                                                             aatgccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTaa        <  1:2414682‑M1/94‑1 (MQ=255)                                                                                                gccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAAcc     >  2:2515380‑M1/4‑100 (MQ=255)                                                                                                 ccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCa    <  1:2136016‑M1/98‑1 (MQ=255)                                                                                                 ccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCa    <  1:2459230‑M1/98‑1 (MQ=255)                                                                                                 ccCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCa    <  2:1556952‑M1/98‑1 (MQ=255)                                                                                                   CATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATg  >  2:2419307‑M1/1‑100 (MQ=255)                                                                                                   CATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATg  <  1:2919327‑M1/100‑1 (MQ=255)                                                                                                   |                                                                                                    TAAGCCACCCGTTCGATACTCCAGCCGTAGGCGCTGGCAATGACGATTGCGCCGACGAATCCCGCAGCTACCAGCAGCCAGCCGGTATCAATGCCCATCATCATCAGCGCGGCGATGATCATAAATGAGACGTAGCTGCCAATCATATAAACCTCGCCGTGGGCGAAGTTGATCATGCCGATAATGCCGTAAACCATG  >  NC_000913/3596117‑3596314

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