BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,344,524	T→C	G81G (GGA→GGG)	yciR ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,344,524	0	T	C	99.0%	299.0 / ‑4.8	100	G81G (GGA→GGG)	yciR	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (51/47); total (53/47)
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.13e-01

Only 100 of 105 total aligned reads displayed.

GGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAG  >  NC_000913/1344434‑1344612

ggAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGccc                                                                                        >  1:4568086/1‑93 (MQ=255)

ggAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGccc                                                                                        >  1:3250224/1‑93 (MQ=255)

ggAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGccc                                                                                        >  1:2616031/1‑93 (MQ=255)

  aaaGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCtc                                                                                      >  1:832280/1‑92 (MQ=255)

  aaaGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTa                                                                                      >  1:483356/1‑93 (MQ=255)

   aaGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTAc                                                                                     <  2:1564790/93‑1 (MQ=255)

     gCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCa                                                                                   >  2:2633303/1‑93 (MQ=255)

     gCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCa                                                                                   <  2:2365228/93‑1 (MQ=255)

     gCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCa                                                                                   <  1:2037498/93‑1 (MQ=255)

     gCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCa                                                                                   >  2:100464/1‑93 (MQ=255)

     gCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCa                                                                                   >  1:3517746/1‑93 (MQ=255)

        cGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGAt                                                                                <  2:3825387/93‑1 (MQ=255)

         gtgtGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATg                                                                               >  2:4431121/1‑93 (MQ=255)

          tgtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGc                                                                              >  1:3893430/1‑93 (MQ=255)

           ttgATAAATCACGAGCACCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGca                                                                             <  1:1932889/92‑1 (MQ=255)

           gtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGca                                                                             >  2:491184/1‑93 (MQ=255)

            tgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGcac                                                                            >  1:847704/1‑93 (MQ=255)

             gATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGcaca                                                                           >  2:1688570/1‑93 (MQ=255)

             gATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGcaca                                                                           >  2:815485/1‑93 (MQ=255)

             gATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGcaca                                                                           <  1:2617545/93‑1 (MQ=255)

             gATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGcaca                                                                           <  2:3863517/93‑1 (MQ=255)

              aTAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAg                                                                          >  2:267612/1‑93 (MQ=255)

               tAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGa                                                                         >  2:3269701/1‑93 (MQ=255)

               aaaaTCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGa                                                                         <  2:2923327/92‑1 (MQ=255)

                    cACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTTGCCCTACCAGATGCACAGAAAGAt                                                                    >  1:3673977/1‑93 (MQ=255)

                    cACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGAt                                                                    <  2:3477231/93‑1 (MQ=255)

                    cACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGAt                                                                    <  2:4508623/93‑1 (MQ=255)

                     aCGATCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATc                                                                   <  1:2749491/93‑1 (MQ=255)

                       gAGCACCCGCCGGTGTGGCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGc                                                                 <  1:2372063/93‑1 (MQ=255)

                       gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGc                                                                 <  2:1368591/93‑1 (MQ=255)

                       gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGc                                                                 >  2:3547241/1‑93 (MQ=255)

                       gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGc                                                                 <  2:472209/93‑1 (MQ=255)

                       gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGc                                                                 >  1:1031891/1‑93 (MQ=255)

                        aGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCt                                                                >  1:102094/1‑93 (MQ=255)

                        aGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCt                                                                <  1:4632072/93‑1 (MQ=255)

                          cAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCttt                                                              <  1:4452094/93‑1 (MQ=255)

                           aaCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTc                                                             >  1:4274875/1‑93 (MQ=255)

                           aaCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTc                                                             <  1:449502/93‑1 (MQ=255)

                           aaCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTc                                                             >  2:2969712/1‑93 (MQ=255)

                           aaCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTc                                                             >  2:629306/1‑93 (MQ=255)

                            aCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTct                                                            >  1:4154532/1‑92 (MQ=255)

                            aCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCa                                                            >  1:2915337/1‑93 (MQ=255)

                            aCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCa                                                            >  2:503282/1‑93 (MQ=255)

                            aCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCa                                                            <  2:3602673/93‑1 (MQ=255)

                             ccgccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCAt                                                           >  1:3233909/1‑93 (MQ=255)

                              cgccgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATc                                                          <  1:491986/93‑1 (MQ=255)

                                 cgGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGAc                                                       <  2:3307007/93‑1 (MQ=255)

                                    tgtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGt                                                    <  1:2246507/93‑1 (MQ=255)

                                     gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTc                                                   <  1:3624386/93‑1 (MQ=255)

                                     gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTc                                                   <  1:3624385/93‑1 (MQ=255)

                                     gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTc                                                   <  1:369511/93‑1 (MQ=255)

                                     gtgtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTc                                                   >  1:3228560/1‑93 (MQ=255)

                                       gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAg                                                 >  1:436547/1‑93 (MQ=255)

                                       gtCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAg                                                 >  2:404208/1‑93 (MQ=255)

                                        tCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGa                                                >  2:839667/1‑93 (MQ=255)

                                         catgccatgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGAc                                               <  1:749233/93‑1 (MQ=255)

                                         catgccatgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGAc                                               <  2:2675172/93‑1 (MQ=255)

                                         catgccatgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCACCATACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGAc                                               <  2:754181/93‑1 (MQ=255)

                                          atgccatgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACt                                              >  1:1356835/1‑93 (MQ=255)

                                               atgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatca                                         >  1:4157005/1‑93 (MQ=255)

                                               atgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatca                                         <  1:1043318/93‑1 (MQ=255)

                                                tgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat                                        >  1:2620702/1‑93 (MQ=255)

                                                tgcAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcat                                        <  1:3390057/93‑1 (MQ=255)

                                                  cAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcatca                                      <  1:813728/93‑1 (MQ=255)

                                                  cAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACtcatcatca                                      <  1:2566182/93‑1 (MQ=255)

                                                     aCGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGc                                   >  2:2507368/1‑93 (MQ=255)

                                                     aCGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGc                                   <  1:3759631/93‑1 (MQ=255)

                                                      cgcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCt                                  <  2:1446219/93‑1 (MQ=255)

                                                      cgcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCt                                  <  2:518967/93‑1 (MQ=255)

                                                         gTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGa                               <  1:3267536/93‑1 (MQ=255)

                                                          ttGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAg                              >  1:1350747/1‑93 (MQ=255)

                                                          ttGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAg                              >  1:1265876/1‑93 (MQ=255)

                                                             ccAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtg                           <  2:2991582/93‑1 (MQ=255)

                                                              cAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtga                          >  2:873843/1‑93 (MQ=255)

                                                               cgCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtgat                         <  2:3675586/92‑1 (MQ=255)

                                                                gCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtgatg                        >  1:3633234/1‑93 (MQ=255)

                                                                gCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGtgatg                        >  2:279224/1‑93 (MQ=255)

                                                                   cATTCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGAcc                     >  2:380695/1‑93 (MQ=255)

                                                                      tCCCGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTc                  <  2:3166700/93‑1 (MQ=255)

                                                                       cccGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCa                 <  2:4442128/93‑1 (MQ=255)

                                                                        ccGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAt                >  1:1571914/1‑93 (MQ=255)

                                                                        ccGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAt                >  2:2943877/1‑93 (MQ=255)

                                                                         cGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAtt               >  1:81279/1‑93 (MQ=255)

                                                                         cGTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGGGATGACCGTCAtt               >  2:3447160/1‑93 (MQ=255)

                                                                          gTTTATTGATTTTTCGGCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGGCAGACTCATCTTCAGGCTGGAGGGGATGGCTGTCCttt              >  2:3905989/1‑93 (MQ=255)

                                                                          gTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAttt              >  1:3002445/1‑93 (MQ=255)

                                                                          gTTTATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAttt              >  2:3649548/1‑93 (MQ=255)

                                                                            ttATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCAtttac            >  2:4669950/1‑91 (MQ=255)

                                                                            ttATTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTcc            >  1:3563134/1‑93 (MQ=255)

                                                                              aTTGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCt          >  2:1813937/1‑93 (MQ=255)

                                                                               ttGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTTCCTg         >  1:1950629/1‑93 (MQ=255)

                                                                               ttGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTg         >  2:1988117/1‑93 (MQ=255)

                                                                                tGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGa        <  1:1068177/93‑1 (MQ=255)

                                                                                tGATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGa        >  1:4437820/1‑93 (MQ=255)

                                                                                 gATTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGAt       >  2:3888225/1‑93 (MQ=255)

                                                                                  aTTTTTCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATg      >  2:1933518/1‑93 (MQ=255)

                                                                                   tttttCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGc     <  2:2955625/93‑1 (MQ=255)

                                                                                   tttttCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGc     <  1:4023436/93‑1 (MQ=255)

                                                                                    ttttCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCg    <  2:359696/93‑1 (MQ=255)

                                                                                      ttCGCCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGag  >  1:4320698/1‑93 (MQ=255)

GGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGATGCACAGAAAGATCGCTTTCATCGACGGTCAGACTCATCATCAGGCTGGAGGTGATGACCGTCATTTCCCTGATGCGAG  >  NC_000913/1344434‑1344612

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