BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	271,281	G→A	G76G (GGG→GGA)	b0257 →	putative transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	271,281	0	G	A	99.6%	830.8 / ‑5.4	230	G76G (GGG→GGA)	b0257	putative transposase
Reads supporting (aligned to +/- strand): ref base G (0/1); new base A (118/111); total (118/112)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

Only 100 of 237 total aligned reads displayed.

TGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGT  >  NC_000913/271191‑271373

tGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACt                                                                                            >  2:1373191/1‑93 (MQ=255)

 gCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTg                                                                                           >  1:4308878/1‑93 (MQ=255)

 gCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTg                                                                                           <  1:1541852/93‑1 (MQ=255)

 gCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCCCTGGCTTGCTGGCAGCCTCATGAAAGAGCTGGGACTg                                                                                           >  1:3471776/1‑93 (MQ=255)

    ggAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTc                                                                                        <  2:2215252/93‑1 (MQ=255)

       gCAAGAAGCATCGCCACAATGGCATCCCAGAGAGAATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGt                                                                                     <  2:6270066/93‑1 (MQ=255)

       gCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGt                                                                                     >  2:7819457/1‑93 (MQ=255)

       gCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGt                                                                                     <  2:3930251/93‑1 (MQ=255)

       gCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGt                                                                                     <  1:3654667/93‑1 (MQ=255)

       gCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGt                                                                                     <  1:3299296/93‑1 (MQ=255)

        cAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGtt                                                                                    >  2:1292508/1‑93 (MQ=255)

         aagaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTg                                                                                   >  1:3843161/1‑93 (MQ=255)

         aagaagCATCGCCACAATGACAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTg                                                                                   >  2:8083250/1‑93 (MQ=255)

          agaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGc                                                                                  <  1:136865/93‑1 (MQ=255)

           gaagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGcc                                                                                 <  2:8888340/93‑1 (MQ=255)

            aagCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCca                                                                                <  1:5210943/93‑1 (MQ=255)

              gCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCcagc                                                                              >  2:5560703/1‑93 (MQ=255)

                 tCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCagcagc                                                                           >  1:1606588/1‑93 (MQ=255)

                 tCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCagcagc                                                                           >  2:9562593/1‑93 (MQ=255)

                  cGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCGGTTGCCagcagcg                                                                          >  1:8431325/1‑92 (MQ=255)

                   gCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCg                                                                         >  1:171189/1‑93 (MQ=255)

                    ccACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGa                                                                        <  1:8890353/93‑1 (MQ=255)

                     cacaATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGAc                                                                       >  1:1016551/1‑93 (MQ=255)

                       caATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTc                                                                     >  1:9095809/1‑93 (MQ=255)

                       caAGGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTc                                                                     <  1:4154195/93‑1 (MQ=255)

                         aTGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCAc                                                                   >  1:5797014/1‑93 (MQ=255)

                          tGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCAcc                                                                  >  2:5400849/1‑93 (MQ=255)

                          tGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCAcc                                                                  <  2:7345820/93‑1 (MQ=255)

                           ggCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCg                                                                 <  2:4992015/93‑1 (MQ=255)

                             cAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGt                                                               <  1:2069137/93‑1 (MQ=255)

                             cAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGt                                                               <  2:931515/93‑1 (MQ=255)

                              aaCCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCggca                                                              >  1:7025799/1‑91 (MQ=255)

                               aCCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGtat                                                             <  1:938109/93‑1 (MQ=255)

                                 ccAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATaa                                                           <  1:106451/93‑1 (MQ=255)

                                 ccAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATaa                                                           <  1:2260716/93‑1 (MQ=255)

                                  cAGAGAGGATACCAGATGGGGCGCTGGCTAGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCAAGCAGCCAACTCACCGGTATAAg                                                          <  1:8940993/93‑1 (MQ=17)

                                   agagagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGc                                                         >  1:2481099/1‑93 (MQ=255)

                                     agagGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGt                                                       <  2:9617842/93‑1 (MQ=255)

                                           tACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGt                                                 >  1:3367006/1‑93 (MQ=255)

                                            aCCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTc                                                <  2:3485495/93‑1 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               >  1:959127/1‑93 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               <  2:334332/93‑1 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               >  1:64085/1‑93 (MQ=255)

                                             ccAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               <  1:3082656/93‑1 (MQ=255)

                                             ccAGATGGGGCGCTGCCTGGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTCCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               <  2:1543432/93‑1 (MQ=17)

                                             ccAGATGGGGAGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCa                                               <  1:1462105/93‑1 (MQ=255)

                                               aGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATg                                             >  2:9774886/1‑93 (MQ=255)

                                               aGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATg                                             <  2:1788652/93‑1 (MQ=255)

                                               aGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGGCATg                                             >  2:42641/1‑93 (MQ=255)

                                                 aTGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAg                                           >  1:2340236/1‑93 (MQ=255)

                                                 aTGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAg                                           >  2:6175667/1‑93 (MQ=255)

                                                   ggggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCa                                         <  2:9487967/93‑1 (MQ=255)

                                                   ggggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCa                                         >  2:4542788/1‑93 (MQ=255)

                                                   ggggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCa                                         <  1:7838992/93‑1 (MQ=255)

                                                    gggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCAc                                        >  1:7593698/1‑93 (MQ=255)

                                                     ggCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACg                                       >  2:4958322/1‑93 (MQ=255)

                                                        gcTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTg                                    <  2:100803/93‑1 (MQ=255)

                                                         cTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTgt                                   >  1:722997/1‑92 (MQ=255)

                                                          tGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCt                                  >  2:9822572/1‑93 (MQ=255)

                                                          tGGCTTGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCt                                  >  1:3777351/1‑93 (MQ=255)

                                                              ttGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATcc                              <  1:1202753/93‑1 (MQ=255)

                                                              ttGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATcc                              <  2:4897532/93‑1 (MQ=255)

                                                               tGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATccc                             <  2:2433580/93‑1 (MQ=255)

                                                               tGCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATccc                             >  2:2786674/1‑93 (MQ=255)

                                                                gCTGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCg                            >  1:5113604/1‑93 (MQ=255)

                                                                  tGGCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGaa                          <  2:1177733/93‑1 (MQ=255)

                                                                  gggaaGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGaa                          <  1:6698930/89‑1 (MQ=255)

                                                                   ggCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAat                         >  2:9155737/1‑93 (MQ=255)

                                                                   ggCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAat                         >  2:8559248/1‑93 (MQ=255)

                                                                    gCAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatc                        <  2:635335/93‑1 (MQ=255)

                                                                     cAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatca                       >  2:7477277/1‑93 (MQ=255)

                                                                     cAGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatca                       >  1:8352104/1‑93 (MQ=255)

                                                                      aGACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcat                      <  2:6146614/93‑1 (MQ=255)

                                                                       gACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcatc                     >  2:1674767/1‑93 (MQ=255)

                                                                       gACTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAatcatc                     >  2:5765596/1‑93 (MQ=255)

                                                                        aCTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCt                    <  2:4589594/93‑1 (MQ=255)

                                                                        aCTCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCt                    >  2:3521774/1‑93 (MQ=255)

                                                                          tCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTg                  >  1:2645449/1‑93 (MQ=255)

                                                                          tCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTg                  >  2:7318809/1‑93 (MQ=255)

                                                                          tCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCACGAATCATCTTg                  >  2:6534596/1‑93 (MQ=255)

                                                                          tCATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTAATGAGCACGTTGCTATCCCGAATCATCTTg                  <  2:8138133/93‑1 (MQ=255)

                                                                           cATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGa                 <  1:6540188/93‑1 (MQ=255)

                                                                           cATGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGa                 <  2:3365643/93‑1 (MQ=255)

                                                                            ctGAAAGGGCTGGTAGTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAg                <  2:3430991/92‑1 (MQ=21)

                                                                            aTGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAg                <  2:8363674/93‑1 (MQ=255)

                                                                             tGAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGc               >  2:7740655/1‑93 (MQ=255)

                                                                              gAAAGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCg              >  2:9899018/1‑93 (MQ=255)

                                                                              gAAAGAGCTGGGACTGGTCAGTTCCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCg              <  2:7217277/93‑1 (MQ=255)

                                                                               aaaGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGa             >  2:7009997/1‑93 (MQ=255)

                                                                               aaaGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGa             <  1:4543592/93‑1 (MQ=255)

                                                                                aaGAGCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGAc            <  1:1996551/93‑1 (MQ=255)

                                                                                 agagCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACa           <  1:9854286/93‑1 (MQ=255)

                                                                                 agagCTGGGACTGGTCAGTTGCCAGCAGCAGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACa           <  1:4887160/93‑1 (MQ=255)

                                                                                  gagCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAg          >  1:6147188/1‑93 (MQ=255)

                                                                                    gCTGGGACTGGTTAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGtt        <  2:3740113/93‑1 (MQ=255)

                                                                                    gCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGtt        <  1:8545535/93‑1 (MQ=255)

                                                                                    gCTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGtt        <  1:8153657/93‑1 (MQ=255)

                                                                                     cTGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTc       <  1:3231912/93‑1 (MQ=255)

                                                                                      tGGGACTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCg      >  2:6770500/1‑93 (MQ=255)

                                                                                       gggaCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGc     >  1:5923875/1‑93 (MQ=255)

TGCCGGAGCAAGAAGCATCGCCACAATGGCAACCCAGAGAGGATACCAGATGGGGCGCTGGCTTGCTGGCAGACTCATGAAAGAGCTGGGGCTGGTCAGTTGCCAGCAGCCGACTCACCGGTATAAGCGTGGCGGTCATGAGCACGTTGCTATCCCGAATCATCTTGAGCGACAGTTCGCCGT  >  NC_000913/271191‑271373

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