BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,388,932	A→G	A52A (GCT→GCC)	b4161 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,388,932	0	A	G	99.6%	694.5 / NA	221	A52A (GCT→GCC)	b4161	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (105/114); minor base C (1/0); total (107/114)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 237 total aligned reads displayed.

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                         >  2:5490476/1‑93 (MQ=255)

aCGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGc                                                                                         >  1:760985/1‑93 (MQ=255)

 cGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCg                                                                                        <  2:6054038/93‑1 (MQ=255)

  gCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGt                                                                                       <  1:6192753/93‑1 (MQ=255)

    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGc                                                                                     >  1:4327695/1‑93 (MQ=255)

    cAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGc                                                                                     >  2:5539778/1‑93 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCTGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   <  2:2107840/93‑1 (MQ=255)

      gACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCAt                                                                                   <  1:2866626/93‑1 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATa                                                                                  >  1:3846816/1‑93 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATa                                                                                  >  1:1725802/1‑93 (MQ=255)

       aCTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATa                                                                                  <  1:639615/93‑1 (MQ=255)

        cTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAc                                                                                 <  1:5073347/93‑1 (MQ=255)

         tACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAcc                                                                                >  1:4035217/1‑93 (MQ=255)

          aCGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCa                                                                               >  2:737771/1‑93 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              <  2:1051923/93‑1 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCaa                                                                              <  2:4242240/93‑1 (MQ=255)

           cgcgGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGGGCATACCaa                                                                              >  1:2612133/1‑93 (MQ=255)

              ggTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGc                                                                           <  2:5783128/93‑1 (MQ=255)

               gTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCg                                                                          <  1:2208140/93‑1 (MQ=255)

               gTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATAACAAAGCg                                                                          >  2:175726/1‑93 (MQ=255)

                tCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCgg                                                                         >  1:5914655/1‑93 (MQ=255)

                 cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  1:1340449/93‑1 (MQ=255)

                 cGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGc                                                                        <  2:4503158/93‑1 (MQ=255)

                   ccGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtg                                                                      >  2:3903585/1‑93 (MQ=255)

                    cGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtga                                                                     <  1:5528069/93‑1 (MQ=255)

                     ggTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgat                                                                    <  2:5372259/93‑1 (MQ=255)

                      gTTGCCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  1:5504032/93‑1 (MQ=255)

                      gTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatg                                                                   <  1:1877452/93‑1 (MQ=255)

                       ttACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCtgatga                                                                  <  2:1296232/93‑1 (MQ=255)

                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAc                                                                 >  1:2097538/1‑93 (MQ=255)

                        tACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAc                                                                 <  1:5158095/93‑1 (MQ=255)

                          ccAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGa                                                               <  2:815324/93‑1 (MQ=255)

                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                              <  2:5353992/93‑1 (MQ=255)

                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAt                                                              >  2:3228917/1‑93 (MQ=255)

                           cAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGAAGAt                                                              >  1:4652806/1‑93 (MQ=255)

                            aGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATa                                                             <  2:2233248/93‑1 (MQ=255)

                             gCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAc                                                            >  1:2960088/1‑93 (MQ=255)

                              cGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGAGAcc                                                           >  1:5261189/1‑93 (MQ=255)

                                aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                         >  2:5820547/1‑93 (MQ=255)

                                aaCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCtt                                                         <  1:4278178/93‑1 (MQ=255)

                                 aCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtc                                                        <  2:2167902/93‑1 (MQ=255)

                                 aCGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtc                                                        <  1:737771/93‑1 (MQ=255)

                                  cGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtca                                                       <  1:5980569/93‑1 (MQ=255)

                                    gATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTtcatc                                                     <  1:3446849/93‑1 (MQ=255)

                                       ggTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAgg                                                  <  1:423462/93‑1 (MQ=255)

                                       ggTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAgg                                                  <  1:3080325/93‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 >  1:5889129/1‑93 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:2935237/93‑1 (MQ=255)

                                        gTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGc                                                 <  2:2598534/93‑1 (MQ=255)

                                         tACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCt                                                <  2:4037633/93‑1 (MQ=255)

                                         tACGGCGAATATTGCAGCGGGGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGAGGCTGAACACCAGACCGTCACTAGGCt                                                >  1:4866051/1‑93 (MQ=255)

                                          aCGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTc                                               >  2:5432372/1‑93 (MQ=255)

                                           cGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCg                                              >  2:1855312/1‑93 (MQ=255)

                                            ggCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGc                                             >  1:5431409/1‑93 (MQ=255)

                                             gCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGcc                                            >  2:1959271/1‑93 (MQ=255)

                                             gCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGcc                                            <  1:4185865/93‑1 (MQ=255)

                                             gCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGcc                                            >  2:3440996/1‑93 (MQ=255)

                                               gAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaa                                          >  1:3295436/1‑93 (MQ=255)

                                                aaTATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCaaa                                         >  1:4811888/1‑93 (MQ=255)

                                                 atatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAAc                                        >  2:4577958/1‑93 (MQ=255)

                                                 atatTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAAc                                        <  1:476622/93‑1 (MQ=255)

                                                   atTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCCCCAAACAg                                      >  1:4535212/1‑93 (MQ=255)

                                                    ttGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGa                                     <  1:1386465/93‑1 (MQ=255)

                                                     tGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGAt                                    <  1:3221429/93‑1 (MQ=255)

                                                      gCAGCGGTCAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGAAGACGATACCTTCATCAGGCTCGACAAACACAtt                                   >  2:6137844/1‑93 (MQ=255)

                                                        aGCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgt                                 <  2:1615792/93‑1 (MQ=255)

                                                         gCGGTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtc                                <  2:2377119/93‑1 (MQ=255)

                                                           ggTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATTAAGATACCGTCATCAGGCTCGCCAAACAGATTgtcgt                              >  2:2535739/1‑93 (MQ=255)

                                                           ggTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgt                              <  1:3420220/93‑1 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                             <  1:884211/93‑1 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                             <  2:193045/93‑1 (MQ=255)

                                                            gTGAACGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTgtcgtc                             >  2:1804866/1‑93 (MQ=255)

                                                               aaCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTa                          <  1:4940043/93‑1 (MQ=255)

                                                               aaCGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGAGGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTa                          <  1:4439781/93‑1 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  2:4409805/93‑1 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  2:488834/93‑1 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        >  1:6221329/1‑93 (MQ=255)

                                                                 cGTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGt                        <  2:2736243/93‑1 (MQ=255)

                                                                  gTCGCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATAACTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTc                       >  2:596560/1‑93 (MQ=255)

                                                                     gCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCggg                    >  1:6293840/1‑93 (MQ=255)

                                                                     gCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCggg                    >  2:5474999/1‑93 (MQ=255)

                                                                     gCCATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCggg                    >  1:479304/1‑93 (MQ=255)

                                                                      ccATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGc                   >  1:5854863/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                  >  2:5160441/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCt                  >  1:2387617/1‑93 (MQ=255)

                                                                       cATCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTAGTCGTAGTCGGGCt                  >  1:2791443/1‑93 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc                >  2:4657493/1‑93 (MQ=255)

                                                                         tCGGCGGATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTtc                <  1:1235884/93‑1 (MQ=255)

                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             <  2:5550876/93‑1 (MQ=255)

                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             >  2:5803149/1‑93 (MQ=255)

                                                                            gcggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTcc             <  1:6083737/93‑1 (MQ=255)

                                                                              ggATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCtt           >  1:2908828/1‑93 (MQ=255)

                                                                               gATTCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTa          <  1:1083281/93‑1 (MQ=255)

                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        <  2:3782895/93‑1 (MQ=255)

                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        >  2:1118278/1‑93 (MQ=255)

                                                                                 ttCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGa        <  2:5358802/93‑1 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGTTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       <  2:4111167/93‑1 (MQ=255)

                                                                                  tCCACATCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGAc       >  2:1212881/1‑93 (MQ=255)

                                                                                      caTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGttt   <  1:3714109/93‑1 (MQ=255)

                                                                                       aTCGGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGtttt  >  1:5532093/1‑93 (MQ=255)

ACGCCAGACTACGCGGTCGCCGGTTACCAGCGAACGGATGGTACGGCGAATATTGCAGCGGTGAACGTCGCCATCGGCGGATTCCACATCAGCGTGCATACCAAAGCGGCTGATGACGATACCTTCATCAGGCTCGCCAAACAGATTGTCGTCGTAGTCGGGCTTCTCCTTAGACGTTTT  >  NC_000913/4388842‑4389021

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