BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,425,773	A→G	G13G (GGT→GGC)	argT ←	ArgT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,425,773	0	A	G	99.4%	499.4 / NA	162	G13G (GGT→GGC)	argT	ArgT
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (57/104); minor base T (0/1); total (57/105)
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.99e-01

Only 100 of 170 total aligned reads displayed.

ATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTATTGTGTTGTGT  >  NC_000913/2425683‑2425862

aTCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGcc                                                                                         >  1:2872437/1‑93 (MQ=255)

 tCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCg                                                                                        >  2:2147116/1‑93 (MQ=255)

 tCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCg                                                                                        >  2:2147115/1‑93 (MQ=255)

 tCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCg                                                                                        >  2:3530994/1‑93 (MQ=255)

  cTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGa                                                                                       >  2:4184636/1‑93 (MQ=255)

  cTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGa                                                                                       <  1:314084/93‑1 (MQ=255)

   tttCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGAc                                                                                      <  1:1090099/93‑1 (MQ=255)

     tCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTa                                                                                    >  2:2690573/1‑93 (MQ=255)

       gATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAAc                                                                                  >  2:4675733/1‑93 (MQ=255)

       gATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAAc                                                                                  >  1:3500798/1‑93 (MQ=255)

       gATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAAc                                                                                  <  1:2898101/93‑1 (MQ=255)

         tGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACaa                                                                                >  1:4445984/1‑93 (MQ=255)

            gaACGTTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGa                                                                             <  2:1532049/93‑1 (MQ=255)

             aaCGGTGCGTAGGTGGTATCGGTTCCGCTACGTACCGTCTCCGGTAGCGCCGCATAGCTGGTAGGCGCTGTGGAGAGGGCGACTCACCAAGAc                                                                            >  2:2262555/1‑93 (MQ=255)

              aCGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACa                                                                           <  2:2416745/93‑1 (MQ=255)

                ggTGCGTATGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACaga                                                                         <  2:1921035/93‑1 (MQ=255)

                   gCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCg                                                                      <  2:899197/93‑1 (MQ=255)

                    cggaGGTGGTATCGGTTCCGATACTTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCga                                                                     <  1:667386/90‑1 (MQ=255)

                     gTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCgag                                                                    <  1:2782773/93‑1 (MQ=255)

                     gTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCgag                                                                    >  1:796948/1‑93 (MQ=255)

                      tAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCgaga                                                                   <  1:998855/93‑1 (MQ=255)

                      tAGGTGGTATCGGTTCCGAAACGCACCGTCTCCGGTAGCGCCGCATAGCTGGAAGGCGCTGTGGAGAGGCCGACTAACAAAGAAAGAGCgaga                                                                   >  2:2332537/1‑93 (MQ=255)

                       aGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGaa                                                                  <  2:4618047/93‑1 (MQ=255)

                       aGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGaa                                                                  >  2:1655011/1‑93 (MQ=255)

                       aGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGaa                                                                  <  1:1393943/93‑1 (MQ=255)

                        ggtggtATCGGTTCCGATATGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAAt                                                                 >  2:4567691/1‑93 (MQ=255)

                        ggtggtATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAAt                                                                 <  1:4313212/93‑1 (MQ=255)

                        ggtggtATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAAt                                                                 <  2:2889705/93‑1 (MQ=255)

                           ggtATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGa                                                              >  2:3066095/1‑93 (MQ=255)

                            gtATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAATCCTCTGTGGAGAGTCCGACTAACAAAGACAGAGCGAGAATCGAc                                                             <  2:3624955/93‑1 (MQ=255)

                             tATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGCGAGGCCGACTAACAAAGACAGAGCGAGAATCGAct                                                            >  1:1753166/1‑93 (MQ=255)

                             tATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGAct                                                            <  1:3355241/93‑1 (MQ=255)

                             tATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGAct                                                            >  2:3856709/1‑93 (MQ=255)

                               tCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGActtc                                                          <  2:1268532/93‑1 (MQ=255)

                               tCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGActtc                                                          <  1:1196633/93‑1 (MQ=255)

                               tCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGActtc                                                          >  1:484423/1‑93 (MQ=255)

                               tCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGActtc                                                          <  1:1031576/93‑1 (MQ=255)

                               tCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGCAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGActtc                                                          >  1:2189980/1‑93 (MQ=255)

                                cGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGActtct                                                         >  2:2453608/1‑93 (MQ=255)

                                  gTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACttcttc                                                       <  2:4331412/93‑1 (MQ=255)

                                   ttCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTca                                                      >  2:4207040/1‑93 (MQ=255)

                                    tCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcat                                                     <  1:959513/93‑1 (MQ=255)

                                       gATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcataca                                                  >  1:4520579/1‑93 (MQ=255)

                                       gATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcataca                                                  >  1:453967/1‑93 (MQ=255)

                                       caTACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcataca                                                  <  1:2011000/92‑1 (MQ=255)

                                        aTACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcatacat                                                 <  1:191938/93‑1 (MQ=255)

                                         tACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTcatacata                                                <  1:3816915/93‑1 (MQ=255)

                                          aCGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATaa                                               <  2:4668001/93‑1 (MQ=255)

                                             tACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAAccc                                            <  2:4430819/93‑1 (MQ=255)

                                             tACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAAccc                                            >  1:4280170/1‑93 (MQ=255)

                                             tACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAAccc                                            >  1:1488798/1‑93 (MQ=255)

                                             tACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAAccc                                            <  2:3513317/93‑1 (MQ=255)

                                              aCCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCt                                           <  1:1714211/93‑1 (MQ=255)

                                                cGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCa                                         <  1:292188/93‑1 (MQ=255)

                                                 gTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCaa                                        <  1:1247519/93‑1 (MQ=255)

                                                 gTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCaa                                        >  1:3689341/1‑93 (MQ=255)

                                                    tgcGGTAGCCCCGCCTAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTg                                     <  1:3062845/91‑1 (MQ=255)

                                                    tcCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTg                                     >  1:2108807/1‑93 (MQ=255)

                                                       ggTATCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTTaaa                                  <  1:240244/93‑1 (MQ=255)

                                                        gTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGaaaa                                 <  1:2054602/93‑1 (MQ=255)

                                                        gTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGaaaa                                 <  2:2245433/93‑1 (MQ=255)

                                                        gTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGaaaa                                 <  1:3649718/93‑1 (MQ=255)

                                                         tAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAAt                                <  1:2687595/93‑1 (MQ=255)

                                                         tAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAAt                                <  2:3801584/93‑1 (MQ=255)

                                                            ggccgcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGc                             <  1:2470475/92‑1 (MQ=255)

                                                             gccgcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCt                            <  1:1604456/93‑1 (MQ=255)

                                                             gccgcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCt                            <  2:430405/93‑1 (MQ=255)

                                                             gccgcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCt                            >  2:3167569/1‑93 (MQ=255)

                                                              ccgcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCtt                           >  2:304232/1‑93 (MQ=255)

                                                               cgcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCttt                          >  1:376266/1‑93 (MQ=255)

                                                                gcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAATTGGCtttt                         <  2:4578696/93‑1 (MQ=255)

                                                                gcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCtttt                         >  1:3818816/1‑93 (MQ=255)

                                                                gcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCtttt                         <  1:3668770/93‑1 (MQ=255)

                                                                gcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCtttt                         >  2:331181/1‑93 (MQ=255)

                                                                gcATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCtttt                         <  2:1476937/93‑1 (MQ=255)

                                                                 cATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCttttt                        <  2:3108633/93‑1 (MQ=255)

                                                                 cATAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCttttt                        <  1:3208455/93‑1 (MQ=255)

                                                                   tAGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTa                      >  2:3997105/1‑93 (MQ=255)

                                                                    aGCTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAt                     <  1:2714236/93‑1 (MQ=255)

                                                                      gtGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATACCCCTCAAGTGAAAATGGCTTTTTTAtgt                   <  2:335140/92‑1 (MQ=255)

                                                                      cTGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAtgt                   <  2:2282765/93‑1 (MQ=255)

                                                                       tGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAtgtt                  <  2:3265795/93‑1 (MQ=255)

                                                                       tGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAtgtt                  <  1:4545388/93‑1 (MQ=255)

                                                                       tGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAtgtt                  >  1:1496523/1‑93 (MQ=255)

                                                                       tGGAAGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTAtgtt                  <  2:1030375/93‑1 (MQ=255)

                                                                          aaGCCGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTtgtg               <  2:3853814/93‑1 (MQ=255)

                                                                             ccGCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAt            <  2:3940449/93‑1 (MQ=255)

                                                                               gCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttg          >  2:2924554/1‑93 (MQ=255)

                                                                               gCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttg          >  1:3943005/1‑93 (MQ=255)

                                                                               gCTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttg          <  2:1276029/93‑1 (MQ=255)

                                                                                cTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgt         <  1:125146/93‑1 (MQ=255)

                                                                                cTGTGGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTGTTATGTTGTGTAttgt         <  2:4046646/93‑1 (MQ=255)

                                                                                  gtgGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgtgt       <  1:2331048/93‑1 (MQ=255)

                                                                                  gtgGAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgtgt       <  1:900869/93‑1 (MQ=255)

                                                                                    ggAGAGGCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgtgttg     <  2:636478/93‑1 (MQ=255)

                                                                                     gagaggCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgtgttgt    >  1:3595969/1‑93 (MQ=255)

                                                                                     gagaggCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgtgttgt    <  2:3616190/93‑1 (MQ=255)

                                                                                     gagaggCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTAttgtgttgt    <  2:2060774/93‑1 (MQ=255)

                                                                                       gaggCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTATtgtgttgtgt  <  1:4038715/93‑1 (MQ=255)

                                                                                       gaggCCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTATtgtgttgtgt  >  1:3659089/1‑93 (MQ=255)

ATCTTTCGATGAGAACGGTGCGTAGGTGGTATCGGTTCCGATACGTACCGTCTCCGGTAGCGCCGCATAGCTGGAAGCCGCTGTGGAGAGACCGACTAACAAAGACAGAGCGAGAATCGACTTCTTCATACATAACCCTCAAGTGAAAATGGCTTTTTTATGTTGTGTATTGTGTTGTGT  >  NC_000913/2425683‑2425862

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