BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,343,880	A→G	V296A (GTG→GCG)	yciR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,343,880	0	A	G	98.2%	521.7 / ‑1.2	167	V296A (GTG→GCG)	yciR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (87/77); total (88/79)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

Only 100 of 176 total aligned reads displayed.

GAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCAACT  >  NC_000913/1343790‑1343969

  ggTGTTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACCCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTc                                                                                       <  1:3926173/93‑1 (MQ=255)

   gtgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcg                                                                                      >  1:2211437/1‑93 (MQ=255)

   gtgtTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACCCCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTcg                                                                                      <  1:4408754/93‑1 (MQ=255)

    tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                                                                                     >  2:4063619/1‑93 (MQ=255)

    tgtTGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTcgc                                                                                     >  2:3704123/1‑93 (MQ=255)

      ttGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcgtcgcgt                                                                                   >  1:211419/1‑93 (MQ=255)

      gtGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGAcgcgtcgcgt                                                                                   <  1:1007347/92‑1 (MQ=255)

       tGGATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGTCAACACCTGGTCATGTTGGAGACAGCTTAAAATAGCCAATGACGCGTCGCGta                                                                                  <  1:277740/93‑1 (MQ=255)

         gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAGCACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                                                                                <  2:2677291/93‑1 (MQ=255)

         gATGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCTAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaat                                                                                <  2:2445167/93‑1 (MQ=255)

          aTGCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGtaata                                                                               >  1:1295272/1‑93 (MQ=255)

            gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                                                                             <  2:2957054/93‑1 (MQ=255)

            gCCAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAAc                                                                             <  1:2593472/93‑1 (MQ=255)

             ccAGTACCAGAAACTCATCCTCACCTGGACGCTCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACt                                                                            <  1:2323880/93‑1 (MQ=255)

              cAGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTg                                                                           <  2:2853839/93‑1 (MQ=255)

               aGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTgg                                                                          >  2:3157887/1‑93 (MQ=255)

               aGTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTgg                                                                          <  1:205192/93‑1 (MQ=255)

                gTACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGt                                                                         <  2:1136503/93‑1 (MQ=255)

                 tACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTc                                                                        <  2:21007/93‑1 (MQ=255)

                 tACCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTc                                                                        <  1:3747184/93‑1 (MQ=255)

                  aCCAGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCa                                                                       <  2:1059050/93‑1 (MQ=255)

                     aGAAACTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACca                                                                    >  1:3078876/1‑93 (MQ=255)

                        aaCTCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaac                                                                 <  2:3345857/93‑1 (MQ=255)

                           tCATCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACcaaacaaa                                                              <  2:1676427/93‑1 (MQ=255)

                             aTCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATg                                                            >  1:1703405/1‑93 (MQ=255)

                             aTCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGAAATAACTGGTCACCAAACAAATg                                                            >  1:4236061/1‑93 (MQ=255)

                              tCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGc                                                           >  2:3433024/1‑93 (MQ=255)

                              tCCTCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGc                                                           >  1:1667789/1‑93 (MQ=255)

                                 tcACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCa                                                        >  1:4144264/2‑93 (MQ=255)

                                 tcACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCa                                                        >  1:3705753/2‑93 (MQ=255)

                                 tcACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCa                                                        >  1:2856652/2‑93 (MQ=255)

                                  cACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCAt                                                       <  2:144148/93‑1 (MQ=255)

                                  cACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCAt                                                       <  2:4399177/93‑1 (MQ=255)

                                  cACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCAt                                                       >  1:1640220/1‑93 (MQ=255)

                                   cccTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATa                                                      <  1:2040419/92‑1 (MQ=255)

                                   aCCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATa                                                      >  1:3997320/1‑93 (MQ=255)

                                   aCCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATa                                                      >  2:2258989/1‑93 (MQ=255)

                                    ccTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAg                                                     >  1:1589744/1‑93 (MQ=255)

                                    ccTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGCCACCAAACAAATGCCCATAg                                                     >  2:171808/1‑93 (MQ=255)

                                     cTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAgg                                                    >  2:1262397/1‑93 (MQ=255)

                                     cTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAgg                                                    <  2:3909428/93‑1 (MQ=255)

                                     cTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAgg                                                    >  2:3115451/1‑93 (MQ=255)

                                      tGGACGCGCCAACATCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGc                                                   >  2:2633634/1‑93 (MQ=255)

                                      tGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGc                                                   >  1:452415/1‑93 (MQ=255)

                                       ggACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcg                                                  <  2:38971/93‑1 (MQ=255)

                                       ggACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcg                                                  <  1:4191302/93‑1 (MQ=255)

                                        gACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgt                                                 <  1:1785940/93‑1 (MQ=255)

                                         aCGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtc                                                <  1:734558/93‑1 (MQ=255)

                                         aCGCGCCAACACCTGGTCATGTTCGAGACAGCTGAAAATAGCCAATGACGCGCCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtc                                                <  1:4566077/93‑1 (MQ=255)

                                          cgcgCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtcg                                               >  2:3464749/1‑93 (MQ=255)

                                           gcgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtcgt                                              <  2:1251243/93‑1 (MQ=255)

                                           gcgcCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGcgtcgt                                              <  2:3449497/93‑1 (MQ=255)

                                              ccAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGa                                           >  1:634332/1‑93 (MQ=255)

                                               cAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGAc                                          >  1:1179288/1‑93 (MQ=255)

                                                aaCACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGAcc                                         >  2:3868622/1‑93 (MQ=255)

                                                aaCACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGAcc                                         >  1:1043135/1‑93 (MQ=255)

                                                 acacCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCt                                        >  2:44139/1‑93 (MQ=255)

                                                     cTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCttttt                                    >  1:3225607/1‑93 (MQ=255)

                                                      tGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTg                                   >  1:2286949/1‑93 (MQ=255)

                                                       ggTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGa                                  >  1:4116813/1‑93 (MQ=255)

                                                       ggTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGa                                  >  1:3760117/1‑93 (MQ=255)

                                                       ggTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGa                                  >  2:3452054/1‑93 (MQ=255)

                                                        gTCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGaa                                 >  2:4729507/1‑93 (MQ=255)

                                                         tCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGaaa                                >  2:4866/1‑93 (MQ=255)

                                                         tCATGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGaaa                                >  2:2566871/1‑93 (MQ=255)

                                                           aTGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAAtt                              >  2:2170746/1‑93 (MQ=255)

                                                            tGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTa                             <  2:2064380/93‑1 (MQ=255)

                                                            tGTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTa                             >  2:3336929/1‑93 (MQ=255)

                                                             gTTCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTAt                            <  2:653303/93‑1 (MQ=255)

                                                              ttCGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATc                           <  2:4105731/93‑1 (MQ=255)

                                                                cGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCa                         >  1:264754/1‑93 (MQ=255)

                                                                cGAGACAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCa                         <  2:555997/93‑1 (MQ=255)

                                                                 gagaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaa                        <  2:151/93‑1 (MQ=255)

                                                                  agaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaaa                       >  2:448759/1‑93 (MQ=255)

                                                                  agaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaaa                       >  2:2000095/1‑93 (MQ=255)

                                                                  agaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaaa                       <  1:1528310/93‑1 (MQ=255)

                                                                  agaCAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCaaa                       <  1:4366165/93‑1 (MQ=255)

                                                                   gaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAAt                      <  2:2345164/93‑1 (MQ=255)

                                                                   gaCAGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAAt                      >  2:3565663/1‑93 (MQ=255)

                                                                      aGCTTCAAATAGCCAATTACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCaa                   <  2:3960838/93‑1 (MQ=255)

                                                                      aGCTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCaa                   <  1:45980/93‑1 (MQ=255)

                                                                        cTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGa                 <  1:383258/93‑1 (MQ=255)

                                                                        cTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGa                 <  1:790273/93‑1 (MQ=255)

                                                                        cTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGa                 <  1:4662454/93‑1 (MQ=255)

                                                                        cTTAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGa                 <  1:2849879/93‑1 (MQ=255)

                                                                          tAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATa               <  1:3280318/93‑1 (MQ=255)

                                                                          tAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATa               <  1:3298894/93‑1 (MQ=255)

                                                                          tAAAATAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATa               <  1:4152460/93‑1 (MQ=255)

                                                                          tAAAAAAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCACAGAGGCGGCGTTGACCTTTTTGAAATTATCCAAATCAAGATa               <  2:1372554/93‑1 (MQ=255)

                                                                            aaaTAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATaca             >  2:138650/1‑93 (MQ=255)

                                                                            aaaTAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATaca             <  2:2262603/93‑1 (MQ=255)

                                                                               tAGCCAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACaa          >  1:3802010/1‑93 (MQ=255)

                                                                                   cAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAAcccc      >  1:1866712/1‑93 (MQ=255)

                                                                                   cAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAAcccc      <  1:1569249/93‑1 (MQ=255)

                                                                                   cAATGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAAcccc      >  2:2742999/1‑93 (MQ=255)

                                                                                    aaTGACAAGACGCGTAATAACTCGGAACCAGACAAATGCACATAGGCGTCGTTGACCTTTTTGAAATTATACACATCAAGATACACAACCCCa     >  2:2825137/1‑93 (MQ=255)

                                                                                     aTGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCaa    <  2:2054999/93‑1 (MQ=255)

                                                                                     aTGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCaa    >  1:3656585/1‑93 (MQ=255)

                                                                                     aTGACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCaa    >  2:1776820/1‑93 (MQ=255)

                                                                                       gACGCGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCACCt  >  1:3190545/1‑93 (MQ=255)

GAGGTGTTGGATGCCAGTACCAGAAACTCATCCCCACCTGGACGCGCCAACACCTGGTCATGTTCGAGACAGCTTAAAATAGCCAATGACACGTCGCGTAATAACTGGTCACCAAACAAATGCCCATAGGCGTCGTTGACCTTTTTGAAATTATCCAAATCAAGATACACAACCCCAACT  >  NC_000913/1343790‑1343969

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