BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,140,817	0	T	C	88.1%	643.9 / 70.4	253	Q925R (CAG→CGG)	rne	RNase E
Reads supporting (aligned to +/- strand): ref base T (12/18); new base C (123/100); total (135/118)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 253 total aligned reads displayed.

GCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGAC  >  NC_000913/1140727‑1140908

 ctgTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGtt                                                                                          >  1:4082995/3‑93 (MQ=255)

  cGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTc                                                                                         <  2:5852421/93‑1 (MQ=255)

  cGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTc                                                                                         <  2:3164358/93‑1 (MQ=255)

  cGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTc                                                                                         >  2:2795645/1‑93 (MQ=255)

  cGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTc                                                                                         >  2:5028684/1‑93 (MQ=255)

  cGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTc                                                                                         <  2:5099181/93‑1 (MQ=255)

   gTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCa                                                                                        >  2:5325928/1‑93 (MQ=255)

    tctcTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAg                                                                                       >  1:1476526/1‑93 (MQ=255)

      tctTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGtt                                                                                     >  2:2586254/1‑93 (MQ=255)

       ctTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTa                                                                                    <  2:54412/93‑1 (MQ=255)

       ctTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTa                                                                                    <  2:6241726/93‑1 (MQ=255)

        ttCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTAc                                                                                   >  2:2616900/1‑93 (MQ=255)

         tCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTAcc                                                                                  <  2:4048741/93‑1 (MQ=255)

          ccTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACcg                                                                                 >  2:3655907/1‑93 (MQ=255)

           cTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACcgc                                                                                <  1:3321750/93‑1 (MQ=255)

           cTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAAATTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACcgc                                                                                >  2:4522657/1‑93 (MQ=255)

            tGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACcgcg                                                                               <  1:4729578/93‑1 (MQ=255)

             gCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCgg                                                                              <  1:2734117/93‑1 (MQ=255)

               ggTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCa                                                                            <  2:2909567/93‑1 (MQ=255)

               ggTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCa                                                                            >  1:2591758/1‑93 (MQ=255)

                       caccGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCAc                                                                    >  2:1242345/1‑93 (MQ=255)

                       caccGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCAc                                                                    <  2:4003186/93‑1 (MQ=255)

                        accGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACt                                                                   >  1:1643779/1‑93 (MQ=255)

                        accGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACt                                                                   >  1:434151/1‑93 (MQ=255)

                         ccGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACtt                                                                  >  1:5047986/1‑93 (MQ=255)

                           ggTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCa                                                                <  1:2950836/93‑1 (MQ=255)

                             ttCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAgg                                                              >  1:3705718/1‑93 (MQ=255)

                             ttCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAgg                                                              >  1:5128778/1‑93 (MQ=255)

                              tCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGa                                                             <  1:6132727/93‑1 (MQ=255)

                              tCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGa                                                             <  1:5233180/93‑1 (MQ=255)

                              tCTGCTTGTTCTGCAACTACCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGa                                                             <  2:4330062/93‑1 (MQ=255)

                               cTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGAt                                                            >  2:1666934/1‑93 (MQ=255)

                               cTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGAt                                                            <  2:2004080/93‑1 (MQ=255)

                               cTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGAt                                                            <  2:3620363/93‑1 (MQ=255)

                                 gCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGc                                                          <  2:1826958/93‑1 (MQ=255)

                                 gCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGc                                                          <  2:1562789/93‑1 (MQ=255)

                                   ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgt                                                        >  2:3832343/1‑93 (MQ=255)

                                   ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgt                                                        >  2:4574734/1‑93 (MQ=255)

                                   ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgt                                                        >  2:3021751/1‑93 (MQ=255)

                                   ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgt                                                        >  1:6469392/1‑93 (MQ=255)

                                   ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgt                                                        <  2:7256201/93‑1 (MQ=255)

                                   ttGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgt                                                        >  1:4398195/1‑93 (MQ=255)

                                    tGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgtc                                                       >  1:3121645/1‑93 (MQ=255)

                                    tGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgtc                                                       <  1:3401667/93‑1 (MQ=255)

                                      ttCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgtcgt                                                     >  2:7179795/1‑93 (MQ=255)

                                      ttCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGcgtcgt                                                     <  1:1650238/93‑1 (MQ=255)

                                       tCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGtt                                                    <  1:232975/93‑1 (MQ=255)

                                         tGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTc                                                  >  1:3027893/1‑93 (MQ=255)

                                           cAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCaa                                                >  2:2094453/1‑93 (MQ=255)

                                           cAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCaa                                                <  2:6758687/93‑1 (MQ=255)

                                           cAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCaa                                                >  1:1371710/1‑93 (MQ=255)

                                             aCTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAAcc                                              >  1:3882421/1‑93 (MQ=255)

                                                tCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACt                                           >  2:7435925/1‑93 (MQ=255)

                                                 ccTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACtt                                          <  1:2821430/93‑1 (MQ=255)

                                                  cTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTc                                         <  2:419020/93‑1 (MQ=255)

                                                  cTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTc                                         >  1:4990399/1‑93 (MQ=255)

                                                  cTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTc                                         <  2:4855779/93‑1 (MQ=255)

                                                   tGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTcc                                        <  1:1984794/93‑1 (MQ=255)

                                                    gggCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCg                                       <  2:4995620/93‑1 (MQ=255)

                                                    gggCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCg                                       >  2:3251346/1‑93 (MQ=255)

                                                    gggCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCg                                       <  1:2054530/93‑1 (MQ=255)

                                                      gCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGt                                     >  1:1027356/1‑93 (MQ=255)

                                                      gCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGt                                     >  1:5773233/1‑93 (MQ=255)

                                                      gCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGt                                     >  1:6366503/1‑93 (MQ=255)

                                                      gCTACGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGt                                     <  1:1280986/93‑1 (MQ=255)

                                                         aCGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGc                                  >  2:1309710/1‑93 (MQ=255)

                                                         aCGGCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGc                                  >  1:6289777/1‑93 (MQ=255)

                                                           ggCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCgg                                <  1:6627589/93‑1 (MQ=255)

                                                           ggCAACATCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCgg                                >  1:6806576/1‑93 (MQ=255)

                                                            gCAACATCAGACTCCGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGt                               >  2:5181776/1‑93 (MQ=255)

                                                                 aTCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGc                          >  2:6035259/1‑93 (MQ=255)

                                                                 aTCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGc                          <  2:1020207/93‑1 (MQ=255)

                                                                 aTCAGACTCGGTAATCACCTCCGGCCGTTCAGTTACATCGGCAGCGGTCACTTCAGGATGCGTCGTTTCAACCAATTCCGGTTGCGGTTCGGc                          >  1:1848667/1‑93 (MQ=255)

                                                                  tCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCa                         <  1:1530722/93‑1 (MQ=255)

                                                                  tCAGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCa                         >  1:1855612/1‑93 (MQ=255)

                                                                    aGACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAc                       <  2:6228697/93‑1 (MQ=255)

                                                                     gACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAcc                      <  1:116797/93‑1 (MQ=255)

                                                                     gACTCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAAcc                      <  1:4174172/93‑1 (MQ=255)

                                                                        tCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGa                   >  2:2240676/1‑93 (MQ=255)

                                                                        tCGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGa                   <  2:1220877/93‑1 (MQ=255)

                                                                         cGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGaa                  >  1:7234653/1‑93 (MQ=255)

                                                                         cGGTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGaa                  <  1:5548896/93‑1 (MQ=255)

                                                                          ggTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAAc                 <  1:7042660/93‑1 (MQ=255)

                                                                           gTAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACg                >  1:6729207/1‑93 (MQ=255)

                                                                            tAATCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACgg               >  1:5995728/1‑93 (MQ=255)

                                                                             aaTCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg              >  2:4142568/1‑93 (MQ=255)

                                                                             aaTCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg              >  2:2288609/1‑93 (MQ=255)

                                                                             aaTCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg              <  1:5434461/93‑1 (MQ=255)

                                                                             aaTCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACggg              <  1:4742860/93‑1 (MQ=255)

                                                                              aTCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACgggg             <  1:6099946/93‑1 (MQ=255)

                                                                              aTCACCTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACgggg             >  1:97052/1‑93 (MQ=255)

                                                                                  ccTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCtt         >  2:1428675/1‑93 (MQ=255)

                                                                                   cTGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTc        >  2:894197/1‑93 (MQ=255)

                                                                                    tGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCt       <  1:70069/93‑1 (MQ=255)

                                                                                    tGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCt       <  1:753462/93‑1 (MQ=255)

                                                                                    tGCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCt       >  1:7284828/1‑93 (MQ=255)

                                                                                     gCGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTa      <  1:925381/93‑1 (MQ=255)

                                                                                      cGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGGACGGGGGCTTCTAc     >  2:4380713/1‑93 (MQ=255)

                                                                                      cGGCCGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTAc     <  1:5669122/93‑1 (MQ=255)

                                                                                         cTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGac  <  2:3450933/93‑1 (MQ=255)

GCCGTCTCTTCCTGCGGTTCAACCACCGGTTCTGCTTGTTCTGCAACTTCCTGGGCTACGGCAACATCAGACTCGGTAATCACCTGCGGCTGTTCAGTTACCGCGGCAGCGATCACTTCAGGATGCGTCGTTTCAACCACTTCCGGTTGCGGTTCGGCAACCTGAACGGGGGCTTCTACGAC  >  NC_000913/1140727‑1140908

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