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%	832.3 / ‑5.0	262	I185T (ATT→ACT)	yihG	putative endonuclease
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (114/147); total (115/147)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

Only 100 of 275 total aligned reads displayed.

GCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043539‑4043718

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

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

    ccTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtgg                                                                                     >  1:4634996/1‑93 (MQ=255)

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

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

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

     cTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGCTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGtggtggt                                                                                    <  1:1487338/93‑1 (MQ=255)

      tttGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAgtggtggtg                                                                                   >  1:4564476/1‑93 (MQ=255)

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

          gtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGAt                                                                               <  2:545321/93‑1 (MQ=255)

          gtggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGAt                                                                               >  2:1719419/1‑93 (MQ=255)

           tggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATg                                                                              >  1:2236884/1‑93 (MQ=255)

           tggCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATg                                                                              >  1:2822124/1‑93 (MQ=255)

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

              caacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAg                                                                           >  2:3341056/1‑93 (MQ=255)

              caacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAg                                                                           <  1:697459/93‑1 (MQ=255)

              caacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAg                                                                           <  1:834574/93‑1 (MQ=255)

               aacaaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAgg                                                                          >  2:1834252/1‑93 (MQ=255)

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

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

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

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

                  aaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGa                                                                       >  2:4636919/1‑93 (MQ=255)

                  aaGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGa                                                                       >  1:3170883/1‑93 (MQ=255)

                     ttttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAAc                                                                    >  2:792044/1‑93 (MQ=255)

                     ttttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAAc                                                                    <  2:3538208/93‑1 (MQ=255)

                      tttGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACt                                                                   >  2:666661/1‑93 (MQ=255)

                        tGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACttt                                                                 <  1:682725/93‑1 (MQ=255)

                         gAAAAGAGGAGTGGGGTTGCTGATGTTTTTCTTGCGTGAAACGGGGGCGTTCAACGGAATTGACAGAGGTGGTGGGACGCACGCGAGACtgtt                                                                >  2:4000066/1‑93 (MQ=255)

                           aaaGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTca                                                              >  2:233252/1‑93 (MQ=255)

                           aaaGTGGAGGGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCGTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTca                                                              >  2:4558450/1‑93 (MQ=255)

                             agtggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTcaca                                                            <  1:1025066/93‑1 (MQ=255)

                              gtggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACaa                                                           >  2:2495941/1‑93 (MQ=255)

                               tggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAg                                                          >  1:2103488/1‑93 (MQ=255)

                                ggagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGa                                                         <  1:654389/93‑1 (MQ=255)

                                 gagtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGaa                                                        >  2:3596982/1‑93 (MQ=255)

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

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

                                  agtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAAc                                                       <  1:3267598/93‑1 (MQ=255)

                                   gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      >  2:3216806/1‑93 (MQ=255)

                                   gtggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACg                                                      >  1:3077619/1‑93 (MQ=255)

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

                                    tggGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGc                                                     >  2:3169297/1‑93 (MQ=255)

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

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

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

                                          tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               <  1:3881644/93‑1 (MQ=255)

                                          tGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTg                                               <  2:27853/93‑1 (MQ=255)

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

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

                                            cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGGAGGCGAAACTTTTCACAAGAACACCGAggcgt                                             >  2:3946293/1‑89 (MQ=255)

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

                                            cTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGt                                             <  2:3731631/93‑1 (MQ=255)

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

                                             tgatgaTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGtt                                            >  2:3054909/1‑93 (MQ=255)

                                              gATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGTGTGGttt                                           >  1:569255/1‑93 (MQ=255)

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

                                                tGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCa                                         >  1:1080676/1‑93 (MQ=255)

                                                tGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCa                                         <  2:243378/93‑1 (MQ=255)

                                                tGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCa                                         >  1:1368242/1‑93 (MQ=255)

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

                                                  tttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAAc                                       >  2:3904210/1‑93 (MQ=255)

                                                  tttttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAAc                                       >  1:2382848/1‑93 (MQ=255)

                                                  ttgttCGTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAAc                                       <  2:1359228/90‑1 (MQ=255)

                                                     ttCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATc                                    <  1:293645/93‑1 (MQ=255)

                                                      tcgtGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGAAGGCGAAACTTTTCACAAGAACGCCGAGTGGGTTCAACATCt                                   >  2:3095717/4‑93 (MQ=255)

                                                      tCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCt                                   >  1:3349853/1‑93 (MQ=255)

                                                        ttGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACACGAACGCCGAGTGGTTTCAACATCttt                                 >  1:416059/1‑93 (MQ=255)

                                                        ttGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCttt                                 <  2:2905447/93‑1 (MQ=255)

                                                         tGCGTGAAACGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGTCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTa                                <  1:3302769/93‑1 (MQ=255)

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

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

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

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

                                                                aaCGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTc                         <  2:4329313/93‑1 (MQ=255)

                                                                  cGGGAGCCTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGc                       <  1:2167041/93‑1 (MQ=255)

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

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

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

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

                                                                        ccTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGa                 >  2:612918/1‑93 (MQ=255)

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

                                                                        ccTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAAAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGa                 <  1:3632564/93‑1 (MQ=255)

                                                                         cTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGAt                <  1:3969734/93‑1 (MQ=255)

                                                                         cTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGAt                <  1:2014682/93‑1 (MQ=255)

                                                                         cTTCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGAt                >  1:112875/1‑93 (MQ=255)

                                                                          ttCAACGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATg               <  1:1888068/93‑1 (MQ=255)

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

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

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

                                                                            aaaCGAAATCGACAGTGGTGGTGGTATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTAAACATCTTTACCACGTCGCTCCGGATgtc             <  2:844205/92‑3 (MQ=255)

                                                                             aaCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACg            >  1:4456340/1‑93 (MQ=255)

                                                                              aCGAAATTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCATGTCGCTCCGGATGACGt           >  2:1275233/1‑93 (MQ=255)

                                                                                 aaaTTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGtaat        >  2:384965/1‑93 (MQ=255)

                                                                                   aTTGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTGAACATCTTTACCACGTCGCTCCGGATGACGtaataa      <  1:2602132/93‑1 (MQ=255)

                                                                                    ttGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACACCTTTACCACGCCGCTCCGGATGACGTaataaa     <  2:2564608/93‑1 (MQ=255)

                                                                                     tGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaat    <  2:3291062/93‑1 (MQ=255)

                                                                                     tGACAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTaataaat    <  1:1565062/93‑1 (MQ=255)

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

                                                                                       aCAGTGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAataaataa  <  1:2278854/93‑1 (MQ=255)

GCCGCCTTTGGTGGCAACAAGTTTTGAAAAGTGGAGTGGGTTTGCTGATGTTTTTCTTGCGTGAAACGGGAGCCTTCAACGAAATTGACAATGGTGGTGGGATGCAGGCGAAACTTTTCACAAGAACGCCGAGTGGTTTCAACATCTTTACCACGTCGCTCCGGATGACGTAATAAATAA  >  NC_000913/4043539‑4043718

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