BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,506,318	(C)_5→6	coding (1134/1227 nt)	yhfW ←	putative mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,506,313	1	.	C	99.2%	890.5 / ‑1.6	261	coding (1139/1227 nt)	yhfW	putative mutase
Reads supporting (aligned to +/- strand): ref base . (2/0); new base C (109/150); total (111/150)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 287 total aligned reads displayed.

AAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGG  >  NC_000913/3506221‑3506400

aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAg                                                                                          <  1:1211352/93‑1 (MQ=255)

 aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑c                                                                                        >  2:6251287/1‑93 (MQ=255)

      aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cccccc                                                                                   >  1:6179673/1‑92 (MQ=255)

      aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑cccccc                                                                                   >  2:7303222/1‑92 (MQ=255)

        aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑ccccccac                                                                                 <  1:1425957/93‑4 (MQ=255)

         gggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑ccccccaca                                                                                >  2:8087574/1‑89 (MQ=255)

           gtgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATc                                                                               <  1:7555117/93‑1 (MQ=255)

              tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCCCATCAGa                                                                            >  1:8605694/1‑93 (MQ=255)

                   tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCg                                                                       <  1:3785826/93‑1 (MQ=255)

                    gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgt                                                                      <  2:3695303/93‑1 (MQ=255)

                      aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtgg                                                                    <  2:5058280/93‑1 (MQ=255)

                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtggt                                                                   <  1:2243067/93‑1 (MQ=255)

                       aaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtggt                                                                   <  2:3165888/93‑1 (MQ=255)

                        aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtggtg                                                                  <  2:1547886/93‑1 (MQ=255)

                        aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCgtggtg                                                                  <  2:6739264/93‑1 (MQ=255)

                           ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTgcgc                                                               >  2:3442453/1‑93 (MQ=255)

                             agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGcac                                                             >  2:2564934/1‑93 (MQ=255)

                               ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCacac                                                           >  2:4376293/1‑93 (MQ=255)

                               ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACAAGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCacac                                                           >  2:6552380/1‑93 (MQ=255)

                                gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACAcc                                                          <  2:4721572/93‑1 (MQ=255)

                                 aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                         <  2:3940620/93‑1 (MQ=255)

                                 aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                         >  1:1979334/1‑93 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                        >  1:3926051/1‑93 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                        >  2:7912755/1‑93 (MQ=255)

                                  ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                        >  1:1146345/1‑93 (MQ=255)

                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCCCACCGAGc                                                      >  1:7429149/1‑93 (MQ=255)

                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                      >  1:8811124/1‑93 (MQ=255)

                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                      <  1:6455752/93‑1 (MQ=255)

                                    aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                      >  1:3200469/1‑93 (MQ=255)

                                        aaaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCg                                                  >  1:5940170/1‑93 (MQ=255)

                                         aaagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGt                                                 >  2:336090/1‑93 (MQ=255)

                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                >  1:8272214/1‑93 (MQ=255)

                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                <  2:3730475/93‑1 (MQ=255)

                                          aagaGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTa                                                <  2:577844/93‑1 (MQ=255)

                                           agagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                               <  1:5052341/93‑1 (MQ=255)

                                           agagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                               >  1:7627573/1‑93 (MQ=255)

                                            gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                              <  2:5459474/93‑1 (MQ=255)

                                            gagCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGc                                              <  1:4249763/93‑1 (MQ=255)

                                             agCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                             <  1:4510850/93‑1 (MQ=255)

                                              gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                            <  1:3134654/93‑1 (MQ=255)

                                              gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                            >  2:1735083/1‑93 (MQ=255)

                                              gCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGa                                            >  1:981224/1‑93 (MQ=255)

                                               cGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                           >  1:8521035/1‑93 (MQ=255)

                                               cGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                           >  1:6272396/1‑93 (MQ=255)

                                               cGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                           <  2:8489386/93‑1 (MQ=255)

                                               cGACCATTTTGTGGCGGGGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGat                                           <  2:8120723/93‑1 (MQ=255)

                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                          >  2:1446541/1‑93 (MQ=255)

                                                gACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatc                                          >  1:7083101/1‑93 (MQ=255)

                                                 aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                         <  1:7907381/93‑1 (MQ=255)

                                                 aCCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                         >  2:7785975/1‑93 (MQ=255)

                                                  ccATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatcat                                        >  2:2307991/1‑93 (MQ=255)

                                                    aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                      <  1:1513633/93‑1 (MQ=255)

                                                    aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                      <  1:1724417/93‑1 (MQ=255)

                                                    aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                      >  2:8218807/1‑93 (MQ=255)

                                                    aTTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                      >  1:7233128/1‑93 (MQ=255)

                                                     ttttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATccc                                     >  2:86167/1‑93 (MQ=255)

                                                      tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                    >  2:207467/1‑93 (MQ=255)

                                                      tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                    <  1:1797959/93‑1 (MQ=255)

                                                      tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                    <  2:8697941/93‑1 (MQ=255)

                                                      tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                    <  2:5704778/93‑1 (MQ=255)

                                                      tttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                    <  1:6706977/93‑1 (MQ=255)

                                                       ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                   <  2:3327803/93‑1 (MQ=255)

                                                       ttGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                   <  1:4632486/93‑1 (MQ=255)

                                                        tgtgGCGGTGGCGCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTg                                  >  2:1908758/1‑93 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACGTCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                 <  1:5434746/93‑1 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                 <  1:2691816/93‑1 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                 <  2:6602333/93‑1 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                 <  1:5746984/93‑1 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                 <  2:250045/93‑1 (MQ=255)

                                                         gtggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                 <  1:1020048/93‑1 (MQ=255)

                                                          tggcggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                >  2:6063763/1‑93 (MQ=255)

                                                              ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGACAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                            >  1:1874920/1‑93 (MQ=255)

                                                              ggtggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATa                            <  1:6910196/93‑1 (MQ=255)

                                                                 ggcgCGCGGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                         >  2:1607605/1‑93 (MQ=255)

                                                                     cgcgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                     <  1:8262628/93‑1 (MQ=255)

                                                                      gcgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCa                    >  2:4534694/1‑93 (MQ=255)

                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                   <  1:3288819/93‑1 (MQ=255)

                                                                       cgGAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCAc                   <  1:2323844/93‑1 (MQ=255)

                                                                        ggAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACt                  <  2:7646107/93‑1 (MQ=255)

                                                                         gAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                 >  2:3478159/1‑93 (MQ=255)

                                                                         gAAAAATTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                 >  1:6598264/1‑93 (MQ=255)

                                                                           aaaaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc               <  1:4070635/93‑1 (MQ=255)

                                                                           aaaaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGc               >  2:3254559/1‑93 (MQ=255)

                                                                             aaTTCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcac             <  1:801771/93‑1 (MQ=255)

                                                                               ttCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacca           <  2:7184282/93‑1 (MQ=255)

                                                                               tgcacacgCGGTAGCCCCCCTCATCAGAAAGCGTGGTTCGCACACCGAGCTGCGTAGCCATCATCCCGTGCTGATAAACCAGCACTGGcacca           <  1:495506/91‑1 (MQ=255)

                                                                                tCACACACGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcaccac          <  1:2824567/93‑1 (MQ=255)

                                                                                 cacacaCGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt         <  1:4771682/93‑1 (MQ=255)

                                                                                 cacacaCGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt         <  1:5316499/93‑1 (MQ=255)

                                                                                 cacacaCGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACt         <  2:640452/93‑1 (MQ=255)

                                                                                  acacacGGTAGCCCCCCACATGAGAAAGCGTGGTGCGCACGCCGAGCTGCGTAGCGATCAACCCTTGCTGATAAACCAGCACTGGCACCACtt        <  1:8300496/93‑1 (MQ=255)

                                                                                  acacacGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACtt        <  2:2965397/93‑1 (MQ=255)

                                                                                   cacacGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc       >  1:8660597/1‑93 (MQ=255)

                                                                                   cacacGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc       >  2:1987085/1‑93 (MQ=255)

                                                                                    acacGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcg      <  2:2781689/93‑1 (MQ=255)

                                                                                      ccGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTcgcg    <  2:6680480/92‑1 (MQ=255)

                                                                                       cGGTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg   <  1:2586636/93‑1 (MQ=255)

                                                                                        ggTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg  >  1:4729098/1‑93 (MQ=255)

                                                                                        ggTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg  >  2:2240050/1‑93 (MQ=255)

                                                                                        ggTAGCCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg  <  2:5005572/93‑1 (MQ=255)

AAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGCGGAAAAATTCACACACGGTAG‑CCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGG  >  NC_000913/3506221‑3506400

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