BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,029,060	T→C	D427G (GAC→GGC)	dcm ←	DNA cytosine methylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,029,060	0	T	C	82.6%	431.6 / 91.9	208	D427G (GAC→GGC)	dcm	DNA cytosine methylase
Reads supporting (aligned to +/- strand): ref base T (13/23); new base C (84/87); total (98/110)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.28e-01

Only 100 of 213 total aligned reads displayed.

CGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGA  >  NC_000913/2028969‑2029149

     gTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGaaa                                                                                     >  1:6107355/1‑93 (MQ=255)

     gTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGaaa                                                                                     >  1:2165479/1‑93 (MQ=255)

     gTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGaaa                                                                                     <  1:3622824/93‑1 (MQ=255)

      tttGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAAc                                                                                    >  1:1315591/1‑93 (MQ=255)

      tttGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAAc                                                                                    <  2:5795340/93‑1 (MQ=255)

      tttGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAAc                                                                                    >  1:7419646/1‑93 (MQ=255)

      tttGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAAc                                                                                    >  1:288480/1‑93 (MQ=255)

       ttGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAAcc                                                                                   <  1:982920/93‑1 (MQ=255)

        tGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCg                                                                                  >  2:2684689/1‑93 (MQ=255)

         gATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCgg                                                                                 >  1:4162189/1‑93 (MQ=255)

          aTTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGa                                                                                <  2:2821588/93‑1 (MQ=255)

           tttttGGCTCAAGCATTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAGCCGGaa                                                                               <  1:493892/93‑1 (MQ=255)

           tttttGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGaa                                                                               <  1:3154359/93‑1 (MQ=255)

            ttttGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAAt                                                                              <  1:4548000/93‑1 (MQ=255)

            ttttGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGCGCCCGAAACCGGAAt                                                                              <  1:2369507/93‑1 (MQ=255)

              ttGGCTCAAGCAGTTTTGCCACCTCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATAc                                                                            <  2:1450523/93‑1 (MQ=255)

              ttGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATAc                                                                            <  1:2725188/93‑1 (MQ=255)

               tGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACg                                                                           <  1:1534470/93‑1 (MQ=255)

               tGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACg                                                                           >  1:6012362/1‑93 (MQ=255)

                ggCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACgg                                                                          <  2:1048887/93‑1 (MQ=255)

                   tCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAAt                                                                       <  1:5071592/93‑1 (MQ=255)

                    cAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAAtt                                                                      <  2:937246/93‑1 (MQ=255)

                     aaGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAAttt                                                                     <  1:384422/93‑1 (MQ=255)

                     aaGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAAttt                                                                     <  1:172076/93‑1 (MQ=255)

                        cAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGc                                                                  <  1:5018865/93‑1 (MQ=255)

                         aGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCt                                                                 <  2:2982439/93‑1 (MQ=255)

                         aGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCt                                                                 <  2:4015610/93‑1 (MQ=255)

                          gTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCtt                                                                >  1:4209231/1‑93 (MQ=255)

                          gTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCtt                                                                >  2:6064075/1‑93 (MQ=255)

                           ttttGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTtc                                                               >  1:7151558/1‑93 (MQ=255)

                           ttttGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTtc                                                               >  2:1194045/1‑93 (MQ=255)

                              tGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTcc                                                            <  2:2529658/93‑1 (MQ=255)

                              tGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTcc                                                            <  2:4986209/93‑1 (MQ=255)

                               gCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTccc                                                           <  1:6346563/93‑1 (MQ=255)

                               gCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTccc                                                           >  2:588624/1‑93 (MQ=255)

                                ccACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCg                                                          <  2:5715795/93‑1 (MQ=255)

                                  aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgc                                                        <  2:7492261/93‑1 (MQ=255)

                                  aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgc                                                        >  1:4825196/1‑93 (MQ=255)

                                  aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgc                                                        <  2:2277058/93‑1 (MQ=255)

                                  aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgc                                                        >  1:250628/1‑93 (MQ=255)

                                  aCCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgc                                                        <  2:6600651/93‑1 (MQ=255)

                                    cgcgGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGgcgc                                                      <  1:4020100/93‑1 (MQ=255)

                                        gCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCa                                                  >  1:7077929/1‑93 (MQ=255)

                                          aaaGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCaaa                                                >  1:1002706/1‑93 (MQ=255)

                                           aaGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCaaaa                                               >  2:6754813/1‑93 (MQ=255)

                                           aaGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCaaaa                                               >  1:7488110/1‑93 (MQ=255)

                                           aaGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCaaaa                                               >  2:3315904/1‑93 (MQ=255)

                                             gACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAAcc                                             >  2:3447068/1‑93 (MQ=255)

                                               ccGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCa                                           >  2:3022395/1‑93 (MQ=255)

                                                cGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCAt                                          <  1:6454490/93‑1 (MQ=255)

                                                  gCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTa                                        <  1:5754418/93‑1 (MQ=255)

                                                  gCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTa                                        <  1:3770736/93‑1 (MQ=255)

                                                  gCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTa                                        <  2:1361584/93‑1 (MQ=255)

                                                  gCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTa                                        >  2:1745431/1‑93 (MQ=255)

                                                   cACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGTTTCTCCCGGCGCTTCAAAACCCATTaa                                       <  2:6174339/93‑1 (MQ=255)

                                                   cACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTaa                                       <  1:6828351/93‑1 (MQ=255)

                                                    aCGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAg                                      >  2:4313896/1‑93 (MQ=255)

                                                     cGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgc                                     <  2:4783038/93‑1 (MQ=255)

                                                     cGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgc                                     <  1:1093384/93‑1 (MQ=255)

                                                       accaccGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgc                                   <  2:4587985/93‑1 (MQ=255)

                                                       accaccGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgc                                   <  2:4196997/93‑1 (MQ=255)

                                                       accaccGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgc                                   >  2:3980614/1‑93 (MQ=255)

                                                         caccGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgcgc                                 >  1:5965167/1‑93 (MQ=255)

                                                         caccGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgcgc                                 >  2:5914711/1‑93 (MQ=255)

                                                          accGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgcgcg                                >  2:1566369/1‑93 (MQ=255)

                                                          accGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgcgcg                                >  2:3773046/1‑93 (MQ=255)

                                                          accGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgcgcg                                >  2:5552814/1‑93 (MQ=255)

                                                          accGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAgcgcgcg                                >  1:5622309/1‑93 (MQ=255)

                                                            cGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCa                              >  2:760231/1‑93 (MQ=255)

                                                            cGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCa                              <  1:1656122/93‑1 (MQ=255)

                                                            cGAGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCa                              <  2:5635964/93‑1 (MQ=255)

                                                              aGTTACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCAtt                            <  2:3261159/93‑1 (MQ=255)

                                                                ttACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTcc                          <  2:5383730/93‑1 (MQ=255)

                                                                ttACCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTcc                          >  2:4545043/1‑93 (MQ=255)

                                                                  aCCGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCg                        >  2:7330271/1‑93 (MQ=255)

                                                                   ccGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGa                       <  1:497704/93‑1 (MQ=255)

                                                                   ccGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGa                       <  1:2065653/93‑1 (MQ=255)

                                                                    cGAACTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAg                      <  1:915099/93‑1 (MQ=255)

                                                                      aaCTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGc                    >  2:1946040/1‑93 (MQ=255)

                                                                        cTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGt                  <  2:2169639/93‑1 (MQ=255)

                                                                        cTGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGt                  <  1:6960432/93‑1 (MQ=255)

                                                                         tGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGtt                 >  1:6629171/1‑93 (MQ=255)

                                                                         tGGCGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGtt                 <  1:2579572/93‑1 (MQ=255)

                                                                            cGATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAAc              <  1:1474453/93‑1 (MQ=255)

                                                                             gATAGGCCTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAAcc             >  2:7164078/1‑93 (MQ=255)

                                                                              aTAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCg            >  2:4832407/1‑93 (MQ=255)

                                                                               tAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGa           <  1:3320826/93‑1 (MQ=255)

                                                                                  gCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGt        >  2:2640744/1‑93 (MQ=255)

                                                                                   ccTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTg       <  2:3215346/93‑1 (MQ=255)

                                                                                    cTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTgg      <  1:7220769/93‑1 (MQ=255)

                                                                                    cTGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTgg      <  2:306577/93‑1 (MQ=255)

                                                                                     tGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa     <  2:5762515/93‑1 (MQ=255)

                                                                                     tGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa     >  1:1426332/1‑93 (MQ=255)

                                                                                     tGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa     >  1:4044256/1‑93 (MQ=255)

                                                                                     tGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa     >  2:3225555/1‑93 (MQ=255)

                                                                                     tGAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGa     <  1:5475313/93‑1 (MQ=255)

                                                                                      gAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGAc    >  2:2868118/1‑93 (MQ=255)

                                                                                      gAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGAc    >  2:37483/1‑93 (MQ=255)

                                                                                      gAGTGCCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGAc    >  2:1882385/1‑93 (MQ=255)

                                                                                        gtgcCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGa  <  1:4257933/93‑1 (MQ=255)

CGCCTGTTTGATTTTTGGCTCAAGCAGTTTTGCCACCGCGGCAAAGACCGGCACGACCACCGAGTTACCGAACTGGCGATAGGCCTGAGTGTCCGAAACCGGAATACGGAATTTCGCTTCTCCCGGCGCTTCAAAACCCATTAAGCGCGCGCATTCCCGAGGCGTTAACCGACGTGGACGA  >  NC_000913/2028969‑2029149

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