BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	394,595	T→C	D306G (GAT→GGT)	yaiH ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	394,595	0	T	C	99.5%	750.9 / ‑5.1	219	D306G (GAT→GGT)	yaiH	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (107/111); total (108/111)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 230 total aligned reads displayed.

GTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTG  >  NC_000913/394505‑394684

  aaTGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGc                                                                                       <  1:1888311/93‑1 (MQ=255)

  aaTGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGc                                                                                       <  1:4127259/93‑1 (MQ=255)

   aTGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCt                                                                                      >  2:3296661/1‑93 (MQ=255)

    tGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCtt                                                                                     >  2:2472160/1‑93 (MQ=255)

       aaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGc                                                                                  >  1:5411350/1‑93 (MQ=255)

       aaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGc                                                                                  >  1:3737684/1‑93 (MQ=255)

       aaCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACACGCTTTGc                                                                                  <  1:6304135/93‑1 (MQ=255)

        aCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGcg                                                                                 >  2:2113953/1‑92 (MQ=255)

        aCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGcc                                                                                 <  1:6190637/93‑1 (MQ=255)

        aCCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGcc                                                                                 <  1:5032481/93‑1 (MQ=255)

          cGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACTCGCTTTGCCgg                                                                               >  1:3266363/1‑93 (MQ=255)

          cGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgg                                                                               <  1:1104527/93‑1 (MQ=255)

          cGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgg                                                                               <  2:4458313/93‑1 (MQ=255)

          cGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgg                                                                               >  2:134651/1‑93 (MQ=255)

          cGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCgg                                                                               >  1:435327/1‑93 (MQ=255)

           gCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCggg                                                                              >  2:6034930/1‑93 (MQ=255)

           gCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCggg                                                                              >  2:5712483/1‑93 (MQ=255)

             caccacCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGa                                                                            >  2:753004/1‑93 (MQ=255)

             caccacCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGa                                                                            <  1:2856041/93‑1 (MQ=255)

                 accGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATc                                                                        >  2:2650190/1‑93 (MQ=255)

                  ccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATcc                                                                       <  2:586369/93‑1 (MQ=255)

                  ccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATcc                                                                       >  2:5379809/1‑93 (MQ=255)

                  ccGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATcc                                                                       <  2:4688624/93‑1 (MQ=255)

                    gCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCAt                                                                     <  2:2733747/93‑1 (MQ=255)

                    gCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCAt                                                                     <  2:3228924/93‑1 (MQ=255)

                    gCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCAt                                                                     >  2:3513626/1‑93 (MQ=255)

                     ccTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATg                                                                    >  2:3707025/1‑93 (MQ=255)

                     ccTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATg                                                                    >  1:4870113/1‑93 (MQ=255)

                      cTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGc                                                                   <  1:4941393/93‑1 (MQ=255)

                       tGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGcc                                                                  <  2:5166970/93‑1 (MQ=255)

                       tGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGcc                                                                  <  2:6158462/93‑1 (MQ=255)

                       tGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGcc                                                                  >  2:2153765/1‑93 (MQ=255)

                        gTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCa                                                                 >  1:144462/1‑93 (MQ=255)

                          cttctGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAAt                                                               >  1:4528730/1‑93 (MQ=255)

                          cttctGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAAt                                                               >  1:890798/1‑93 (MQ=255)

                           ttctGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATc                                                              <  1:1607204/93‑1 (MQ=255)

                              tGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACt                                                           >  2:846782/1‑93 (MQ=255)

                              tGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACt                                                           >  1:3385163/1‑93 (MQ=255)

                               gAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACtt                                                          <  2:5947351/93‑1 (MQ=255)

                               gAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACtt                                                          >  2:3467331/1‑93 (MQ=255)

                                aataacCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACttt                                                         <  1:2840939/93‑1 (MQ=255)

                                   aatCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGt                                                      >  2:5550477/1‑93 (MQ=255)

                                   aatCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGt                                                      <  1:4123555/93‑1 (MQ=255)

                                   aatCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGt                                                      >  1:5205204/1‑93 (MQ=255)

                                       cccGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAAc                                                  <  1:589347/93‑1 (MQ=255)

                                        ccGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACt                                                 >  2:682923/1‑93 (MQ=255)

                                        ccGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACt                                                 <  2:3577295/93‑1 (MQ=255)

                                         cGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTg                                                >  2:1562357/1‑93 (MQ=255)

                                         cGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTg                                                >  2:1777896/1‑93 (MQ=255)

                                          ggACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgc                                               <  2:3845027/93‑1 (MQ=255)

                                          ggACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgc                                               <  1:1043051/93‑1 (MQ=255)

                                          agacgacCCTCTATGGGGGCCATATATCCCCAACCTAAACCAAGCGCGCCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgc                                               <  1:5776282/92‑1 (MQ=255)

                                           gacgacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcg                                              <  2:489268/93‑1 (MQ=255)

                                           gacgacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcg                                              >  1:2084970/1‑93 (MQ=255)

                                            acgacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgc                                             >  1:2989554/1‑93 (MQ=255)

                                            acgacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgc                                             <  1:2554802/93‑1 (MQ=255)

                                            acgacCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgc                                             >  1:643611/1‑93 (MQ=255)

                                               acCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcg                                          >  2:1918741/1‑93 (MQ=255)

                                               acCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcg                                          >  1:1725030/1‑93 (MQ=255)

                                               acCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcg                                          >  2:2528753/1‑93 (MQ=255)

                                               acCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcg                                          >  1:6109135/1‑93 (MQ=255)

                                                cccTCTTTCGGGGCCATATATACCCAACCTAACCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcgc                                         <  1:2350445/93‑1 (MQ=255)

                                                cccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcgc                                         <  2:4482280/93‑1 (MQ=255)

                                                cccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTgcgcgcgc                                         <  1:999569/93‑1 (MQ=255)

                                                 ccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCt                                        >  2:1726691/1‑93 (MQ=255)

                                                 ccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCt                                        >  1:4558064/1‑93 (MQ=255)

                                                 ccTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCt                                        >  1:4533480/1‑93 (MQ=255)

                                                    ctTTCGGGGCCATATATACCCACCCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGAt                                     <  1:2670422/93‑1 (MQ=255)

                                                     tttCGGGGCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATa                                    >  2:3119764/1‑93 (MQ=255)

                                                         ggggCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATc                                >  2:4760372/1‑93 (MQ=255)

                                                          gggCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCa                               >  1:545033/1‑93 (MQ=255)

                                                           ggCCATATATACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATTaa                              <  2:1627601/93‑1 (MQ=255)

                                                            gCCATATATCCCCACCCTAACCCAAGCGCACCCGCTATGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAAc                             <  1:2259491/93‑1 (MQ=255)

                                                            gCCATATATACCCAACCTAACCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAAc                             <  1:3637415/93‑1 (MQ=255)

                                                               atatatACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTc                          >  2:4606593/1‑93 (MQ=255)

                                                               atatatACCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTc                          <  2:4416982/93‑1 (MQ=255)

                                                                tatataCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCt                         <  2:3068946/93‑1 (MQ=255)

                                                                tatataCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCt                         <  1:5321043/93‑1 (MQ=255)

                                                                 atataCCCAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTg                        <  2:2939080/93‑1 (MQ=255)

                                                                     aCCCAACCTAAACCAAGAGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTACGCGCGCTGATAAATCAACGTCTGCATg                    >  1:3062428/1‑93 (MQ=255)

                                                                      cccAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGc                   >  2:1990935/1‑93 (MQ=255)

                                                                        cAACCTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGa                 >  2:3641395/1‑93 (MQ=255)

                                                                           ccTGAACCAAGCGCACCCGCTTTGCGCGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCt              <  2:6368041/93‑1 (MQ=255)

                                                                           ccTAAACCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCt              >  2:5039311/1‑93 (MQ=255)

                                                                              aaaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGcc           <  2:4064795/93‑1 (MQ=255)

                                                                              aaaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGcc           <  1:117425/93‑1 (MQ=255)

                                                                               aaCCAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCt          <  2:4407039/93‑1 (MQ=255)

                                                                                 ccAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTgg        >  1:5099895/1‑93 (MQ=255)

                                                                                  cAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGt       <  2:1135673/93‑1 (MQ=255)

                                                                                  cAAGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGt       >  1:4648008/1‑93 (MQ=255)

                                                                                   aaGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGtt      >  2:433531/1‑93 (MQ=255)

                                                                                    aGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTa     >  1:5640395/1‑93 (MQ=255)

                                                                                    aGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTa     <  2:2171404/93‑1 (MQ=255)

                                                                                    aGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTa     >  2:5747960/1‑93 (MQ=255)

                                                                                    aGCGCACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTa     >  1:4282778/1‑93 (MQ=255)

                                                                                     gcgcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTAc    <  1:4813332/93‑1 (MQ=255)

                                                                                     gcgcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTAc    <  1:6043265/93‑1 (MQ=255)

                                                                                     gcgcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTAc    <  2:3359571/93‑1 (MQ=255)

                                                                                       gcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTg  <  1:4048644/93‑1 (MQ=255)

                                                                                       gcACCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTg  <  1:4932907/93‑1 (MQ=255)

GTAATGAAACCGCCACCACCGCCTGTCTTCTGAATAATCCCCGGACGACCCTCTTTCGGGGCCATATATACCCAACCTAAACCAAGCGCATCCGCTTTGCCGGGGACATCCATGCCAATCACTTTCGTAAACTGCGCGCGCTGATAAATCAACGTCTGCATGCGATCTGCCTGGTTACTG  >  NC_000913/394505‑394684

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