BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,953,377	T→C	T217A (ACG→GCG)	recB ←	RecB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,953,377	0	T	C	99.7%	1039.1 / ‑5.1	321	T217A (ACG→GCG)	recB	RecB
Reads supporting (aligned to +/- strand): ref base T (0/1); new base C (142/178); total (142/179)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 340 total aligned reads displayed.

TCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCC  >  NC_000913/2953287‑2953466

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         >  2:3092086/1‑93 (MQ=255)

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         <  1:6844077/93‑1 (MQ=255)

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         >  2:9143067/1‑93 (MQ=255)

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         >  1:5004858/1‑93 (MQ=255)

tCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGctt                                                                                         >  2:3792513/1‑93 (MQ=255)

 cAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtc                                                                                        <  2:366833/93‑1 (MQ=255)

 cAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtc                                                                                        <  2:1721046/93‑1 (MQ=255)

   gcgcATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtcat                                                                                      <  2:4436316/93‑1 (MQ=255)

     gcATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTtcatca                                                                                    >  1:8521315/1‑93 (MQ=255)

       aTCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCatcatc                                                                                  <  2:5335604/93‑1 (MQ=255)

       aTCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCatcatc                                                                                  >  2:1452621/1‑93 (MQ=255)

        tCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCg                                                                                 >  1:7674122/1‑93 (MQ=255)

        tCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCg                                                                                 >  2:454252/1‑93 (MQ=255)

         cgaGTTCACCCCCTGCGTGGCGCCAGTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCgg                                                                                <  2:7007653/91‑1 (MQ=255)

         ccAGTTCACCCGCTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCgg                                                                                <  1:9416277/93‑1 (MQ=255)

         ccAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCgg                                                                                >  1:2859777/1‑93 (MQ=255)

           aGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggc                                                                              >  1:1732103/1‑93 (MQ=255)

             ttCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggcgg                                                                            >  2:3271759/1‑93 (MQ=255)

              tCACCCACTGCGTCGCGTCACTGCTGTTTTGCCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCACCGCTTCATCATCGggcggc                                                                           >  1:7522167/1‑93 (MQ=255)

              tCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGggcggc                                                                           <  1:3853673/93‑1 (MQ=255)

               cACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGgcggcg                                                                          <  2:2532337/93‑1 (MQ=255)

                 cccACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGt                                                                        <  1:1391897/93‑1 (MQ=255)

                     cTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAACCAGCGCTTCATCATCGGGCGGCGGTGCtt                                                                    >  1:9541074/1‑93 (MQ=255)

                      tGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCttt                                                                   >  1:479816/1‑93 (MQ=255)

                      tGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCttt                                                                   >  1:9021836/1‑93 (MQ=255)

                       gCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTg                                                                  <  1:9636716/93‑1 (MQ=255)

                          tCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATa                                                               <  1:4571354/93‑1 (MQ=255)

                          tCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATa                                                               >  2:9123440/1‑93 (MQ=255)

                            gcgcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAAc                                                             <  2:5088981/93‑1 (MQ=255)

                             cgcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAAcc                                                            >  2:6109796/1‑93 (MQ=255)

                              gcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCg                                                           <  2:6397156/93‑1 (MQ=255)

                              gcCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCg                                                           >  2:3925960/1‑93 (MQ=255)

                               ccACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCgg                                                          >  2:6705716/1‑93 (MQ=255)

                               ccACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCgg                                                          <  2:6232910/93‑1 (MQ=255)

                                 aCTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcg                                                        <  1:1568731/93‑1 (MQ=255)

                                  ctgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcgc                                                       >  1:1003436/1‑93 (MQ=255)

                                  ctgctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGgcgc                                                       >  2:8026767/1‑93 (MQ=255)

                                    gctgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCtt                                                     >  1:9712488/1‑93 (MQ=255)

                                     ctgTTTTACTGTATCAATACGCGCCACATTTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGGGGCGGTGCTTTGATAACCGGCGCTTc                                                    <  2:8093483/93‑1 (MQ=255)

                                      tgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCg                                                   >  2:9368058/1‑93 (MQ=255)

                                      tgTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCg                                                   >  2:1157972/1‑93 (MQ=255)

                                       gTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGc                                                  >  2:8060362/1‑93 (MQ=255)

                                       gTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGc                                                  >  1:3033630/1‑93 (MQ=255)

                                       gTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGc                                                  >  2:7068779/1‑93 (MQ=255)

                                        ttttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGcc                                                 >  1:9060704/1‑93 (MQ=255)

                                         tttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCt                                                <  2:6795047/93‑1 (MQ=255)

                                         tttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCt                                                <  1:8830906/93‑1 (MQ=255)

                                         tttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCt                                                >  1:7939489/1‑93 (MQ=255)

                                         tttACCGTAACAATCTGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCt                                                <  2:1634104/93‑1 (MQ=255)

                                          ttACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCtt                                               <  2:4918439/93‑1 (MQ=255)

                                            aCCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGc                                             <  1:7395655/93‑1 (MQ=255)

                                            aCCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGc                                             >  2:54309/1‑93 (MQ=255)

                                             ccGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATCACCGGCGCTTCGCCTTGCa                                            >  1:8134269/1‑93 (MQ=255)

                                             ccGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCa                                            <  1:5952543/93‑1 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         <  2:5878622/93‑1 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         <  1:3355766/93‑1 (MQ=255)

                                                tATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGAt                                         >  2:1289996/1‑93 (MQ=255)

                                                  tCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATaa                                       >  2:8334468/1‑93 (MQ=255)

                                                  tCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATaa                                       >  2:750438/1‑93 (MQ=255)

                                                  tCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATaa                                       >  1:1904685/1‑93 (MQ=255)

                                                  tCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATaa                                       <  1:9233464/93‑1 (MQ=255)

                                                    aaTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACg                                     >  1:8159913/1‑93 (MQ=255)

                                                     aTACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGa                                    <  1:6777025/93‑1 (MQ=255)

                                                      tACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAt                                   >  2:9621148/1‑93 (MQ=255)

                                                       aCGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGAtt                                  <  1:4394880/93‑1 (MQ=255)

                                                        cgcgCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTa                                 <  2:2741869/93‑1 (MQ=255)

                                                        cgcgCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTa                                 >  1:1109607/1‑93 (MQ=255)

                                                        cgcgCCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTa                                 <  1:5423193/93‑1 (MQ=255)

                                                         gcgcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTaa                                <  1:4819316/93‑1 (MQ=255)

                                                         gcgcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTaa                                <  1:4847285/93‑1 (MQ=255)

                                                          cgcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAat                               <  1:1883062/93‑1 (MQ=255)

                                                           gcCACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCAGCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAata                              <  1:5119198/93‑1 (MQ=255)

                                                            ccACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAatat                             <  1:2207175/93‑1 (MQ=255)

                                                            ccACAATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAatat                             >  1:7969016/1‑93 (MQ=255)

                                                             cacaATTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATc                            >  1:3626267/1‑93 (MQ=255)

                                                                aaTTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATcgcg                         >  2:6531026/1‑93 (MQ=255)

                                                                aaTTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATcgcg                         <  2:6325094/93‑1 (MQ=255)

                                                                 aTTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCgcgc                        <  1:1935479/93‑1 (MQ=255)

                                                                 aTTTGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCgcgc                        >  1:7156456/1‑93 (MQ=255)

                                                                    tGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagc                     <  1:987726/93‑1 (MQ=255)

                                                                    tGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagc                     <  2:2310755/93‑1 (MQ=255)

                                                                    tGCGCGTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCgcagc                     <  1:8896456/93‑1 (MQ=255)

                                                                      ggggTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCaa                   <  1:5877002/90‑1 (MQ=255)

                                                                      cgcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCaa                   <  2:3873129/93‑1 (MQ=255)

                                                                       gcgTGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAAc                  <  2:2052245/93‑1 (MQ=255)

                                                                          tGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGcc               <  1:2084044/93‑1 (MQ=255)

                                                                          tGACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGcc               >  2:7843733/1‑93 (MQ=255)

                                                                           gACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCt              >  1:6439320/1‑93 (MQ=255)

                                                                           gACGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCt              <  2:7391710/93‑1 (MQ=255)

                                                                             cGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGc            <  1:5709395/93‑1 (MQ=255)

                                                                             cGGGAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGc            <  2:4798635/93‑1 (MQ=255)

                                                                               ggAAGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCgg          >  2:689819/1‑93 (MQ=255)

                                                                                  aGCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGccc       <  2:5770537/93‑1 (MQ=255)

                                                                                   gCCAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCt      >  1:833473/1‑93 (MQ=255)

                                                                                     cAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCttt    >  1:2641128/1‑93 (MQ=255)

                                                                                     cAGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCttt    >  1:7839649/1‑93 (MQ=255)

                                                                                      aGCGCTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGTCCCTTTc   >  1:2398428/1‑93 (MQ=255)

                                                                                       gcgcTTCATCATCGGGCGGGGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTcc  <  2:2923368/93‑1 (MQ=255)

                                                                                       gcgcTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTcc  >  1:2802867/1‑93 (MQ=255)

                                                                                       gcgcTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTcc  <  2:954526/93‑1 (MQ=255)

TCAGCGCATCCAGTTCACCCACTGCGTCGCGCCACTGCTGTTTTACCGTATCAATACGCGCCACAATTTGCGCGTGACGGGAAGCCAGCGTTTCATCATCGGGCGGCGGTGCTTTGATAACCGGCGCTTCGCCTTGCAGATAACGATTAATATCGCGCAGCAACGCCTGCGGCCCTTTCC  >  NC_000913/2953287‑2953466

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