BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,162,448	A→G	T411A (ACA→GCA)	btuB →	BtuB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,162,448	0	A	G	84.4%	632.5 / 105.3	275	T411A (ACA→GCA)	btuB	BtuB
Reads supporting (aligned to +/- strand): ref base A (24/19); new base G (118/114); total (142/133)
Fisher's exact test for biased strand distribution p-value = 6.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.63e-01

Only 100 of 279 total aligned reads displayed.

ACGATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCG  >  NC_000913/4162356‑4162540

  gATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCa                                                                                            >  2:5644612/1‑93 (MQ=255)

   aTAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCAt                                                                                           >  1:1936699/1‑93 (MQ=255)

    tAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATc                                                                                          >  1:1001630/1‑93 (MQ=255)

    tAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATc                                                                                          <  1:2825738/93‑1 (MQ=255)

     aaCTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCt                                                                                         >  2:3261195/1‑93 (MQ=255)

        tCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTtat                                                                                      <  2:3719762/93‑1 (MQ=255)

        tCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTtat                                                                                      >  2:264359/1‑93 (MQ=255)

         cacaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTtata                                                                                     <  1:2923671/93‑1 (MQ=255)

          acaGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATaa                                                                                    >  2:6536786/1‑93 (MQ=255)

           caGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAg                                                                                   >  1:3170965/1‑93 (MQ=255)

            aGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAgg                                                                                  >  1:1477715/1‑93 (MQ=255)

            aGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAgg                                                                                  >  1:6677400/1‑93 (MQ=255)

             gTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGc                                                                                 >  1:3802564/1‑93 (MQ=255)

             gTTGGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCAAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGc                                                                                 >  1:6974107/1‑93 (MQ=255)

              tttGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCa                                                                                <  2:3419528/93‑1 (MQ=255)

              tttGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCa                                                                                <  2:547777/93‑1 (MQ=255)

               ttGGTCGTCATAGAGCCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCAc                                                                               <  2:2803808/93‑1 (MQ=255)

                tGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCAcc                                                                              >  1:2226185/1‑93 (MQ=255)

                 ggTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCa                                                                             >  2:4581520/1‑93 (MQ=255)

                      tcATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCt                                                                        <  1:959393/93‑1 (MQ=255)

                      tcATGGAACCTGGCAAACCAGCGCCGGTTGGAAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCt                                                                        >  1:6183547/1‑93 (MQ=255)

                        aTGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgg                                                                      <  1:4051537/93‑1 (MQ=255)

                        aTGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgg                                                                      <  2:173256/93‑1 (MQ=255)

                        aTGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgg                                                                      <  2:4936118/93‑1 (MQ=255)

                         tGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTggg                                                                     <  1:2429510/93‑1 (MQ=255)

                         tGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTggg                                                                     >  1:5421368/1‑93 (MQ=255)

                         tGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTggg                                                                     <  1:3522296/93‑1 (MQ=255)

                          ggAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTgggg                                                                    >  2:7363731/1‑93 (MQ=255)

                              ccTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAAc                                                                <  1:2891568/93‑1 (MQ=255)

                                tGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTg                                                              >  2:4773478/1‑93 (MQ=255)

                                tGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTg                                                              >  1:3024659/1‑93 (MQ=255)

                                 ggCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGt                                                             <  2:5875059/93‑1 (MQ=255)

                                  gCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTa                                                            >  2:7297288/1‑93 (MQ=255)

                                  gCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTa                                                            <  2:6746740/93‑1 (MQ=255)

                                  gCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTCATAAGGCACCAAATCTGGGGCAACTGTa                                                            >  2:2292421/1‑93 (MQ=255)

                                    aaaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATg                                                          <  2:7016280/93‑1 (MQ=255)

                                     aaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATgg                                                         >  1:3657506/1‑93 (MQ=255)

                                     aaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATgg                                                         <  2:7228557/93‑1 (MQ=255)

                                     aaCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATgg                                                         >  1:6514725/1‑93 (MQ=255)

                                      aCCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGc                                                        <  1:7305936/93‑1 (MQ=255)

                                            gcCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTAc                                                  >  1:1356163/1‑93 (MQ=255)

                                                gTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGaa                                              >  1:1810936/1‑93 (MQ=255)

                                                 ttGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGaaa                                             <  2:487347/93‑1 (MQ=255)

                                                   gggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATc                                           <  2:6318220/93‑1 (MQ=255)

                                                   gggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATc                                           >  1:4597858/1‑93 (MQ=255)

                                                   gggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATc                                           <  1:3416146/93‑1 (MQ=255)

                                                    ggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATcc                                          <  1:5119494/93‑1 (MQ=255)

                                                    ggAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATcc                                          <  1:4876765/93‑1 (MQ=255)

                                                     gAATTCATCGAAGGATATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCg                                         >  2:5180122/1‑93 (MQ=255)

                                                       aTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGaa                                       >  2:6176521/1‑93 (MQ=255)

                                                       aTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGaa                                       <  2:4110856/93‑1 (MQ=255)

                                                       aTTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGaa                                       <  1:5460360/93‑1 (MQ=255)

                                                            tCGAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTgg                                  >  2:2991072/1‑93 (MQ=255)

                                                             cGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGa                                 <  2:3949003/93‑1 (MQ=255)

                                                              gAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAc                                >  1:6391192/1‑93 (MQ=255)

                                                              gAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAc                                >  1:1245483/1‑93 (MQ=255)

                                                              gAAGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAc                                <  1:2725175/93‑1 (MQ=255)

                                                                aGGTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGAccc                              >  1:1165251/1‑93 (MQ=255)

                                                                 ggTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCg                             >  2:1094297/1‑93 (MQ=255)

                                                                 ggTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCg                             >  1:4719316/1‑93 (MQ=255)

                                                                  gTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCgg                            >  2:1863172/1‑93 (MQ=255)

                                                                  gTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCgg                            >  2:3888464/1‑93 (MQ=255)

                                                                  gTTATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCgg                            <  1:739522/93‑1 (MQ=255)

                                                                   ttATCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGga                           <  2:304241/93‑1 (MQ=255)

                                                                    tATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGgag                          >  1:5511506/1‑93 (MQ=255)

                                                                    tATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGgag                          >  1:937921/1‑93 (MQ=255)

                                                                      tCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGaa                        >  2:1166364/1‑93 (MQ=255)

                                                                      tCGCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGaa                        >  1:5684330/1‑93 (MQ=255)

                                                                        gCTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGaaaa                      >  1:2470859/1‑93 (MQ=255)

                                                                         cTTCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAg                     <  1:6047193/93‑1 (MQ=255)

                                                                         cTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAg                     <  2:3170965/93‑1 (MQ=255)

                                                                          ttCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGc                    >  2:7455658/1‑93 (MQ=255)

                                                                          ttCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGc                    <  1:3675913/93‑1 (MQ=255)

                                                                           tCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCa                   >  2:4263901/1‑93 (MQ=255)

                                                                           tCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCa                   >  2:4396506/1‑93 (MQ=255)

                                                                           tCATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCa                   >  1:516616/1‑93 (MQ=255)

                                                                            cATTGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCaa                  >  1:5802725/1‑93 (MQ=255)

                                                                               tGCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACa               <  2:5039376/93‑1 (MQ=255)

                                                                                gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg              <  2:1488545/93‑1 (MQ=255)

                                                                                gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg              <  1:5995163/93‑1 (MQ=255)

                                                                                gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg              >  1:5031624/1‑93 (MQ=255)

                                                                                gCTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAg              >  1:2003861/1‑93 (MQ=255)

                                                                                 cTTCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGt             <  1:848218/93‑1 (MQ=255)

                                                                                  ttCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTg            <  1:4973870/93‑1 (MQ=255)

                                                                                  ttCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTg            >  1:2783250/1‑93 (MQ=255)

                                                                                  ttCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACGGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAAccgcg            >  1:7165840/1‑89 (MQ=255)

                                                                                  ttCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTg            >  1:110022/1‑93 (MQ=255)

                                                                                   tCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTgg           <  2:3610049/93‑1 (MQ=255)

                                                                                   tCCTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTgg           <  1:7109480/93‑1 (MQ=255)

                                                                                     cTACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGa         >  1:996514/1‑93 (MQ=255)

                                                                                      tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCt                                    <  2:1388180/65‑1 (MQ=255)

                                                                                      tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCt                                    >  1:1388180/1‑65 (MQ=255)

                                                                                      tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa        <  2:6328941/93‑1 (MQ=255)

                                                                                      tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa        <  1:6960759/93‑1 (MQ=255)

                                                                                      tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa        >  1:6232186/1‑93 (MQ=255)

                                                                                      tACGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa        <  1:4501575/93‑1 (MQ=255)

                                                                                      tACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa        <  1:1605698/93‑1 (MQ=255)

                                                                                      tACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGaa        >  2:798548/1‑93 (MQ=255)

                                                                                        cGGGGCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAgg      <  2:4164618/93‑1 (MQ=255)

                                                                                         ggggCATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGgc     >  2:5715784/1‑93 (MQ=255)

ACGATAACTCACAGTTTGGTCGTCATGGAACCTGGCAAACCAGCGCCGGTTGGGAATTCATCGAAGGTTATCGCTTCATTGCTTCCTACGGGACATCTTATAAGGCACCAAATCTGGGGCAACTGTATGGCTTCTACGGAAATCCGAATCTGGACCCGGAGAAAAGCAAACAGTGGGAAGGCGCG  >  NC_000913/4162356‑4162540

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