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.8%	1455.4 / NA	409	R382M (AGG→ATG)	yiaY	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (223/185); minor base T (0/1); total (223/186)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 419 total aligned reads displayed.

GACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGCATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAAT  >  NC_000913/3754220‑3754399

    gCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                                                                                     >  1:5628181/1‑93 (MQ=255)

     cTATCACGAATTTACGGGCAAGGGTTGCAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                    <  2:4196788/93‑1 (MQ=255)

     cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                    >  1:892600/1‑93 (MQ=255)

     cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                    >  1:7598193/1‑93 (MQ=255)

     cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                    <  2:286526/93‑1 (MQ=255)

     cTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGc                                                                                    >  2:6997814/1‑93 (MQ=255)

      tATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCg                                                                                   <  2:1888706‑M1/93‑2 (MQ=255)

      tATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCg                                                                                   >  1:6970902‑M1/1‑92 (MQ=255)

      tATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCg                                                                                   >  1:2697155‑M1/1‑92 (MQ=255)

       aTCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgc                                                                                  >  1:8748006‑M1/1‑91 (MQ=255)

       aTCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgc                                                                                  >  1:1836483‑M1/1‑91 (MQ=255)

         cACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcga                                                                                <  2:6509022‑M1/93‑5 (MQ=255)

         cACGAATGTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcga                                                                                <  2:1296286‑M1/93‑5 (MQ=255)

          aCGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgat                                                                               >  1:1193664‑M1/1‑88 (MQ=255)

           cGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgata                                                                              <  2:7543700‑M1/93‑7 (MQ=255)

           cGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgata                                                                              <  2:171623‑M1/93‑7 (MQ=255)

              aTTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaat                                                                           <  1:3365175‑M1/93‑10 (MQ=255)

               tttAGGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                          <  1:1462705‑M1/93‑11 (MQ=255)

               tttACGGGCAAGGGTTACAGAGATAGTTAAAAAAGCGGCTAACAATTGGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatc                                                                          <  2:1552313‑M1/93‑11 (MQ=255)

                 tACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgc                                                                        <  2:7970622‑M1/93‑13 (MQ=255)

                   cGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgcca                                                                      <  1:1914192‑M1/93‑15 (MQ=255)

                    gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccac                                                                     <  2:2305161‑M1/93‑16 (MQ=255)

                    gggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccac                                                                     >  2:83473‑M1/1‑78 (MQ=255)

                     ggCAAGGGTTACATAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccaca                                                                    <  2:4456273‑M1/93‑17 (MQ=255)

                     ggCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccaca                                                                    >  1:5342752‑M1/1‑77 (MQ=255)

                      gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaa                                                                   >  1:3993822‑M1/1‑76 (MQ=255)

                      gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaa                                                                   >  2:6722276‑M1/1‑76 (MQ=255)

                      gCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaa                                                                   <  1:3732105‑M1/93‑18 (MQ=255)

                        aaGGGTTATAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatt                                                                 >  2:6298169‑M1/1‑74 (MQ=255)

                         aGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaattt                                                                >  1:5307164‑M1/1‑73 (MQ=255)

                          gggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttc                                                               >  2:5824186‑M1/1‑72 (MQ=255)

                          gggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttc                                                               >  1:4370313‑M1/1‑72 (MQ=255)

                           ggTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttct                                                              <  1:7761999‑M1/93‑23 (MQ=255)

                            gTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttctt                                                             <  2:8462513‑M1/93‑24 (MQ=255)

                            gTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttctt                                                             <  2:4266829‑M1/93‑24 (MQ=255)

                               aCAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgt                                                          <  1:3323230‑M1/93‑27 (MQ=255)

                               aCAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgt                                                          >  2:5231311‑M1/1‑67 (MQ=255)

                                cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtg                                                         >  2:534175‑M1/1‑66 (MQ=255)

                                cAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtg                                                         >  1:3646230‑M1/1‑66 (MQ=255)

                                 agagATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtga                                                        >  1:7679090‑M1/1‑65 (MQ=255)

                                  gagaTAGTTTAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgag                                                       <  1:2440421‑M1/93‑30 (MQ=255)

                                  gagaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgag                                                       >  1:2944491‑M1/1‑64 (MQ=255)

                                  gagaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGATTCgcgatcaatcgccacacactcttcgtgag                                                       >  1:5031874‑M1/1‑64 (MQ=255)

                                  aagatAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgag                                                       <  1:4785233‑M1/92‑30 (MQ=255)

                                   agaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagt                                                      >  2:207427‑M1/1‑63 (MQ=255)

                                   agaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagt                                                      <  2:5458686‑M1/93‑31 (MQ=255)

                                    gaTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagtt                                                     <  1:6195543‑M1/93‑32 (MQ=255)

                                     aTAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttg                                                    <  2:6046937‑M1/93‑33 (MQ=255)

                                       aGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcc                                                  >  2:6819371‑M1/1‑59 (MQ=255)

                                       aGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcc                                                  >  1:1316939‑M1/1‑59 (MQ=255)

                                        gTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcct                                                 <  1:8268453‑M1/93‑36 (MQ=255)

                                         ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                <  1:6016819‑M1/93‑37 (MQ=255)

                                         ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                <  2:5590085‑M1/93‑37 (MQ=255)

                                         ttGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                >  1:494504‑M1/1‑57 (MQ=255)

                                         ttGAAAAAGCGGCTAAAAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                <  1:2114455‑M1/93‑37 (MQ=255)

                                         ttGAAAAAACAGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctg                                                <  1:5222374‑M1/93‑37 (MQ=255)

                                           gAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgga                                              >  2:3217574‑M1/1‑55 (MQ=255)

                                            aaagaGCGGCTCACATTTTGCCGGCCGTTGTGGAAATGATGATTACATCACTGCgcgatatatcgccacaatttcttcgtgagttgcctggct                                             <  1:7838080‑M1/93‑40 (MQ=255)

                                             aaaaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatc                                            <  2:4746954‑M1/93‑41 (MQ=255)

                                             aaaaGCGGCTAACAATTTGCCAGCCGCTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatc                                            <  2:3516139‑M1/93‑41 (MQ=255)

                                              gaagCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcg                                           >  1:7443919‑M1/2‑52 (MQ=255)

                                               aaGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgg                                          >  1:7288837‑M1/1‑51 (MQ=255)

                                                    gCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttag                                     <  1:6432818‑M1/93‑48 (MQ=255)

                                                    gCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttag                                     >  1:5920776‑M1/1‑46 (MQ=255)

                                                     cTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttaga                                    >  2:7263182‑M1/1‑45 (MQ=255)

                                                       aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctgtatcggtttagtaa                                  >  1:4626324‑M1/1‑43 (MQ=255)

                                                       aaCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaa                                  >  1:633781‑M1/1‑43 (MQ=255)

                                                        aCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaa                                 <  2:865275‑M1/93‑52 (MQ=255)

                                                        aCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaa                                 >  2:7385849‑M1/1‑42 (MQ=255)

                                                        aCAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaa                                 <  1:1072126‑M1/93‑52 (MQ=255)

                                                         cAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaag                                >  2:3068897‑M1/1‑41 (MQ=255)

                                                         cAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaag                                >  2:7602851‑M1/1‑41 (MQ=255)

                                                         cAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaag                                <  2:3891492‑M1/93‑53 (MQ=255)

                                                          aaTTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagc                               <  2:4056681‑M1/93‑54 (MQ=255)

                                                             ttGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccac                            >  1:8045171‑M1/1‑37 (MQ=255)

                                                              tGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccaca                           <  1:698955‑M1/93‑58 (MQ=255)

                                                               gCCAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacag                          >  1:3214364‑M1/1‑35 (MQ=255)

                                                                ccAGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacagg                         <  1:6600211‑M1/93‑60 (MQ=255)

                                                                  aGCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca                       >  2:4934205‑M1/1‑32 (MQ=255)

                                                                  aGCCGTGGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggca                       <  2:2652538‑M1/93‑62 (MQ=255)

                                                                   gCCGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcat                      >  1:357459‑M1/1‑31 (MQ=255)

                                                                   gCCGTTGTGGAAATGATGATTACATCGCTGCacgataaatcgccacaatttcttcgtgagttgcctggatcgggttaggaaagccacaggcac                      >  1:2794316‑M1/1‑31 (MQ=255)

                                                                    ccGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc                     <  2:8221392‑M1/93‑64 (MQ=255)

                                                                    ccGTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatc                     <  1:3884462‑M1/93‑64 (MQ=255)

                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                   <  1:1491946‑M1/93‑66 (MQ=255)

                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                   >  2:5072198‑M1/1‑28 (MQ=255)

                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                   <  1:6103094‑M1/93‑66 (MQ=255)

                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                   >  2:3811695‑M1/1‑28 (MQ=255)

                                                                      gTTGTGGAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctt                   >  2:2542606‑M1/1‑28 (MQ=255)

                                                                           ggAAATGATGATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagg              <  2:3702360‑M1/93‑71 (MQ=255)

                                                                               atgatgATTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcat          >  1:4158443‑M1/1‑19 (MQ=255)

                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt         >  1:172765‑M1/1‑18 (MQ=255)

                                                                                tgatgaTTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatt         >  2:8739607‑M1/1‑18 (MQ=255)

                                                                                 gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatta        >  2:7128395‑M1/1‑17 (MQ=255)

                                                                                 gatgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcatta        <  1:6732855‑M1/93‑77 (MQ=255)

                                                                                   tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagt      <  2:6383875‑M1/93‑79 (MQ=255)

                                                                                   tgatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagt      <  1:133369‑M1/93‑79 (MQ=255)

                                                                                    gatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtc     >  2:5915564‑M1/1‑14 (MQ=255)

                                                                                    gatTACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtc     >  2:3386147‑M1/1‑14 (MQ=255)

                                                                                       tACATCGCTGCgcgataaatcgccacaatttcttcgtgagttgcctggatcgggttagtaaagccacaggcatctttcagggcattagtcgcc  >  1:3116326‑M1/1‑11 (MQ=255)

GACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACCTCGCTGCATAAATCGCCACAATTTCTTCGTGAGTTGCCTGGATCGGGTTAGTAAAGCCACAGGCATCTTTCAGGGCATTAGTCGCCAAT  >  NC_000913/3754220‑3754399

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