BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,035,428	A→G	noncoding (330/2905 nt)	rrlA →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,035,428	0	A	G	100.0%	389.4 / NA	164	noncoding (330/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (76/88); total (76/88)

Only 100 of 174 total aligned reads displayed.

CAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG  >  NC_000913/4035338‑4035506

cAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCa                                                                              >  2:1764940/1‑93 (MQ=18)

cAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCa                                                                              >  1:996769/1‑93 (MQ=18)

cAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCa                                                                              >  1:3485791/1‑93 (MQ=18)

 aGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAg                                                                             >  2:4152942/1‑93 (MQ=18)

  gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGc                                                                            >  2:3287849/1‑93 (MQ=18)

  gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGAGATACAGGGTGGCAGc                                                                            >  2:1355167/1‑93 (MQ=16)

   tAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGcc                                                                           <  1:3483425/93‑1 (MQ=18)

    aGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGccc                                                                          <  2:350315/93‑1 (MQ=18)

      cggcgAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCg                                                                        <  1:1233086/93‑1 (MQ=18)

       ggcgAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGt                                                                       >  1:3983904/1‑93 (MQ=18)

          gagcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTaca                                                                    <  1:3630271/93‑1 (MQ=18)

          gagcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTaca                                                                    <  1:1171904/93‑1 (MQ=18)

          gagcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTaca                                                                    <  2:705106/93‑1 (MQ=21)

          gagcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTaca                                                                    <  2:1740973/93‑1 (MQ=18)

           agcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTacac                                                                   >  1:886073/1‑93 (MQ=18)

           agcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTacac                                                                   <  1:4420243/93‑1 (MQ=21)

            gcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTAcaca                                                                  >  2:3667997/1‑93 (MQ=18)

             cgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaa                                                                 <  2:30718/93‑1 (MQ=21)

              gaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaaa                                                                <  2:582714/93‑1 (MQ=18)

              gaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaaa                                                                <  2:4133061/93‑1 (MQ=18)

                aCGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaaaaa                                                              >  1:649397/1‑93 (MQ=18)

                   gggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGc                                                           >  1:1363197/1‑93 (MQ=18)

                    ggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGca                                                          <  1:4270/93‑1 (MQ=21)

                    ggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGca                                                          <  1:1492122/93‑1 (MQ=21)

                     gAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcac                                                         <  1:223290/93‑1 (MQ=18)

                      agcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcaca                                                        <  2:3940869/93‑1 (MQ=21)

                      agcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcaca                                                        <  2:4678780/93‑1 (MQ=21)

                       gcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACAt                                                       <  1:1915260/93‑1 (MQ=21)

                        cagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATg                                                      >  1:2292054/1‑93 (MQ=18)

                         agcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGc                                                     >  1:1923117/1‑93 (MQ=18)

                          gcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCt                                                    >  2:4044097/1‑93 (MQ=18)

                          gcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCt                                                    >  2:285335/1‑93 (MQ=18)

                           cccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtg                                                   >  2:1610503/1‑93 (MQ=18)

                            ccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtgt                                                  <  1:2758604/93‑1 (MQ=18)

                             cAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtgtg                                                 >  2:3519015/1‑93 (MQ=18)

                              agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                >  2:332489/1‑93 (MQ=18)

                              agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                <  1:1286053/93‑1 (MQ=18)

                               gagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAg                                               >  2:917385/1‑93 (MQ=18)

                                 gCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCt                                             >  1:3435476/1‑93 (MQ=18)

                                 gCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCt                                             <  1:2498267/93‑1 (MQ=18)

                                 gCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCt                                             >  1:2476218/1‑93 (MQ=21)

                                  ccTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTc                                            >  1:3607915/1‑93 (MQ=18)

                                    tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          >  1:3207607/1‑93 (MQ=18)

                                    tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          >  1:4606126/1‑93 (MQ=18)

                                     gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGAt                                         >  1:1858541/1‑93 (MQ=18)

                                      aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        >  1:2157175/1‑93 (MQ=18)

                                       aTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGa                                       >  1:3040648/1‑93 (MQ=18)

                                       aTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGa                                       >  1:2108637/1‑93 (MQ=18)

                                        tCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAg                                      <  1:3724584/93‑1 (MQ=18)

                                        tCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAg                                      <  2:1542245/93‑1 (MQ=18)

                                          aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  1:4658595/93‑1 (MQ=18)

                                          aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  1:3904477/93‑1 (MQ=18)

                                           gtgtgtgtgtTAGTGGACGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   <  2:1284990/93‑1 (MQ=16)

                                           gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   >  1:2796700/1‑93 (MQ=18)

                                           gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   <  1:4561495/93‑1 (MQ=18)

                                           gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   >  1:1490050/1‑93 (MQ=18)

                                           gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   >  1:4142271/1‑93 (MQ=18)

                                           gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGGGAGCTGGATGAGTAg                                   >  2:3916646/1‑93 (MQ=255)

                                            tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  <  1:571191/93‑1 (MQ=18)

                                            tgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgg                                  <  2:3981608/93‑1 (MQ=18)

                                             gtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAggg                                 <  1:2434200/93‑1 (MQ=18)

                                             gtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAggg                                 <  1:2353502/93‑1 (MQ=18)

                                              tgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                                <  1:1020376/93‑1 (MQ=18)

                                              tgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                                >  1:2256391/1‑93 (MQ=18)

                                               gtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCg                               >  1:4493274/1‑93 (MQ=18)

                                                tgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCgg                              >  1:4091077/1‑93 (MQ=18)

                                                tgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCgg                              <  2:234210/93‑1 (MQ=18)

                                                 gtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCggg                             <  1:3615969/93‑1 (MQ=18)

                                                  tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGa                            >  2:3759907/1‑93 (MQ=18)

                                                  tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGGAGGGCGGGa                            >  1:4236656/1‑93 (MQ=16)

                                                   gtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac                           >  1:1124405/1‑93 (MQ=21)

                                                    ttAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGaca                          >  1:2575579/1‑93 (MQ=18)

                                                    ttAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGaca                          <  2:2643381/93‑1 (MQ=18)

                                                      aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCTGGTCACg                        >  1:1038470/1‑93 (MQ=255)

                                                      aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACg                        <  2:3171757/93‑1 (MQ=18)

                                                       gTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGt                       >  2:521898/1‑93 (MQ=18)

                                                            aGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTAt                  <  2:2951746/93‑1 (MQ=18)

                                                             gCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc                 >  2:853970/1‑93 (MQ=18)

                                                             gCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc                 >  2:2222922/1‑93 (MQ=18)

                                                             gCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATc                 <  2:3859824/93‑1 (MQ=18)

                                                               gTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCt               <  2:2601079/93‑1 (MQ=18)

                                                                tCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg              >  2:2013694/1‑93 (MQ=18)

                                                                tCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg              <  2:2402610/93‑1 (MQ=18)

                                                                 cTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGt             <  2:1150376/93‑1 (MQ=18)

                                                                  tGGAAAGGCGTGCGATACAGGGTGGTAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc            >  1:2608131/1‑93 (MQ=255)

                                                                  tGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGGCACGTGGTATCCTGTc            >  2:896041/1‑93 (MQ=16)

                                                                   ggAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCt           >  1:1769313/1‑93 (MQ=18)

                                                                   ggAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCt           >  2:1387837/1‑93 (MQ=18)

                                                                   ggAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCt           <  2:2977377/93‑1 (MQ=18)

                                                                    gAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTg          >  2:3508855/1‑93 (MQ=18)

                                                                     aaaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGa         >  1:4143952/1‑93 (MQ=18)

                                                                       aGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAat       >  1:387768/1‑93 (MQ=18)

                                                                       aGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAat       >  2:1833375/1‑93 (MQ=21)

                                                                        ggCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAata      <  2:4635914/93‑1 (MQ=18)

                                                                        ggCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAata      <  2:854864/93‑1 (MQ=18)

                                                                         gCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAatat     <  2:1917877/93‑1 (MQ=18)

                                                                           gTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgg   <  2:2558832/93‑1 (MQ=17)

                                                                           gTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgg   >  1:3649215/1‑93 (MQ=17)

                                                                            tGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggg  <  1:3055225/93‑1 (MQ=14)

                                                                            tGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggg  <  1:3506585/93‑1 (MQ=14)

CAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG  >  NC_000913/4035338‑4035506

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