BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,972,641	G→A	R145H (CGC→CAC)	wecD →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,972,641	0	G	A	99.7%	1361.7 / ‑6.0	377	R145H (CGC→CAC)	wecD	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (1/0); new base A (162/214); total (163/214)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 384 total aligned reads displayed.

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAAG  >  NC_000913/3972553‑3972728

aaTGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGTTCACGGt                                                                                     <  1:2809079/93‑1 (MQ=255)

    cGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTaaaa                                                                                 >  2:4494388/1‑93 (MQ=255)

       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                              >  2:6176753/1‑93 (MQ=255)

       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                              >  2:3816470/1‑93 (MQ=255)

       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                              <  2:9566160/93‑1 (MQ=255)

       cAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAaacaa                                                                              <  1:8265267/93‑1 (MQ=255)

         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                            >  1:5881282/1‑93 (MQ=255)

         gATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACt                                                                            <  1:3433519/93‑1 (MQ=255)

             ggcgAATTGGCCTGCGGGCTGGACGCGGGGCAGGTGCTGAGCTGATGCAAACGTCGCTAAACTGGGCGGATGCTCACGGTAAAACAACCTTGc                                                                        <  2:9092143/92‑1 (MQ=255)

             cgcgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGc                                                                        >  2:2646370/1‑93 (MQ=255)

               cgAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCgg                                                                      >  1:4085398/1‑93 (MQ=255)

                gAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCggg                                                                     <  1:2020143/93‑1 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                    <  1:2307634/93‑1 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                    >  1:6759537/1‑93 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                    <  2:3491068/93‑1 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                    <  1:8677048/93‑1 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                    <  2:6405935/93‑1 (MQ=255)

                 aaTTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGt                                                                    <  2:8601217/93‑1 (MQ=255)

                  aTTGGCCTGATGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCATGGTAAAACAACTTTGCGGGTg                                                                   <  2:2032175/93‑1 (MQ=255)

                   ttGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                                  >  1:4314791/1‑93 (MQ=255)

                   atGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTgg                                                                  <  1:5168411/92‑1 (MQ=255)

                    tGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGc                                                                 <  1:9458622/93‑1 (MQ=255)

                      gCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGa                                                               >  2:501718/1‑93 (MQ=255)

                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                            >  1:1099395/1‑93 (MQ=255)

                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                            >  2:7202754/1‑93 (MQ=255)

                         tgctgGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGAccc                                                            <  2:5644296/93‑1 (MQ=255)

                          gctggctgGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCa                                                           >  2:1320859/1‑93 (MQ=255)

                           ctggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaa                                                          <  1:6042506/93‑1 (MQ=255)

                            tggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaa                                                         >  1:7702242/1‑93 (MQ=255)

                             ggctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCaaag                                                        >  1:5398729/1‑92 (MQ=255)

                              gctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATg                                                       <  2:7376704/93‑1 (MQ=255)

                               ctggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATgg                                                      <  1:9983163/93‑1 (MQ=255)

                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                    >  2:7552687/1‑93 (MQ=255)

                                 ggACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGc                                                    <  2:7996418/93‑1 (MQ=255)

                                  gACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCa                                                   <  2:8109669/93‑1 (MQ=255)

                                   aCGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCaa                                                  <  1:9711680/93‑1 (MQ=255)

                                     gcgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAaca                                                <  1:3713221/93‑1 (MQ=255)

                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                               <  2:1433628/93‑1 (MQ=255)

                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                               <  1:1015665/93‑1 (MQ=255)

                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                               <  2:7006067/93‑1 (MQ=255)

                                      cgGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAacac                                               <  1:6633740/93‑1 (MQ=255)

                                        gTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccg                                             >  2:2824743/1‑93 (MQ=255)

                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                            >  1:9540395/1‑93 (MQ=255)

                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                            <  1:8163824/93‑1 (MQ=255)

                                         tGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgc                                            >  1:4873163/1‑93 (MQ=255)

                                           cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                          >  1:6899561/1‑93 (MQ=255)

                                           cAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACAccgccg                                          <  1:7589571/93‑1 (MQ=255)

                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                         <  1:8564495/93‑1 (MQ=255)

                                            aGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACcgccgc                                         <  2:4338873/93‑1 (MQ=255)

                                             ggTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCcgcg                                        >  2:2123270/1‑93 (MQ=255)

                                               tGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCt                                      >  2:6591719/1‑93 (MQ=255)

                                                gctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCtt                                     <  1:5564635/93‑1 (MQ=255)

                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                    <  1:4624254/93‑1 (MQ=255)

                                                 ctgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTa                                    >  1:3978788/1‑93 (MQ=255)

                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                                   <  1:4500574/93‑1 (MQ=255)

                                                  tgagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaa                                   >  2:7470479/1‑93 (MQ=255)

                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCGAAATGGGCAACACGGCCGCGCTTaaa                                  >  1:3118178/1‑93 (MQ=255)

                                                   gagctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTaaa                                  <  1:5010076/93‑1 (MQ=255)

                                                    agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                                 <  1:2249858/93‑1 (MQ=255)

                                                    agctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAAc                                 <  1:6483899/93‑1 (MQ=255)

                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                                >  2:4588472/1‑93 (MQ=255)

                                                     gctgaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGCCGACCCAAATGGGCAACACCGCCGCGCTTAAACg                                >  1:1612143/1‑93 (MQ=255)

                                                        gaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                             >  2:1885376/1‑93 (MQ=255)

                                                        gaTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGata                             <  1:7486312/93‑1 (MQ=255)

                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                            <  1:7920544/93‑1 (MQ=255)

                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                            <  1:1692789/93‑1 (MQ=255)

                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                            <  1:8015430/93‑1 (MQ=255)

                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                            <  2:6758900/93‑1 (MQ=255)

                                                         aTGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatac                            <  1:7282756/93‑1 (MQ=255)

                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                           >  2:9356646/1‑93 (MQ=255)

                                                          tGCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGataca                           <  1:7090865/93‑1 (MQ=255)

                                                           gCAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacat                          >  2:412672/1‑93 (MQ=255)

                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                         >  1:1398637/1‑93 (MQ=255)

                                                            cAAACGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacata                         <  1:7857531/93‑1 (MQ=255)

                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                        <  1:9269690/93‑1 (MQ=255)

                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                        >  2:5457642/1‑93 (MQ=255)

                                                             aaaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGatacatac                        <  2:2439375/93‑1 (MQ=255)

                                                              aaCGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGAtacataca                       >  1:6320514/1‑93 (MQ=255)

                                                                cGGCGCTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACaaa                     <  2:7344369/93‑1 (MQ=255)

                                                                   cgcTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtg                  <  2:9855078/93‑1 (MQ=255)

                                                                    gcTAAACTGGGCGTATGCTCACGGTGAAACAACTTTGGGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGGTACATACAAAgcgg                 >  2:5893560/1‑93 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt                >  1:9030937/1‑93 (MQ=255)

                                                                     cTAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggt                <  2:3391234/93‑1 (MQ=255)

                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAGTGGGCAACACCGCCGTGCTTAAACGATACATACAAAgtggtg               <  1:5279465/93‑1 (MQ=255)

                                                                      tAAACTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAgtggtg               >  2:884585/1‑93 (MQ=255)

                                                                       aaaCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGc              >  2:2303704/1‑93 (MQ=255)

                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCGGGGCTTAAAGGAAACATACAAAGGGGGGCGa            >  1:9844606/1‑93 (MQ=255)

                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa            >  1:9606063/1‑93 (MQ=255)

                                                                         aCTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGAAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGa            >  1:1194782/1‑93 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa           <  1:5819252/93‑1 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa           <  1:1728783/93‑1 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa           >  2:6258713/1‑93 (MQ=255)

                                                                          cTGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGaa           <  2:9675617/93‑1 (MQ=255)

                                                                           tGGGCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAAt          <  1:9809958/93‑1 (MQ=255)

                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt        <  2:1753633/93‑1 (MQ=255)

                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGt        >  2:4724041/1‑93 (MQ=255)

                                                                             ggCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCACTTAAACGATACATACAAAGTGGCGCGTAtct        >  2:8433519/1‑91 (MQ=255)

                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa       <  1:5997104/93‑1 (MQ=255)

                                                                              gCGTATGCTCACGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTa       <  2:92620/93‑1 (MQ=255)

                                                                                   tgctcgCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAAg  >  1:8236614/4‑93 (MQ=255)

AATGCGACAGATGCGCGAATTGGCCTGCTGGCTGGACGCGGTGCAGGTGCTGAGCTGATGCAAACGGCGCTAAACTGGGCGTATCGTCGCGGTAAAACAACTTTGCGGGTGGCGACCCAAATGGGCAACACCGCCGCGCTTAAACGATACATACAAAGTGGTGCGAATGTAGAAAG  >  NC_000913/3972553‑3972728

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