BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,520,082	A→G	H527H (CAT→CAC)	yjhG ←	putative dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,520,082	0	A	G	99.4%	529.4 / ‑4.3	166	H527H (CAT→CAC)	yjhG	putative dehydratase
Reads supporting (aligned to +/- strand): ref base A (0/0); major base G (72/93); minor base C (1/0); total (73/93)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

Only 100 of 168 total aligned reads displayed.

GGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAC  >  NC_000913/4519992‑4520171

 gggCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGc                                                                                        <  2:2718584/93‑1 (MQ=255)

 gggCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGc                                                                                        <  1:2296886/93‑1 (MQ=255)

   gCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCtt                                                                                      <  2:907928/93‑1 (MQ=255)

   gCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCtt                                                                                      <  1:3784459/93‑1 (MQ=255)

    ccTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTa                                                                                     >  1:3050553/1‑93 (MQ=255)

     cTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTAc                                                                                    >  1:3802912/1‑93 (MQ=255)

        cGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGt                                                                                 <  2:2252603/93‑1 (MQ=255)

        cGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGt                                                                                 >  1:886917/1‑93 (MQ=255)

        cGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGt                                                                                 <  2:4073120/93‑1 (MQ=255)

         ggCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTa                                                                                >  2:3219688/1‑93 (MQ=255)

         cgCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTa                                                                                <  1:3187339/92‑1 (MQ=255)

          gCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTAt                                                                               >  2:3462072/1‑93 (MQ=255)

          gCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTAt                                                                               <  2:141579/93‑1 (MQ=255)

           ccAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATg                                                                              <  1:171628/93‑1 (MQ=255)

           ccAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATg                                                                              >  1:2230778/1‑93 (MQ=255)

            cAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGa                                                                             <  2:4291092/93‑1 (MQ=255)

            cAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGa                                                                             >  2:2515324/1‑93 (MQ=255)

            cAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGa                                                                             <  1:1258897/93‑1 (MQ=255)

              gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGGTTAACGAAACGTGCTTACCGTATGACa                                                                           >  1:781201/1‑93 (MQ=255)

              gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                           >  2:3306144/1‑93 (MQ=255)

              gcgcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACa                                                                           >  1:1806920/1‑93 (MQ=255)

               cgcTTCTGGCCCCACGTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAg                                                                          <  1:1300783/93‑1 (MQ=255)

                gcTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAGCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGa                                                                         <  1:347295/93‑1 (MQ=255)

                 cTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGAt                                                                        >  1:4331798/1‑93 (MQ=255)

                 cTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGAt                                                                        <  2:1420035/93‑1 (MQ=255)

                    cTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCt                                                                     <  1:1126771/93‑1 (MQ=255)

                     tGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCtt                                                                    >  1:2946558/1‑93 (MQ=255)

                      ggCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTc                                                                   <  2:3042071/93‑1 (MQ=255)

                      ggCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTc                                                                   <  2:961405/93‑1 (MQ=255)

                      ggCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTc                                                                   <  1:4133517/93‑1 (MQ=255)

                       gCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCa                                                                  <  1:3273012/93‑1 (MQ=255)

                        ccccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAg                                                                 <  1:1173830/93‑1 (MQ=255)

                        ccccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAg                                                                 <  1:197347/93‑1 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                >  2:1475923/1‑93 (MQ=255)

                         cccACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAgg                                                                <  1:552652/93‑1 (MQ=255)

                              aTGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACt                                                           >  2:1723512/1‑93 (MQ=255)

                               tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                          >  1:3620044/1‑93 (MQ=255)

                               tGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTg                                                          <  1:3097371/93‑1 (MQ=255)

                                ggCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTgg                                                         <  1:3559954/93‑1 (MQ=255)

                                  ccGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTa                                                       >  1:2857190/1‑93 (MQ=255)

                                  ccGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTa                                                       <  1:70951/93‑1 (MQ=255)

                                   cGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTaa                                                      >  2:1271798/1‑93 (MQ=255)

                                   cGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTaa                                                      <  1:479633/93‑1 (MQ=255)

                                    gATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAc                                                     <  1:3856927/93‑1 (MQ=255)

                                     aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                    >  2:2957068/1‑93 (MQ=255)

                                     aTGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAAcc                                                    <  1:3004537/93‑1 (MQ=255)

                                      tGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCt                                                   >  2:1380433/1‑93 (MQ=255)

                                      tGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGCAACCt                                                   >  2:4388496/1‑93 (MQ=255)

                                       gCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTg                                                  <  2:2255431/93‑1 (MQ=255)

                                         ccgcgCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGt                                                <  2:556184/92‑1 (MQ=255)

                                            cgcCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAgg                                             >  2:2088527/1‑93 (MQ=255)

                                              ccAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGtt                                           >  2:1660085/1‑93 (MQ=255)

                                                cgtcgAAACGCCCAAGACACGTGCCTCGGTGGTTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttc                                         <  2:3567193/89‑1 (MQ=255)

                                                aGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttc                                         <  2:3358392/93‑1 (MQ=255)

                                                 gTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttct                                        >  2:4256065/1‑93 (MQ=255)

                                                 gTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttct                                        <  2:775739/93‑1 (MQ=255)

                                                 gTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttct                                        >  1:1372868/1‑93 (MQ=255)

                                                   aGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTttcttc                                      <  2:3379462/93‑1 (MQ=255)

                                                    gAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGCAGGTTTCGTcc                                     >  1:1012118/1‑93 (MQ=255)

                                                      aaCGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCAt                                   <  2:2575765/93‑1 (MQ=255)

                                                        cGCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATcc                                 <  2:989046/93‑1 (MQ=255)

                                                         gCCCGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATccc                                <  1:4467454/93‑1 (MQ=255)

                                                          cccGAGAATCGTGCATCGGTGAGTAACGAAACGTGCTTACCGTATGACAGATGCTACAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCt                               <  2:1834342/93‑1 (MQ=255)

                                                          cccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCt                               >  2:4649398/1‑93 (MQ=255)

                                                           ccGAGAATCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTg                              <  2:1385875/93‑1 (MQ=255)

                                                           ccGAGAAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTg                              >  2:3743937/1‑93 (MQ=255)

                                                               gaAACGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTAcc                          <  1:10116/93‑1 (MQ=255)

                                                                aaaCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCt                         >  2:1964067/1‑93 (MQ=255)

                                                                aaaCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCt                         <  2:3686843/93‑1 (MQ=255)

                                                                 aaCGTGCATCGGTGGTTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTg                        >  1:4389369/1‑93 (MQ=255)

                                                                  aCGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGa                       >  1:841237/1‑93 (MQ=255)

                                                                   cGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGaa                      <  2:1899319/93‑1 (MQ=255)

                                                                   cGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGaa                      <  1:38899/93‑1 (MQ=255)

                                                                   cGTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGaa                      <  2:2896085/93‑1 (MQ=255)

                                                                    gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                     >  2:2125220/1‑93 (MQ=255)

                                                                    gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                     <  1:1227212/93‑1 (MQ=255)

                                                                    gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                     <  2:3793955/93‑1 (MQ=255)

                                                                    gTGCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAg                     >  2:2363710/1‑93 (MQ=255)

                                                                      gCATCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGt                   <  1:4655082/93‑1 (MQ=255)

                                                                        aTCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTcc                 >  1:31063/1‑93 (MQ=255)

                                                                        aTCGGTGATTAACGAAACCTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTca                 >  2:1080896/1‑92 (MQ=255)

                                                                         tCGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCa                >  1:2447118/1‑93 (MQ=255)

                                                                          cGGTGATTAACGAAGCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               <  1:1272227/93‑1 (MQ=255)

                                                                          cGGTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCaa               <  1:3651699/93‑1 (MQ=255)

                                                                           ggTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAc              <  2:4310424/93‑1 (MQ=255)

                                                                           ggTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAc              <  1:3997126/93‑1 (MQ=255)

                                                                           ggTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAc              >  2:910981/1‑93 (MQ=255)

                                                                           ggTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAAc              <  2:1022543/93‑1 (MQ=255)

                                                                            gTGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACg             <  1:4453384/93‑1 (MQ=255)

                                                                             tGATTAACGAAACGTGGTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc            >  2:159658/1‑93 (MQ=255)

                                                                             tGATTAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGc            >  1:2591233/1‑93 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         >  1:2632257/1‑93 (MQ=255)

                                                                                ttAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaa         >  2:2827880/1‑93 (MQ=255)

                                                                                 tAACGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaat        <  2:180557/93‑1 (MQ=255)

                                                                                  aaCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaata       >  2:3714772/1‑93 (MQ=255)

                                                                                  aaCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaata       >  1:996601/1‑93 (MQ=255)

                                                                                   aCGAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCaataa      <  2:2728657/93‑1 (MQ=255)

                                                                                     gAAACGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATg    <  2:542390/93‑1 (MQ=255)

                                                                                       aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc  >  2:2204441/1‑93 (MQ=255)

                                                                                       aaCGTGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAc  >  2:1368747/1‑93 (MQ=255)

GGGGCCTCCGGCCAGCGCTTCTGGCCCCACATGGCCGATGCACGCGCCAGTAGAAACGCCCGAGAAACGTGCATCGGTGATTAACGAAACATGCTTACCGTATGACAGATGCTTCAGGGCACTGGTAACCTGGTAGGTTTCTTCCATCCCTGTACCTGAAGGTCCAACGCCAATAATGAC  >  NC_000913/4519992‑4520171

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