BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,754,310	C→A	R382M (AGG→ATG)	yiaY ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,754,310	0	C	A	99.7%	3053.7 / NA	764	R382M (AGG→ATG)	yiaY	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (0/2); new base A (406/355); total (407/357)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 778 total aligned reads displayed.

TACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGCATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACC  >  NC_000913/3754217‑3754402

 aCGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCg                                                                                           <  1:5340914/96‑1 (MQ=255)

 aCGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCg                                                                                           >  2:8662410/1‑96 (MQ=255)

     cAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                       >  1:4890709/1‑96 (MQ=255)

      aGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCg                                                                                      <  1:9869608‑M1/96‑2 (MQ=255)

      aGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCg                                                                                      <  2:11733219‑M1/96‑2 (MQ=255)

       gCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgc                                                                                     <  1:6781479‑M1/96‑3 (MQ=255)

        cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcg                                                                                    <  2:8892083‑M1/96‑4 (MQ=255)

        cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcg                                                                                    >  1:2058046‑M1/1‑93 (MQ=255)

        cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcg                                                                                    >  1:5825722‑M1/1‑93 (MQ=255)

           tCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgata                                                                                 <  1:11263101‑M1/96‑7 (MQ=255)

            cACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataa                                                                                >  1:7215600‑M1/1‑89 (MQ=255)

               gAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                             <  1:2283257‑M1/96‑11 (MQ=255)

               gAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                             <  2:6745803‑M1/96‑11 (MQ=255)

               gAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                             <  1:14242456‑M1/96‑11 (MQ=255)

               gAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                             <  2:3505908‑M1/96‑11 (MQ=255)

               gAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                             <  2:12445393‑M1/96‑11 (MQ=255)

               gAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                             >  2:7619618‑M1/1‑86 (MQ=255)

                   ttACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgcca                                                                         >  2:5150310‑M1/1‑82 (MQ=255)

                    tACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccac                                                                        >  2:6994904‑M1/1‑81 (MQ=255)

                     aCGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccaca                                                                       <  2:7845776‑M1/96‑17 (MQ=255)

                      cGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaa                                                                      >  1:5540251‑M1/1‑79 (MQ=255)

                      cGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGACGTTGTGGAAATGATGATTACATCGCTGAgcgataaatcgccaaaa                                                                      >  2:7430296‑M1/1‑79 (MQ=255)

                        gccaaGGGTTTCAGACATAGTTGAAAAATCGGCTAAGATTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatc                                                                    <  1:13761313‑M1/94‑20 (MQ=255)

                         gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaattt                                                                   >  1:653524‑M1/1‑76 (MQ=255)

                         gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaattt                                                                   >  2:2758924‑M1/1‑76 (MQ=255)

                         gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaattt                                                                   >  2:8404968‑M1/1‑76 (MQ=255)

                         gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGGTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaattt                                                                   >  1:1017659‑M1/1‑76 (MQ=255)

                            aGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttctt                                                                >  2:6495266‑M1/1‑73 (MQ=255)

                            aGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttctt                                                                <  2:3975743‑M1/96‑24 (MQ=255)

                             gggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGGTTACATCGCTGCgcgataaatcgccacaatttcttc                                                               >  1:12710669‑M1/1‑72 (MQ=255)

                             gggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGGGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttctcc                                                               >  2:7258725‑M1/1‑72 (MQ=255)

                              ggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcg                                                              <  2:5495680‑M1/96‑26 (MQ=255)

                              ggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcg                                                              >  1:13478280‑M1/1‑71 (MQ=255)

                               gTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgt                                                             <  1:77447‑M1/96‑27 (MQ=255)

                                 tACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtga                                                           >  2:2835253‑M1/1‑68 (MQ=255)

                                 tACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtga                                                           >  1:13189708‑M1/1‑68 (MQ=255)

                                  aCAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgag                                                          <  1:4547859‑M1/96‑30 (MQ=255)

                                   cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagt                                                         <  1:1459151‑M1/96‑31 (MQ=255)

                                    agagATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagtt                                                        <  1:9491633‑M1/96‑32 (MQ=255)

                                      agaTAGTTGAAAAAGCGGCTAACAATTTACCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgc                                                      >  2:13916945‑M1/1‑63 (MQ=255)

                                       gaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcc                                                     >  1:12500535‑M1/1‑62 (MQ=255)

                                        aTAGTTAAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcct                                                    <  1:3098611‑M1/96‑36 (MQ=255)

                                         tAGTTGAAAAAGCGGCTACCAATGTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                   <  2:12606431‑M1/96‑37 (MQ=255)

                                         tAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                   <  2:13478280‑M1/96‑37 (MQ=255)

                                          ggtTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgcgttgcctgg                                                  <  2:2085892‑M1/95‑38 (MQ=255)

                                           gTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgga                                                 <  1:10347972‑M1/96‑39 (MQ=255)

                                            ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAGTGATGATTACATCGCTGCgcgataaatcgccacgatttcttcgtgagttgcctgatc                                                >  2:10144209‑M1/1‑57 (MQ=255)

                                            ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggat                                                <  1:2853596‑M1/96‑40 (MQ=255)

                                             tGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTAGATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatc                                               >  1:3748114‑M1/1‑56 (MQ=255)

                                             tGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatc                                               <  2:5273138‑M1/96‑41 (MQ=255)

                                             tGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatc                                               <  2:1078718‑M1/96‑41 (MQ=255)

                                              gAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcg                                              >  2:1101073‑M1/1‑55 (MQ=255)

                                              gAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcg                                              >  1:182200‑M1/1‑55 (MQ=255)

                                               aaaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgg                                             <  1:12272477‑M1/96‑43 (MQ=255)

                                               aaaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgg                                             <  1:6989600‑M1/96‑43 (MQ=255)

                                                aaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcggg                                            >  2:9077185‑M1/1‑53 (MQ=255)

                                                aaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcggg                                            >  2:945186‑M1/1‑53 (MQ=255)

                                                 aaaGCGGCTGACAATTTGCCAGCCGTTGTGGGAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggt                                           >  1:9090464‑M1/1‑52 (MQ=255)

                                                  aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggtt                                          >  1:8160505‑M1/1‑51 (MQ=255)

                                                  aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggtt                                          <  2:13742741‑M1/96‑46 (MQ=255)

                                                   aGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggtta                                         <  2:3841140‑M1/96‑47 (MQ=255)

                                                      ggCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagta                                      >  1:11954646‑M1/1‑47 (MQ=255)

                                                      ggCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagta                                      >  2:10370825‑M1/1‑47 (MQ=255)

                                                       gCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaa                                     >  2:2367531‑M1/1‑46 (MQ=255)

                                                            cAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcca                                >  2:1149017‑M1/1‑41 (MQ=255)

                                                            cAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagcca                                >  2:9260985‑M1/1‑41 (MQ=255)

                                                              aTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccaca                              <  2:5916975‑M1/96‑58 (MQ=255)

                                                              aTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccaca                              <  2:13338655‑M1/96‑58 (MQ=255)

                                                                  gCCAGTCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggagcgggttagtaaagccacaggaa                          <  1:11149700‑M1/96‑62 (MQ=255)

                                                                  gCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca                          >  1:4252902‑M1/1‑35 (MQ=255)

                                                                   ccAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcat                         >  1:10024068‑M1/1‑34 (MQ=255)

                                                                    cAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc                        >  1:9388102‑M1/1‑33 (MQ=255)

                                                                    cAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc                        >  1:7607479‑M1/1‑33 (MQ=255)

                                                                    cAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc                        >  2:6111320‑M1/1‑33 (MQ=255)

                                                                      gCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                      >  2:11221547‑M1/1‑31 (MQ=255)

                                                                       ccGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatcttt                     >  1:345603‑M1/1‑30 (MQ=255)

                                                                           tgtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagg                 <  1:8480091‑M1/96‑71 (MQ=255)

                                                                           tgtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcacg                 <  1:8068549‑M1/96‑71 (MQ=255)

                                                                            gtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcaggg                >  1:1632005‑M1/1‑25 (MQ=255)

                                                                            gtgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatatttcagag                <  1:10498158‑M1/96‑72 (MQ=255)

                                                                             tgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggctagggttagaaaagccacaggcatctttcagggc               >  1:11491381‑M1/1‑24 (MQ=255)

                                                                             tgGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggc               <  1:12166121‑M1/96‑73 (MQ=255)

                                                                              ggAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggca              >  2:11960722‑M1/1‑23 (MQ=255)

                                                                                aaaTGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt            >  2:4587311‑M1/1‑21 (MQ=255)

                                                                                aaaTGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt            >  1:2718368‑M1/1‑21 (MQ=255)

                                                                                aaaTGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt            <  2:4254462‑M1/96‑76 (MQ=255)

                                                                                 aaTGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattc           <  2:14363301‑M1/96‑77 (MQ=255)

                                                                                 aaTGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatta           >  1:6212013‑M1/1‑20 (MQ=255)

                                                                                  atgatgATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattag          >  1:8703728‑M1/1‑19 (MQ=255)

                                                                                   tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagt         >  1:10967052‑M1/1‑18 (MQ=255)

                                                                                    gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtc        <  2:9707006‑M1/96‑80 (MQ=255)

                                                                                    gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtc        >  2:1556414‑M1/1‑17 (MQ=255)

                                                                                    gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagac        <  1:11654640‑M1/96‑80 (MQ=255)

                                                                                     atgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtag       >  2:13008035‑M1/1‑16 (MQ=255)

                                                                                      tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgc      >  2:12184833‑M1/1‑15 (MQ=255)

                                                                                      tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgc      <  1:7917481‑M1/96‑82 (MQ=255)

                                                                                       gatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgcc     <  2:3189990‑M1/96‑83 (MQ=255)

                                                                                        atTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgcca    <  2:1755474‑M1/96‑84 (MQ=255)

                                                                                         ttACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgccaa   <  2:801209‑M1/96‑85 (MQ=255)

                                                                                          tACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgccaat  >  2:1274905‑M1/1‑11 (MQ=255)

TACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGCATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAATACC  >  NC_000913/3754217‑3754402

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