BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,201,958	G→A	W354* (TGG→TAG)	hydG →	HydG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,201,958	0	G	A	99.7%	1132.1 / ‑6.1	311	W354* (TGG→TAG)	hydG	HydG
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (159/150); total (160/151)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 317 total aligned reads displayed.

GCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGC  >  NC_000913/4201868‑4202047

  cATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCg                                                                                       >  2:4358531/1‑93 (MQ=255)

  cATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCg                                                                                       >  1:2292155/1‑93 (MQ=255)

   aTTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCgg                                                                                      >  1:7639841/1‑93 (MQ=255)

    ttttCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCggg                                                                                     >  1:6120296/1‑93 (MQ=255)

    ttttCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCggg                                                                                     >  1:6385061/1‑93 (MQ=255)

          gCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAatat                                                                               >  1:1816156/1‑93 (MQ=255)

          gCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAatat                                                                               <  2:5729190/93‑1 (MQ=255)

          gCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAatat                                                                               <  1:6909340/93‑1 (MQ=255)

           cAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATAtt                                                                              <  2:2808447/93‑1 (MQ=255)

            aGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTc                                                                             <  2:1419465/93‑1 (MQ=255)

            aGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTc                                                                             <  2:256732/93‑1 (MQ=255)

             gcgcTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCg                                                                            <  1:3290578/93‑1 (MQ=255)

             gcgcTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCg                                                                            >  2:396570/1‑93 (MQ=255)

             gcgcTTTGCCGAGCGTAATCGCAAGGCGGGAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATAttgg                                                                            >  2:3162861/1‑91 (MQ=255)

                cTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGa                                                                         <  2:6253862/93‑1 (MQ=255)

                cTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGa                                                                         >  2:2763322/1‑93 (MQ=255)

                cTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGa                                                                         <  1:781745/93‑1 (MQ=255)

                 tttGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAg                                                                        <  2:636447/93‑1 (MQ=255)

                  ttGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGc                                                                       >  1:6457262/1‑93 (MQ=255)

                  ttGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGc                                                                       <  1:4583266/93‑1 (MQ=255)

                    gCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTg                                                                     <  2:7910114/93‑1 (MQ=255)

                      cGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGa                                                                   >  1:5066924/1‑93 (MQ=255)

                       gAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGaa                                                                  >  1:1102802/1‑93 (MQ=255)

                       gAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGaa                                                                  >  1:4931566/1‑93 (MQ=255)

                       gAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGaa                                                                  >  2:4020275/1‑93 (MQ=255)

                         gCGTAATCGCAAGGCGGTAAAAGGTTTTCCGCCCCAGGAGATGGATCTGTTGATTCATTACGTCTAGCCGGGAAATATTCGTGAGCTGGaaaa                                                                <  1:3516814/93‑1 (MQ=255)

                         gCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGacaa                                                                >  2:7766597/1‑93 (MQ=255)

                         gCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGaaaa                                                                <  2:3456655/93‑1 (MQ=255)

                         gCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGaaaa                                                                <  1:6651807/93‑1 (MQ=255)

                          catAATCGCAAGGCGGTAAAAGGTTATACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAc                                                               <  2:7362495/91‑1 (MQ=255)

                            tAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgc                                                             >  1:4232533/1‑93 (MQ=255)

                            tAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgc                                                             >  2:5557995/1‑93 (MQ=255)

                             aaTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgcg                                                            >  1:7436290/1‑93 (MQ=255)

                             aaTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgcg                                                            >  2:4887732/1‑93 (MQ=255)

                             aaTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgcg                                                            <  1:6955964/93‑1 (MQ=255)

                             aaTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgcg                                                            >  1:5589836/1‑93 (MQ=255)

                             aaTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAAcgcg                                                            <  1:5742939/93‑1 (MQ=255)

                              aTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCgg                                                           >  1:7888150/1‑93 (MQ=255)

                              aTCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCgg                                                           <  2:877028/93‑1 (MQ=255)

                               tCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGt                                                          >  2:369629/1‑93 (MQ=255)

                                 gCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTgg                                                        >  2:2491631/1‑93 (MQ=255)

                                  cAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGa                                                       >  2:4697325/1‑93 (MQ=255)

                                   aaGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGaa                                                      >  1:1005601/1‑93 (MQ=255)

                                   aaGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGaa                                                      <  2:1734714/93‑1 (MQ=255)

                                   aaGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGaa                                                      <  1:5342418/93‑1 (MQ=255)

                                   aaGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGaa                                                      <  1:94552/93‑1 (MQ=255)

                                      gCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACgg                                                   <  1:1222243/93‑1 (MQ=255)

                                        ggTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGc                                                 >  2:3366437/1‑93 (MQ=255)

                                          tAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAg                                               <  1:5844358/93‑1 (MQ=255)

                                           aaaaGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAgt                                              <  2:1508289/93‑1 (MQ=255)

                                           aaaaGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAgt                                              <  2:2188536/93‑1 (MQ=255)

                                           aaaaGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAgt                                              <  1:7765813/93‑1 (MQ=255)

                                            aaaGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAgtg                                             <  1:2360780/93‑1 (MQ=255)

                                             aaGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAgtgg                                            <  2:4551988/93‑1 (MQ=255)

                                                gTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGtgc                                         <  1:4308090/93‑1 (MQ=255)

                                                gTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGtgc                                         <  2:5834480/93‑1 (MQ=255)

                                                gTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGtgc                                         <  1:6065689/93‑1 (MQ=255)

                                                 ttttACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGtgct                                        >  1:6240039/1‑93 (MQ=255)

                                                  tttACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGtgctg                                       >  1:477322/1‑93 (MQ=255)

                                                    tACGCCCCAGGCGCTGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTgctgct                                     <  1:5434058/93‑1 (MQ=255)

                                                         cccAGGCGTTGGATCTGTTGATTCGTTACGCCTCGCCGGGAGATATTCGTGAGCTGGAAACCGCGGTGGAACGGGCAGTGGTGCTGCTGACCg                                <  1:7017856/93‑1 (MQ=255)

                                                          ccAGGCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCgg                               >  2:5950062/1‑93 (MQ=255)

                                                              gggATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGaa                           <  2:3128605/91‑1 (MQ=255)

                                                              gCGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGaa                           <  1:4376610/93‑1 (MQ=255)

                                                               cGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAat                          <  1:3608212/93‑1 (MQ=255)

                                                               cGATGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAat                          <  1:4687316/93‑1 (MQ=255)

                                                                 aTGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatat                        <  2:2607198/93‑1 (MQ=255)

                                                                 aTGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatat                        >  1:152536/1‑93 (MQ=255)

                                                                 aTGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatat                        <  1:4468948/93‑1 (MQ=255)

                                                                  tGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatata                       <  2:5468206/93‑1 (MQ=255)

                                                                  tGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatata                       >  2:5056245/1‑93 (MQ=255)

                                                                  tGGATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatata                       >  1:2236538/1‑93 (MQ=255)

                                                                   ggATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAatatat                      <  2:7186028/93‑1 (MQ=255)

                                                                    gATCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATAtt                     <  2:1081651/93‑1 (MQ=255)

                                                                      tCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTc                   >  2:5516363/1‑93 (MQ=255)

                                                                      tCTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTc                   >  2:2087095/1‑93 (MQ=255)

                                                                       cTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTcc                  <  2:3916487/93‑1 (MQ=255)

                                                                       cTGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTcc                  <  1:3862217/93‑1 (MQ=255)

                                                                        tGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCg                 >  1:4239616/1‑93 (MQ=255)

                                                                        tGTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCg                 >  1:5987810/1‑93 (MQ=255)

                                                                         gTTGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGa                >  1:4053836/1‑93 (MQ=255)

                                                                          ttGATTCATTACGACTAGCCGTGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGaa               <  2:6777900/93‑1 (MQ=255)

                                                                           tGATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAc              >  2:797054/1‑93 (MQ=255)

                                                                            gATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAcg             >  2:5270762/1‑93 (MQ=255)

                                                                            gATTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAcg             >  2:1859069/1‑93 (MQ=255)

                                                                             aTTCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAcgc            >  1:282057/1‑93 (MQ=255)

                                                                              ttCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAAcgcg           >  2:2879126/1‑93 (MQ=255)

                                                                               tCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGa          <  2:7583831/93‑1 (MQ=255)

                                                                               tCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGa          >  2:1206187/1‑93 (MQ=255)

                                                                               tCATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGa          <  1:3753438/93‑1 (MQ=255)

                                                                                cATTACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAg         >  1:3127794/1‑93 (MQ=255)

                                                                                  ttACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCTAGCt       >  1:3047905/1‑93 (MQ=255)

                                                                                  ttACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCt       >  2:6850811/1‑93 (MQ=255)

                                                                                  ttACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCt       <  1:1888820/93‑1 (MQ=255)

                                                                                  ttACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCt       <  1:832316/93‑1 (MQ=255)

                                                                                  ttACGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCt       >  2:2708901/1‑93 (MQ=255)

                                                                                    aCGACTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGc     <  1:4979500/93‑1 (MQ=255)

                                                                                       aCTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGc  >  1:2632053/1‑93 (MQ=255)

                                                                                       aCTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGc  <  1:4607469/93‑1 (MQ=255)

                                                                                       aCTAGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGc  <  2:2294702/93‑1 (MQ=255)

GCCATTTTCTGCAGCGCTTTGCCGAGCGTAATCGCAAGGCGGTAAAAGGTTTTACGCCCCAGGCGATGGATCTGTTGATTCATTACGACTGGCCGGGAAATATTCGTGAGCTGGAAAACGCGGTGGAACGGGCAGTGGTGCTGCTGACCGGGGAATATATTTCCGAACGCGAGCTGCCGC  >  NC_000913/4201868‑4202047

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