BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,485,021	A→G	S46G (AGT→GGT)	b4262 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,485,021	0	A	G	99.1%	1070.9 / ‑2.2	342	S46G (AGT→GGT)	b4262	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (164/175); total (165/177)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 360 total aligned reads displayed.

TTCACCACCATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGAGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTGCTTGGG  >  NC_000913/4484931‑4485110

ttCACCACTATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGt                                                                                         <  2:1671484/93‑1 (MQ=255)

ttCACCACCATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGt                                                                                         <  1:7467927/93‑1 (MQ=255)

  caccacCATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTa                                                                                       >  2:4472165/1‑93 (MQ=255)

   accaccATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTac                                                                                      >  2:4251759/1‑93 (MQ=255)

   accaccATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGAGTTac                                                                                      >  2:2731542/1‑93 (MQ=255)

     caccaTCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTacga                                                                                    <  2:7069161/93‑1 (MQ=255)

       ccaTCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTccgacg                                                                                  >  2:3247204/1‑93 (MQ=255)

       ccaTCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTacgacg                                                                                  >  2:7697648/1‑93 (MQ=255)

       ccaTCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTacgacg                                                                                  >  1:1960817/1‑93 (MQ=255)

         atcatGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGAcgcg                                                                                >  2:216310/1‑93 (MQ=255)

           catGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGtt                                                                              >  1:3813502/1‑93 (MQ=255)

            atgatgACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTa                                                                             <  1:6237427/93‑1 (MQ=255)

            atgatgACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTa                                                                             >  1:2703784/1‑93 (MQ=255)

            atgatgACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTa                                                                             <  1:7134929/93‑1 (MQ=255)

             tgatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAg                                                                            >  2:5469694/1‑93 (MQ=255)

             tgatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAg                                                                            <  2:7272601/93‑1 (MQ=255)

             tgatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAg                                                                            >  2:2796449/1‑93 (MQ=255)

             tgatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAg                                                                            >  2:7588902/1‑93 (MQ=255)

             tgatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAg                                                                            >  1:1253820/1‑93 (MQ=255)

             tgatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAg                                                                            <  2:2453390/93‑1 (MQ=255)

              gatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAgg                                                                           <  1:2496536/93‑1 (MQ=255)

              gatgaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAgg                                                                           <  1:7292995/93‑1 (MQ=255)

                 gaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGgcgc                                                                        <  1:5791408/93‑1 (MQ=255)

                 gaCACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGgcgc                                                                        >  2:616322/1‑93 (MQ=255)

                  acacTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCa                                                                       <  2:4017040/93‑1 (MQ=255)

                     cTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGa                                                                    >  1:1106330/1‑93 (MQ=255)

                      tGTTCATGCTGGTTTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGaa                                                                   <  2:1836511/93‑1 (MQ=255)

                      tGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGaa                                                                   >  1:986584/1‑93 (MQ=255)

                      tGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGaa                                                                   <  1:7431161/93‑1 (MQ=255)

                      tGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGaa                                                                   <  1:605380/93‑1 (MQ=255)

                       gTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAAt                                                                  >  2:3865048/1‑93 (MQ=255)

                          cATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGta                                                               <  2:1359876/93‑1 (MQ=255)

                            tGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGtata                                                             >  2:4052666/1‑93 (MQ=255)

                            tGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGtata                                                             >  2:883035/1‑93 (MQ=255)

                            tGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGtata                                                             <  1:406200/93‑1 (MQ=255)

                            tGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGgata                                                             >  1:3678648/1‑93 (MQ=255)

                              cTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATAcc                                                           <  1:1031784/93‑1 (MQ=255)

                              cTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATAcc                                                           <  1:1431810/93‑1 (MQ=255)

                               tGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCt                                                          >  1:4347109/1‑93 (MQ=255)

                                ggTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCtt                                                         <  2:4361485/93‑1 (MQ=255)

                                ggTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCtt                                                         <  2:3407838/93‑1 (MQ=255)

                                ggTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCtt                                                         <  1:758448/93‑1 (MQ=255)

                                 gtgtCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTg                                                        >  1:5908647/1‑93 (MQ=255)

                                 gtgtCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTg                                                        >  2:2188161/1‑93 (MQ=255)

                                  tgtCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGc                                                       <  1:3523495/93‑1 (MQ=255)

                                  ggtCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGc                                                       <  2:7495839/92‑1 (MQ=255)

                                   gtCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCt                                                      <  1:3009474/93‑1 (MQ=255)

                                    tCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTg                                                     <  2:1964162/93‑1 (MQ=255)

                                     cGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGa                                                    <  1:7911033/93‑1 (MQ=255)

                                     cGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGa                                                    >  1:1559051/1‑93 (MQ=255)

                                        tGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCg                                                 <  1:4980627/93‑1 (MQ=255)

                                        tGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCg                                                 >  2:5248645/1‑93 (MQ=255)

                                         gTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGt                                                <  2:1113380/93‑1 (MQ=255)

                                         gTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGt                                                <  2:5858387/93‑1 (MQ=255)

                                         gTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGt                                                <  2:6397031/93‑1 (MQ=255)

                                            gggCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATCCCTTGCTGAGCGTGcc                                             >  2:2555756/1‑93 (MQ=255)

                                              gCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGCTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGa                                           >  1:3313765/1‑93 (MQ=255)

                                               cATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGaa                                          >  2:3626309/1‑93 (MQ=255)

                                                aTTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGaaa                                         >  1:4938568/1‑93 (MQ=255)

                                                aTTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGaaa                                         <  1:41746/93‑1 (MQ=255)

                                                 ttATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATATCTTGCTGAGCGTGCCGAAAg                                        >  1:2848521/1‑93 (MQ=255)

                                                 ttATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAg                                        >  2:556557/1‑93 (MQ=255)

                                                 ttATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAg                                        >  1:4539265/1‑93 (MQ=255)

                                                  tATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGa                                       >  2:1543600/1‑93 (MQ=255)

                                                  tATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATATATACCTTGCTGAGCGTGCCGAAAGa                                       <  1:5672310/93‑1 (MQ=255)

                                                   aTCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAt                                      >  1:5334593/1‑93 (MQ=255)

                                                   aTCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAt                                      <  1:2188044/93‑1 (MQ=255)

                                                    tCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAtg                                     <  2:3319154/93‑1 (MQ=255)

                                                    tCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAtg                                     <  2:7786629/93‑1 (MQ=255)

                                                     cAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAtgt                                    <  2:313378/93‑1 (MQ=255)

                                                     cAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAtgt                                    <  2:6122814/93‑1 (MQ=255)

                                                      aaGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGAtgtg                                   <  2:4776941/93‑1 (MQ=255)

                                                        gTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCa                                 <  2:6588968/93‑1 (MQ=255)

                                                        gTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCa                                 <  1:508918/93‑1 (MQ=255)

                                                          ttGTCGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGa                               <  2:5695686/93‑1 (MQ=255)

                                                              cGATCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGAtctt                           <  1:5851113/93‑1 (MQ=255)

                                                                aTCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGAtcttct                         >  1:3245351/1‑93 (MQ=255)

                                                                 tCAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATcttctt                        >  1:7534661/1‑93 (MQ=255)

                                                                 gaaGCTGAAAAAAGCCGGGCAGGGGGGTTACCACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATcttctt                        <  1:925450/91‑1 (MQ=255)

                                                                  cAGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCttcttc                       <  2:4941905/93‑1 (MQ=255)

                                                                   aGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTcc                      <  1:7758778/93‑1 (MQ=255)

                                                                   aGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTcc                      <  1:638908/93‑1 (MQ=255)

                                                                   aGCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTcc                      >  1:669407/1‑93 (MQ=255)

                                                                    gCTGCAAACCGCCGGGCAGGGGGATTACGCCGCGTTCGGCACAGGGATGTATACCTTGCTAAGCGTGCCGACAGATGTGCAGCTCTTCTTcca                     <  1:7021654/93‑2 (MQ=255)

                                                                    gCTGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTccc                     <  1:5500953/93‑1 (MQ=255)

                                                                      tGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGa                   <  2:3572931/93‑1 (MQ=255)

                                                                      tGAAAAAAGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGa                   >  1:6935259/1‑93 (MQ=255)

                                                                         aaaaaGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATgg                <  1:1344060/93‑1 (MQ=255)

                                                                           aaaGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATggcg              >  1:1851384/1‑93 (MQ=255)

                                                                           aaaGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATggcg              <  2:6843326/93‑1 (MQ=255)

                                                                            aaGCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATggcgg             >  1:469583/1‑93 (MQ=255)

                                                                              gCCGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGct           >  2:442687/1‑93 (MQ=255)

                                                                               ccGGGCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGctc          <  2:549579/93‑1 (MQ=255)

                                                                                 gggCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTg        >  1:4141113/1‑93 (MQ=255)

                                                                                 gggCAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTg        >  1:2259334/1‑93 (MQ=255)

                                                                                    cAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTGCtt     <  2:3032968/93‑1 (MQ=255)

                                                                                    cAGGGGGGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTGCtt     >  1:2616831/1‑93 (MQ=255)

                                                                                      ggggggTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTGCTTgg   >  2:6511958/1‑93 (MQ=255)

                                                                                      ggggggTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCAGATGGCGGCTCTGCTTgg   >  1:1533243/1‑93 (MQ=255)

                                                                                       ggggGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTGCTTggg  <  2:1238043/93‑1 (MQ=255)

TTCACCACCATCATGATGACACTGTTCATGCTGGTGTCGCTGTCGGGCATTATCAAGTTTGTCGATCAGCTGAAAAAAGCCGGGCAGGGGAGTTACGACGCGTTAGGCGCAGGAATGTATACCTTGCTGAGCGTGCCGAAAGATGTGCAGATCTTCTTCCCGATGGCGGCTCTGCTTGGG  >  NC_000913/4484931‑4485110

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