BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,043,629	A→G	I185T (ATT→ACT)	yihG ←	putative endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,043,629	0	A	G	99.6%	1822.4 / ‑4.3	557	I185T (ATT→ACT)	yihG	putative endonuclease
Reads supporting (aligned to +/- strand): ref base A (1/0); major base G (239/315); minor base T (0/2); total (240/317)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 582 total aligned reads displayed.

GCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043539‑4043718

gccgccGTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACTGGAGCCTTCAACGAAATTGACAGtg                                                                                         <  2:2593593/93‑1 (MQ=255)

 ccgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtgg                                                                                        >  1:4605263/1‑93 (MQ=255)

 ccgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtgg                                                                                        <  1:618772/93‑1 (MQ=255)

  cgccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggt                                                                                       <  2:6129491/93‑1 (MQ=255)

    ccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtgg                                                                                     <  2:4537133/93‑1 (MQ=255)

     cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggt                                                                                    <  2:6346064/93‑1 (MQ=255)

     cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggt                                                                                    >  2:8217041/1‑93 (MQ=255)

     cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggt                                                                                    <  2:1726059/93‑1 (MQ=255)

       ttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggtgg                                                                                  >  1:3648923/1‑93 (MQ=255)

       ttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggtgg                                                                                  <  1:6708881/93‑1 (MQ=255)

       ttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggtgg                                                                                  <  2:5774392/93‑1 (MQ=255)

        tggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTggg                                                                                 <  1:7153103/93‑1 (MQ=255)

         ggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGa                                                                                <  2:2738230/93‑1 (MQ=255)

         ggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGa                                                                                <  2:2863029/93‑1 (MQ=255)

         ggtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGa                                                                                >  1:1287926/1‑93 (MQ=255)

            ggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGc                                                                             >  2:4053907/1‑93 (MQ=255)

            ggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGc                                                                             <  2:1420638/93‑1 (MQ=255)

             gCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCa                                                                            <  1:3576476/93‑1 (MQ=255)

             gCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCa                                                                            <  1:1342226/93‑1 (MQ=255)

              caacaaGTTGTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAg                                                                           <  1:3747143/93‑1 (MQ=255)

               aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          <  2:1220038/93‑1 (MQ=255)

               aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          <  2:8606038/93‑1 (MQ=255)

                acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         >  1:3358750/1‑93 (MQ=255)

                acaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGc                                                                         <  2:1486208/93‑1 (MQ=255)

                 caaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCg                                                                        >  1:5290084/1‑93 (MQ=255)

                 caaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCg                                                                        >  1:1914514/1‑93 (MQ=255)

                        tGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACttt                                                                 >  1:7611355/1‑93 (MQ=255)

                          aaaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTc                                                               <  2:9233334/93‑1 (MQ=255)

                          aaaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTc                                                               <  1:8628850/93‑1 (MQ=255)

                          aaaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTc                                                               <  1:439897/93‑1 (MQ=255)

                           aaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTca                                                              <  2:8457360/93‑1 (MQ=255)

                                tgagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGa                                                         <  1:2893161/92‑1 (MQ=255)

                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  2:1049208/93‑1 (MQ=255)

                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  1:1170395/93‑1 (MQ=255)

                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        <  1:5506384/93‑1 (MQ=255)

                                  agttgGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCTCAAGAAc                                                       <  2:9854611/93‑1 (MQ=255)

                                  agtggGGTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAAc                                                       <  2:7302197/93‑1 (MQ=255)

                                   gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTAACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      >  2:3578367/1‑93 (MQ=255)

                                    tggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGc                                                     <  1:8807980/93‑1 (MQ=255)

                                     gggTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGcc                                                    <  2:6053995/93‑1 (MQ=255)

                                     gggTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGcc                                                    <  2:372959/93‑1 (MQ=255)

                                      ggTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCg                                                   <  2:6852070/93‑1 (MQ=255)

                                       gTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGa                                                  <  1:412834/93‑1 (MQ=255)

                                        tttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAg                                                 >  2:7634214/1‑93 (MQ=255)

                                        tttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAg                                                 <  1:8519950/93‑1 (MQ=255)

                                        tttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAg                                                 <  2:3426785/93‑1 (MQ=255)

                                         ttGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGt                                                <  2:5553564/93‑1 (MQ=255)

                                          tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               >  2:331051/1‑93 (MQ=255)

                                           gCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTgg                                              >  2:6237585/1‑93 (MQ=255)

                                            cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             >  1:4226524/1‑93 (MQ=255)

                                            cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             >  2:1890773/1‑93 (MQ=255)

                                            cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             >  2:4808259/1‑93 (MQ=255)

                                              gATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGttt                                           <  1:7769491/93‑1 (MQ=255)

                                              gATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGttt                                           >  2:8144771/1‑93 (MQ=255)

                                               aTGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTc                                          <  2:9411316/93‑1 (MQ=255)

                                                 gTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCaa                                        <  2:6757104/93‑1 (MQ=255)

                                                  tttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAAc                                       <  1:2412407/93‑1 (MQ=255)

                                                   ttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACa                                      <  2:7983882/93‑1 (MQ=255)

                                                   ttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACa                                      >  2:3044129/1‑93 (MQ=255)

                                                          gCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAGCGCCGAGTGGTTTCAACATCTTTAc                               <  2:7757866/93‑1 (MQ=255)

                                                          gCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTAc                               <  2:5145777/93‑1 (MQ=255)

                                                            gTGAAACGTGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  1:5261751/93‑1 (MQ=255)

                                                            gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  2:4368863/93‑1 (MQ=255)

                                                            gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  2:564734/93‑1 (MQ=255)

                                                            gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  1:8095108/93‑1 (MQ=255)

                                                            gTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCa                             <  1:5147762/93‑1 (MQ=255)

                                                              gAAACGGGAGCCTTCAACGAAATTGACAGCGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACg                           >  1:6298896/1‑93 (MQ=255)

                                                               aaaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACgc                          >  1:6097458/1‑92 (MQ=255)

                                                               aaaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGt                          <  1:212096/93‑1 (MQ=255)

                                                                aaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTc                         >  1:6842265/1‑93 (MQ=255)

                                                                  cGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGc                       >  1:7427291/1‑93 (MQ=255)

                                                                   gggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCt                      <  1:5042575/93‑1 (MQ=255)

                                                                   gggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCt                      <  1:5840172/93‑1 (MQ=255)

                                                                   gggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCt                      <  1:9175592/93‑1 (MQ=255)

                                                                    ggAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTc                     >  1:10036098/1‑93 (MQ=255)

                                                                     gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATGTTTACCCCGTCGCTcc                    >  2:5346264/1‑93 (MQ=255)

                                                                     gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTcc                    <  2:876749/93‑1 (MQ=255)

                                                                     gAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTcc                    <  2:7171447/93‑1 (MQ=255)

                                                                       gCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCgg                  >  1:9807880/1‑93 (MQ=255)

                                                                       gCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCgg                  <  1:8277448/93‑1 (MQ=255)

                                                                        ccTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGa                 <  2:3920041/93‑1 (MQ=255)

                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  2:9131102/1‑93 (MQ=255)

                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               >  2:4691926/1‑93 (MQ=255)

                                                                           tCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGa              >  1:3934028/1‑93 (MQ=255)

                                                                           tCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGa              <  1:4769757/93‑1 (MQ=255)

                                                                           tCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGa              <  1:2689371/93‑1 (MQ=255)

                                                                            cAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGAc             >  1:403096/1‑93 (MQ=255)

                                                                            cAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGAc             <  2:591199/93‑1 (MQ=255)

                                                                             aaCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACg            <  1:696760/93‑1 (MQ=255)

                                                                              aCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGt           <  1:9954328/93‑1 (MQ=255)

                                                                               cGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGta          >  2:7824126/1‑93 (MQ=255)

                                                                               cGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGta          <  2:968312/93‑1 (MQ=255)

                                                                                gAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaa         >  1:4911697/1‑93 (MQ=255)

                                                                                gAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaa         <  1:961806/93‑1 (MQ=255)

                                                                                 aaaTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGCTTCAACATCTTTACCACGTCGCTCCGGATGACGtcat        >  1:1573220/1‑93 (MQ=255)

                                                                                  aaTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaata       >  2:6556234/1‑93 (MQ=255)

                                                                                    ttGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaa     >  2:222490/1‑93 (MQ=255)

                                                                                    ttGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaa     <  2:1283933/93‑1 (MQ=255)

                                                                                     tGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaat    >  2:4178685/1‑93 (MQ=255)

                                                                                      gACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaata   <  2:8913084/93‑1 (MQ=255)

GCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043539‑4043718

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