Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	257,908	Δ776 bp	100%		insB9–[crl]	insB9, insA9, [crl]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	257908–258676	258683	8–776	92 [0]	[0] 94	insB9–[crl]	insB9,insA9,[crl]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 257907	0 (0.000)	74 (1.170)	53/160	NT	100%	intergenic (+8/+16)	crl/insB9	RNA polymerase holoenzyme assembly factor Crl/IS1 transposase B
	?	NC_000913	258684 =	0 (0.000)					pseudogene (9/331 nt)	crl	RNA polymerase holoenzyme assembly factor Crl

CATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTT  >  NC_000913/257809‑258003                                                                                                   |                                                                                                 cATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA                                                                                                  <  1:947917‑M1/99‑1 (MQ=255) cATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTA                                                                                                  >  2:766791‑M1/1‑99 (MQ=255)  aTCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAa                                                                                                 <  2:1486489‑M1/99‑2 (MQ=255)   tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAag                                                                                                <  2:96822‑M1/99‑3 (MQ=255)   tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAag                                                                                                >  1:741083‑M1/1‑97 (MQ=255)   tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAag                                                                                                >  2:928592‑M1/1‑97 (MQ=255)   tCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAag                                                                                                >  2:404882‑M1/1‑97 (MQ=255)    cacaACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                               <  1:337601‑M1/99‑4 (MQ=255)    cacaACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagt                                                                                               <  1:1717652‑M1/99‑4 (MQ=255)      caacaGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgc                                                                                             <  1:1990526‑M1/99‑6 (MQ=255)         caGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaa                                                                                          <  2:1685843‑M1/99‑9 (MQ=255)            gagaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATCGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagat                                                                                       <  2:916052‑M1/99‑12 (MQ=255)             agaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagata                                                                                      <  2:1434387‑M1/99‑13 (MQ=255)              gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                                                     >  1:1887185‑M1/1‑86 (MQ=255)              gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                                                     >  1:1975656‑M1/1‑86 (MQ=255)              gaTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataa                                                                                     >  1:1978649‑M1/1‑86 (MQ=255)               aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                                                                    <  2:1993802‑M1/99‑15 (MQ=255)               aTAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataat                                                                                    <  2:208796‑M1/99‑15 (MQ=255)                tAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatc                                                                                   >  2:1420354‑M1/1‑84 (MQ=255)                  gCAATGACGTTACCGAGTGGACACGAGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCTGTATATTCGTGAATGTAagtgcaaagatgatcga                                                                                 <  2:365129‑M1/99‑18 (MQ=255)                  gCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcga                                                                                 >  2:428708‑M1/1‑82 (MQ=255)                   cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                <  1:1383476‑M1/99‑19 (MQ=255)                   cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                <  2:1558571‑M1/99‑19 (MQ=255)                   cAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgat                                                                                <  1:1625329‑M1/99‑19 (MQ=255)                    aaTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgatt                                                                               <  2:754212‑M1/99‑20 (MQ=255)                     aTGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattc                                                                              >  1:965829‑M1/1‑79 (MQ=255)                      tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                                                                             <  2:1946839‑M1/99‑22 (MQ=255)                      tGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattct                                                                             >  1:910914‑M1/1‑78 (MQ=255)                        aCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttt                                                                           <  2:1174270‑M1/99‑24 (MQ=255)                        aCGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttt                                                                           <  2:1492895‑M1/99‑24 (MQ=255)                         cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                                                                          <  1:897118‑M1/99‑25 (MQ=255)                         cGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttt                                                                          <  2:93352‑M1/99‑25 (MQ=255)                          gTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattcttttt                                                                         >  2:575018‑M1/1‑74 (MQ=255)                           ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                                                                        >  1:426583‑M1/1‑73 (MQ=255)                           ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                                                                        >  2:1587254‑M1/1‑73 (MQ=255)                           ttACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGGATATTCGTGAAGGTAagtgcaaagataatcgattctttttc                                                                        >  2:1885399‑M1/1‑73 (MQ=255)                             aCCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcga                                                                      >  1:243030‑M1/1‑71 (MQ=255)                              ccGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgat                                                                     <  1:109172‑M1/99‑30 (MQ=255)                               cGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgatt                                                                    >  1:1731657‑M1/1‑69 (MQ=255)                                 aGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgt                                                                  >  2:413263‑M1/1‑67 (MQ=255)                                  gTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtc                                                                 <  1:683459‑M1/99‑34 (MQ=255)                                   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct                                                                <  2:784143‑M1/99‑35 (MQ=255)                                   tGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtct                                                                <  1:1977071‑M1/99‑35 (MQ=255)                                      acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggc                                                             >  2:695688‑M1/1‑62 (MQ=255)                                      acacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggc                                                             >  1:453632‑M1/1‑62 (MQ=255)                                       cacCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggct                                                            <  2:1461858‑M1/99‑39 (MQ=255)                                        acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                           <  1:321904‑M1/99‑40 (MQ=255)                                        acCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctg                                                           >  2:1324603‑M1/1‑60 (MQ=255)                                          ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgta                                                         <  1:230589‑M1/99‑42 (MQ=255)                                          ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgta                                                         <  1:1875335‑M1/99‑42 (MQ=255)                                          ccGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgta                                                         <  1:467394‑M1/99‑42 (MQ=255)                                            gAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatg                                                       >  1:933298‑M1/1‑56 (MQ=255)                                             aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                      <  1:1269665‑M1/99‑45 (MQ=255)                                             aaGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgc                                                      <  1:2023560‑M1/99‑45 (MQ=255)                                                agCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtc                                                   >  2:925604‑M1/1‑52 (MQ=255)                                                 gCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                                                  >  1:1646321‑M1/1‑51 (MQ=255)                                                 gCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtca                                                  >  2:1745797‑M1/1‑51 (MQ=255)                                                  cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                                                 >  2:1452555‑M1/1‑50 (MQ=255)                                                  cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                                                 <  1:144164‑M1/99‑50 (MQ=255)                                                  cAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaa                                                 >  2:936779‑M1/1‑50 (MQ=255)                                                    gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg                                               >  1:208759‑M1/1‑48 (MQ=255)                                                    gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg                                               >  2:619780‑M1/1‑48 (MQ=255)                                                    gATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacg                                               >  2:112351‑M1/1‑48 (MQ=255)                                                        gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaa                                           >  1:478556‑M1/1‑44 (MQ=255)                                                        gATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaa                                           >  2:639005‑M1/1‑44 (MQ=255)                                                          tCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaac                                         <  2:1169852‑M1/99‑58 (MQ=255)                                                            aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccg                                       <  1:108707‑M1/99‑60 (MQ=255)                                                            aaaaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccg                                       >  1:355984‑M1/1‑40 (MQ=255)                                                              aaaaTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggc                                     >  1:475895‑M1/1‑38 (MQ=255)                                                                 aTTTACCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcacc                                  >  1:1675566‑M1/1‑35 (MQ=255)                                                                     aCCGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaa                              <  2:449196‑M1/99‑69 (MQ=255)                                                                       cGCACTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagt                            >  1:385244‑M1/1‑29 (MQ=255)                                                                          aCTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcg                         >  2:360078‑M1/1‑26 (MQ=255)                                                                           cTAGGCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgt                        >  1:203790‑M1/1‑25 (MQ=255)                                                                              ggCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaa                     >  2:1780279‑M1/1‑22 (MQ=255)                                                                               gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaat                    >  1:776073‑M1/1‑21 (MQ=255)                                                                               gCCCGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaat                    >  2:1268502‑M1/1‑21 (MQ=255)                                                                                cccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaatt                   >  1:534430‑M1/1‑20 (MQ=255)                                                                                 ccGTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattc                  >  1:1685709‑M1/1‑19 (MQ=255)                                                                                   gTATATTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctg                <  2:1243779‑M1/99‑83 (MQ=255)                                                                                    tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg               <  2:1966143‑M1/99‑84 (MQ=255)                                                                                    tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg               <  2:1776526‑M1/99‑84 (MQ=255)                                                                                    tataTTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgg               >  1:681668‑M1/1‑16 (MQ=255)                                                                                      tatTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctgggg             >  2:228333‑M1/1‑14 (MQ=255)                                                                                       atTCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggc            <  1:619791‑M1/99‑87 (MQ=255)                                                                                        ttCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggct           >  2:696155‑M1/1‑12 (MQ=255)                                                                                        ttCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggct           >  2:711654‑M1/1‑12 (MQ=255)                                                                                        ttCGTGAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggct           >  2:1846298‑M1/1‑12 (MQ=255)                                                                                        ttCGTGAAGGTAagtgcaaagataatcgattctttctcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggct           >  2:1206378‑M1/1‑12 (MQ=255)                                                                                             gAAGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgg      >  1:139354‑M1/1‑7 (MQ=255)                                                                                              aaGGTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtgga     <  2:333127‑M1/99‑94 (MQ=255)                                                                                                 gTAagtgcaaagataatcgattctttttcgattgtctggctgtatgcgtcaacgtgaaaccggcaccggaagtgcgtgaattctggggctggtggatgg  <  2:1566273‑M1/99‑97 (MQ=255)                                                                                                   |                                                                                                 CATCACAACAGGAGATAGCAATGACGTTACCGAGTGGACACCCGAAGAGCAGATTGATCAAAAAATTTACCGCACTAGGCCCGTATATTCGTGAAGGTAGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTT  >  NC_000913/257809‑258003

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