BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,052,844	T→C	I5M (ATA→ATG)	glnG ←	GlnG

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,052,844	0	T	C	98.1%	1142.4 / 3.3	376	I5M (ATA→ATG)	glnG	GlnG
Reads supporting (aligned to +/- strand): ref base T (1/6); new base C (167/202); total (168/208)
Fisher's exact test for biased strand distribution p-value = 1.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 388 total aligned reads displayed.

GTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATT  >  NC_000913/4052754‑4052933

gTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCAt                                                                                         >  1:6452515/1‑93 (MQ=255)

gTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCAt                                                                                         <  1:2809978/93‑1 (MQ=255)

 ttCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATc                                                                                        >  1:1855924/1‑93 (MQ=255)

   ctcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATccc                                                                                      >  1:315480/1‑93 (MQ=255)

   ctcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATccc                                                                                      >  2:7933707/1‑93 (MQ=255)

    tcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCt                                                                                     >  1:7056418/1‑93 (MQ=255)

    tcAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCt                                                                                     >  1:4874216/1‑93 (MQ=255)

     cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTc                                                                                    <  1:5550914/93‑1 (MQ=255)

     cAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTc                                                                                    >  2:2220337/1‑93 (MQ=255)

      aaaCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCg                                                                                   >  1:5075698/1‑93 (MQ=255)

        aCGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGtt                                                                                 >  1:5010199/1‑93 (MQ=255)

          gtcgtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTga                                                                               >  2:6066188/1‑92 (MQ=255)

          gtcgtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGc                                                                               >  2:3926962/1‑93 (MQ=255)

          gtcgtACAGATTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGc                                                                               <  2:7565187/93‑1 (MQ=255)

           tcgtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCa                                                                              >  1:2105923/1‑93 (MQ=255)

             gtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATa                                                                            >  1:6007071/1‑93 (MQ=255)

             gtACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCAACATCGACTACCCAGCCCATCCCTCGATGCATa                                                                            >  2:4040528/1‑93 (MQ=255)

              tACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaa                                                                           >  1:7204032/1‑93 (MQ=255)

               aCAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaaa                                                                          <  2:82811/93‑1 (MQ=255)

               aCAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATaaa                                                                          <  1:3863772/93‑1 (MQ=255)

                cAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAAc                                                                         <  2:7458559/93‑1 (MQ=255)

                cAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTAATCGACTACCCAGACCATCCGTCGTTGCATaaag                                                                         >  1:5428422/1‑92 (MQ=255)

                  ggTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGt                                                                       >  2:1542244/1‑93 (MQ=255)

                    ttAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCa                                                                     <  2:7646393/93‑1 (MQ=255)

                    ttAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCa                                                                     <  2:6851008/93‑1 (MQ=255)

                      aaaCCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCAcc                                                                   <  1:5271061/93‑1 (MQ=255)

                          cTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTa                                                               <  1:6089576/93‑1 (MQ=255)

                          cTGCCCAAGAGAGCGCACGTTCAAGCACCCACCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTa                                                               <  2:5419669/93‑1 (MQ=255)

                           tGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACTCAGACCATCCCTCGTTGCATAAACGTCACCTTTAt                                                              >  2:26312/1‑93 (MQ=255)

                           tGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTAt                                                              >  2:4629769/1‑93 (MQ=255)

                            ggcccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTAtt                                                             <  1:1230902/91‑1 (MQ=255)

                               ccAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTcc                                                          <  1:2943703/93‑1 (MQ=255)

                                cAGCGAGCGCACGTTCAAGCACCCAACGGATGGTACTATCGTCATCGACTACCCAGACAATCCCTCGTTGCATAAACGTCACCTTTATTTCCt                                                         >  2:6749419/1‑93 (MQ=255)

                                cAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCt                                                         <  2:431609/93‑1 (MQ=255)

                                 agcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTg                                                        <  1:7324712/93‑1 (MQ=255)

                                 agcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTg                                                        >  2:1355481/1‑93 (MQ=255)

                                  gcgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGa                                                       >  2:6472697/1‑93 (MQ=255)

                                   cgagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGAt                                                      <  1:1051156/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  1:5212090/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  1:4197009/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  1:5883866/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  1:951756/93‑1 (MQ=255)

                                    gagcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATa                                                     <  2:4824921/93‑1 (MQ=255)

                                     agcgCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAg                                                    <  1:7859596/93‑1 (MQ=255)

                                        gcACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCa                                                 <  1:1172934/93‑1 (MQ=255)

                                         cACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAg                                                <  1:3452483/93‑1 (MQ=255)

                                          aCGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAgg                                               <  2:1026924/93‑1 (MQ=255)

                                           cGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGt                                              >  2:7279606/1‑93 (MQ=255)

                                            gTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTa                                             >  2:1743009/1‑93 (MQ=255)

                                            gTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTa                                             <  2:1424520/93‑1 (MQ=255)

                                             ttCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTaa                                            >  2:5469618/1‑93 (MQ=255)

                                              tCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTaaa                                           >  2:4233554/1‑93 (MQ=255)

                                               cAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCCCCTTTATTCCCTGATAGGCAGGTAAAc                                          <  1:847530/93‑1 (MQ=255)

                                               cAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAc                                          <  1:4015706/93‑1 (MQ=255)

                                                aaGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAAcc                                         >  1:3698305/1‑93 (MQ=255)

                                                 aGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCg                                        <  2:7372537/93‑1 (MQ=255)

                                                   cACCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgag                                      <  2:5603714/93‑1 (MQ=255)

                                                    aCCCAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCgaga                                     >  1:2173207/1‑93 (MQ=255)

                                                     cccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGaa                                    >  1:1931029/1‑93 (MQ=255)

                                                      ccAACGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAAc                                   <  1:3544837/93‑1 (MQ=255)

                                                        agcGGATGGAACGATCGTCATCGACTGCCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTc                                 <  1:7559642/91‑1 (MQ=255)

                                                        aaCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTc                                 <  1:2412516/93‑1 (MQ=255)

                                                        aaCGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTc                                 >  1:1288753/1‑93 (MQ=255)

                                                          cGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCgg                               >  2:2884099/1‑93 (MQ=255)

                                                          cGGATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCgg                               >  1:3277002/1‑93 (MQ=255)

                                                           ggATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGt                              <  2:4544/93‑1 (MQ=255)

                                                            gATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                             <  2:7792584/93‑1 (MQ=255)

                                                            gATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                             <  1:4561690/93‑1 (MQ=255)

                                                            gATGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTa                             <  2:4942824/93‑1 (MQ=255)

                                                             aTGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTAt                            >  2:6596556/1‑93 (MQ=255)

                                                             aTGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTAt                            <  1:5309984/93‑1 (MQ=255)

                                                             aTGGAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTAt                            >  1:4736974/1‑93 (MQ=255)

                                                                gAACTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGcc                         >  1:249719/1‑93 (MQ=255)

                                                                 aaCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGccc                        <  1:6954260/93‑1 (MQ=255)

                                                                 aaCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGccc                        <  1:1747497/93‑1 (MQ=255)

                                                                 aaCTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGccc                        >  2:5256338/1‑93 (MQ=255)

                                                                   cTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTg                      <  1:6637293/93‑1 (MQ=255)

                                                                   cTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTg                      <  1:6393229/93‑1 (MQ=255)

                                                                   cTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTg                      <  2:3658195/93‑1 (MQ=255)

                                                                   cTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTg                      <  1:621955/93‑1 (MQ=255)

                                                                   cTATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGATATGCCCTg                      >  2:4363481/1‑93 (MQ=255)

                                                                    tATCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTgg                     >  2:7195346/1‑93 (MQ=255)

                                                                     aTCGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGc                    >  2:6950184/1‑93 (MQ=255)

                                                                       cGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCa                  >  2:1925492/1‑93 (MQ=255)

                                                                       cGTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCa                  >  2:328816/1‑93 (MQ=255)

                                                                        gTCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCaa                 <  1:2253586/93‑1 (MQ=255)

                                                                         tCATCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAAc                <  2:3087326/93‑1 (MQ=255)

                                                                           aTCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTg              >  1:7748534/1‑93 (MQ=255)

                                                                           aTCGACTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTg              <  1:746699/93‑1 (MQ=255)

                                                                               aCTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTCCCTGATAGTCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaa          <  2:243756/93‑1 (MQ=255)

                                                                               aCTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTaa          <  1:819151/93‑1 (MQ=255)

                                                                               aCTACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGACAGGCAGGTAAACCGAGAACTCGGTATGCCGTGGCCAACTGGGaa          >  1:4265115/1‑93 (MQ=255)

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt        >  2:3470585/1‑93 (MQ=255)

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt        <  2:5770478/93‑1 (MQ=255)

                                                                                 tACCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAt        >  1:4498521/1‑93 (MQ=255)

                                                                                  aCCCAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAAtt       <  1:1555795/93‑1 (MQ=255)

                                                                                   cccAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGACCTCGGTATGCCCTGGCCAACTGGTAAATTc      <  2:7091411/93‑1 (MQ=255)

                                                                                    ccAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCa     <  2:3576651/93‑1 (MQ=255)

                                                                                     cAGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCaa    <  1:3627522/93‑1 (MQ=255)

                                                                                      aGACCATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAAt   <  1:5651914/93‑1 (MQ=255)

GTTCTCAAACGTCGTACAGGTTAAACCTGCCCCAGCGAGCGCACGTTCAAGCACCCAACGGATGGAACTATCGTCATCGACTACCCAGACTATCCCTCGTTGCATAAACGTCACCTTTATTTCCTGATAGGCAGGTAAACCGAGAACTCGGTATGCCCTGGCCAACTGGTAAATTCAATT  >  NC_000913/4052754‑4052933

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