BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,141,015	T→C	Q859R (CAG→CGG)	rne ←	RNase E

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,141,015	0	T	C	99.6%	851.9 / NA	249	Q859R (CAG→CGG)	rne	RNase E
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (146/101); total (148/101)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 264 total aligned reads displayed.

ACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACG  >  NC_000913/1140925‑1141104

acaacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGcgccgc                                                                                         >  2:27300/1‑93 (MQ=255)

 caacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCgccgct                                                                                        <  2:5932096/93‑1 (MQ=255)

  aacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctc                                                                                       <  2:1965742/93‑1 (MQ=255)

  aacTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctc                                                                                       <  1:3565552/93‑1 (MQ=255)

   acTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctct                                                                                      >  2:3785798/1‑93 (MQ=255)

   acTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGctct                                                                                      >  2:6299287/1‑93 (MQ=255)

    cTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCtt                                                                                     <  2:813611/93‑1 (MQ=255)

    cTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCtt                                                                                     <  2:2797604/93‑1 (MQ=255)

      ggCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCa                                                                                   >  2:423105/1‑93 (MQ=255)

        cgcgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAAc                                                                                 <  1:4820327/93‑1 (MQ=255)

         gcgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAAct                                                                                >  2:5413262/1‑92 (MQ=255)

         gcgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAAcc                                                                                >  2:5086496/1‑93 (MQ=255)

         gcgcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAAcc                                                                                >  1:1379148/1‑93 (MQ=255)

           gcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCCACCTg                                                                              >  2:4885213/1‑93 (MQ=255)

           gcTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTg                                                                              <  1:5899925/93‑1 (MQ=255)

            cTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGt                                                                             >  1:3391089/1‑93 (MQ=255)

            cTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGt                                                                             >  2:4087395/1‑93 (MQ=255)

             tAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTa                                                                            >  2:2390723/1‑93 (MQ=255)

              aacaacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTAc                                                                           >  1:2927136/1‑93 (MQ=255)

                caacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACAt                                                                         >  1:6292386/1‑93 (MQ=255)

                caacCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACAt                                                                         <  2:2611490/93‑1 (MQ=255)

                  acCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCt                                                                       >  2:5142969/1‑93 (MQ=255)

                  acCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCt                                                                       >  2:5704962/1‑93 (MQ=255)

                  acCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCt                                                                       <  1:2470949/93‑1 (MQ=255)

                  acCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCCCTCTTCAACCTGTACATCt                                                                       >  2:4002475/1‑93 (MQ=255)

                   ccGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGTCGTTCTTCAACCTGTACAACtt                                                                      >  2:3522961/1‑93 (MQ=255)

                    cGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTg                                                                     <  1:280050/93‑1 (MQ=255)

                     ggTTCGATAGCGGCGGCGACAGGGCCCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGc                                                                    <  2:5581409/93‑1 (MQ=255)

                      gTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCg                                                                   <  1:5681988/93‑1 (MQ=255)

                        tCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACCTCTTGCGGa                                                                 >  2:1657629/1‑93 (MQ=255)

                          gATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACg                                                               >  2:1155334/1‑93 (MQ=255)

                           aTAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGt                                                              >  2:1486196/1‑93 (MQ=255)

                           aTAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGt                                                              <  1:2070647/93‑1 (MQ=255)

                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTa                                                             >  2:6427799/1‑93 (MQ=255)

                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTa                                                             >  1:435486/1‑93 (MQ=255)

                            tAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTa                                                             >  1:4155620/1‑93 (MQ=255)

                             aGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTAc                                                            >  2:1191295/1‑93 (MQ=255)

                              gctgcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACa                                                           >  2:2189411/4‑93 (MQ=255)

                              gcggcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACa                                                           >  1:5949316/1‑93 (MQ=255)

                                ggcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAAt                                                         <  2:2041478/93‑1 (MQ=255)

                                ggcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAAt                                                         >  1:2377379/1‑93 (MQ=255)

                                 gcggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAAtt                                                        >  1:815784/1‑93 (MQ=255)

                                  cggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTg                                                       >  1:1775299/1‑93 (MQ=255)

                                  cggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTg                                                       >  2:4309552/1‑93 (MQ=255)

                                   ggcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTgg                                                      <  2:2258513/93‑1 (MQ=255)

                                    gcgACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGa                                                     <  2:5367608/93‑1 (MQ=255)

                                      gACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATa                                                   >  1:2663454/1‑93 (MQ=255)

                                        cAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGc                                                 <  1:4659120/93‑1 (MQ=255)

                                        cAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGc                                                 >  2:2106031/1‑93 (MQ=255)

                                        cAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGc                                                 <  2:2877744/93‑1 (MQ=255)

                                          gggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCgg                                               >  2:6190058/1‑93 (MQ=255)

                                          gggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCgg                                               >  2:6029835/1‑93 (MQ=255)

                                          gggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCgg                                               >  1:4227444/1‑93 (MQ=255)

                                           ggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGa                                              >  2:5546878/1‑93 (MQ=255)

                                           ggACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGa                                              <  1:154942/93‑1 (MQ=255)

                                            gACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGAt                                             >  1:4521683/1‑93 (MQ=255)

                                             agcTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATc                                            >  1:4055532/3‑93 (MQ=255)

                                              ccTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATcc                                           <  2:3891188/93‑1 (MQ=255)

                                                  aGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGAc                                       >  1:4854781/1‑93 (MQ=255)

                                                  aGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGAc                                       <  2:584739/93‑1 (MQ=255)

                                                  aGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGAc                                       <  2:1354265/93‑1 (MQ=255)

                                                   gTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACt                                      >  1:6346664/1‑93 (MQ=255)

                                                    tCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACtt                                     <  1:2562487/93‑1 (MQ=255)

                                                    tCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACtt                                     <  2:4171651/93‑1 (MQ=255)

                                                    tCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACtt                                     <  1:3729008/93‑1 (MQ=255)

                                                    tCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACtt                                     <  2:4762393/93‑1 (MQ=255)

                                                     cACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACttt                                    <  2:3144672/93‑1 (MQ=255)

                                                       ccATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGc                                  >  1:5274464/1‑93 (MQ=255)

                                                        cATCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGcc                                 >  2:1475418/1‑93 (MQ=255)

                                                         aTCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCa                                <  1:2589444/93‑1 (MQ=255)

                                                         aTCGGCTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCa                                <  1:6452480/93‑1 (MQ=255)

                                                            ggCTGCACATGTACTTCTTCCTGGTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCagag                             >  2:1116109/1‑93 (MQ=255)

                                                              cTGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGc                           >  1:2501622/1‑93 (MQ=255)

                                                               tGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGcc                          >  2:2296327/1‑93 (MQ=255)

                                                               tGCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGcc                          >  1:859994/1‑93 (MQ=255)

                                                                gCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCa                         >  2:3358734/1‑93 (MQ=255)

                                                                gCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCa                         <  1:6193643/93‑1 (MQ=255)

                                                                gCACATGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCa                         <  2:1583300/93‑1 (MQ=255)

                                                                 cacaTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAg                        >  2:4954954/1‑93 (MQ=255)

                                                                   caTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGtt                      >  1:1129996/1‑93 (MQ=255)

                                                                   caTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGtt                      >  1:3165003/1‑93 (MQ=255)

                                                                   caTGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGtt                      <  1:5042759/93‑1 (MQ=255)

                                                                     tGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTcc                    <  2:4134204/93‑1 (MQ=255)

                                                                     tGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTcc                    <  1:4494214/93‑1 (MQ=255)

                                                                     tGTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTcc                    <  1:6096479/93‑1 (MQ=255)

                                                                      gTACTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCg                   <  2:2962811/93‑1 (MQ=255)

                                                                        aCTTCTTCCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGga                 >  2:319880/1‑93 (MQ=255)

                                                                          ttcttcCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGgaga               <  2:4434823/93‑1 (MQ=255)

                                                                           tcttcCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCGGAGGCCAGTTCCGGAGAc              <  2:907703/93‑1 (MQ=255)

                                                                           tcttcCTGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGAc              >  2:1008333/1‑93 (MQ=255)

                                                                                 tGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCAc        >  2:6073910/1‑93 (MQ=255)

                                                                                 tGTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCAc        >  1:1302907/1‑93 (MQ=255)

                                                                                  gTTCGCGCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACg       >  2:1893055/1‑93 (MQ=255)

                                                                                     cgcgCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTa    >  2:2112703/1‑93 (MQ=255)

                                                                                     cgcgCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTa    <  1:3483096/93‑1 (MQ=255)

                                                                                     cgcgCCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTa    <  1:1025022/93‑1 (MQ=255)

                                                                                      gcgcCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGCGACGCGCACGCTAc   >  1:3866886/1‑93 (MQ=255)

                                                                                      gcgcCGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTAc   >  2:3224589/1‑93 (MQ=255)

                                                                                       cgccgcTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACg  >  2:5312262/1‑93 (MQ=255)

                                                                                       cgccgcTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACg  <  1:669078/93‑1 (MQ=255)

ACAACTGGCGCGCTAACAACCGGTTCGATAGCGGCGGCGACAGGGACCTCAGTCACCATCGGCTGCACATGTACTTCTTCCTGTTCGCGCTGCTCTTCAACCTGTACATCTTGCGGACGTACAATTGGATAGCGGATCCAGACTTTGCCAGAGGCCAGTTCCGGAGACGCGCACGCTACG  >  NC_000913/1140925‑1141104

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: `‑`