BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,355	A→G	H314R (CAC→CGC)	metL →	aspartokinase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,355	0	A	G	99.6%	779.8 / ‑4.8	242	H314R (CAC→CGC)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (113/128); total (114/128)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 261 total aligned reads displayed.

AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444

aaCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCcgcg                                                                                         >  1:3168673/1‑93 (MQ=255)

aaCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCcgcg                                                                                         >  2:2804148/1‑93 (MQ=255)

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        <  2:4145549/93‑1 (MQ=255)

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        <  2:3606796/93‑1 (MQ=255)

  cTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgat                                                                                       >  1:5782078/1‑93 (MQ=255)

     ggcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatgat                                                                                    <  1:970052/92‑1 (MQ=255)

       cTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGt                                                                                  >  2:5752658/1‑93 (MQ=255)

        tGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTc                                                                                 >  2:5471810/1‑93 (MQ=255)

        tGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTc                                                                                 >  2:1831252/1‑93 (MQ=255)

             cTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGttt                                                                            <  2:155468/93‑1 (MQ=255)

             cTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGttt                                                                            >  1:3341515/1‑93 (MQ=255)

             cTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGttt                                                                            <  2:4687589/93‑1 (MQ=255)

              tACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttg                                                                           >  2:6267208/1‑93 (MQ=255)

              tACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttg                                                                           <  2:3365623/93‑1 (MQ=255)

               acacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgt                                                                          >  1:4137595/1‑92 (MQ=255)

                cacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgat                                                                         >  2:647951/1‑93 (MQ=255)

                cacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgat                                                                         <  2:275362/93‑1 (MQ=255)

                cacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgat                                                                         <  1:5924678/93‑1 (MQ=255)

                  cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattg                                                                       >  2:6570186/1‑93 (MQ=255)

                  cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattg                                                                       <  2:2395843/93‑1 (MQ=255)

                  cGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattg                                                                       >  2:1640575/1‑93 (MQ=255)

                   gCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattga                                                                      <  2:6414884/93‑1 (MQ=255)

                    ccGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAg                                                                     <  1:2467602/93‑1 (MQ=255)

                    ccGGATAAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAg                                                                     <  1:5594769/93‑1 (MQ=255)

                     cGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGt                                                                    <  2:903835/93‑1 (MQ=255)

                     cGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGt                                                                    <  1:4398214/93‑1 (MQ=255)

                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   >  1:1620679/1‑93 (MQ=255)

                      ggATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGtt                                                                   <  1:2830902/93‑1 (MQ=255)

                       gATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGttt                                                                  <  2:4480567/93‑1 (MQ=255)

                       gATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGttt                                                                  <  2:2553913/93‑1 (MQ=255)

                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTc                                                                 <  2:4296525/93‑1 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                <  2:5410020/93‑1 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                <  1:4447117/93‑1 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                <  2:4955235/93‑1 (MQ=255)

                           aaGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAgg                                                              >  1:2672078/1‑93 (MQ=255)

                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAggg                                                             >  1:1271052/1‑92 (MQ=255)

                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             <  1:2988580/93‑1 (MQ=255)

                            aGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             >  1:2856707/1‑93 (MQ=255)

                             ggTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTg                                                            <  1:2613722/93‑1 (MQ=255)

                             ggTGCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTg                                                            <  2:6048479/93‑1 (MQ=255)

                              gTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGc                                                           >  2:3210181/1‑93 (MQ=255)

                                  cACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTgcgcgc                                                       >  1:504189/1‑93 (MQ=255)

                                  cACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGc                                                       >  2:3096979/1‑93 (MQ=255)

                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      <  1:6289592/93‑1 (MQ=255)

                                    cgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     >  1:5188237/1‑93 (MQ=255)

                                     gcgcATTGCACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGCTTGAGTTTCAGGGGCCCGCCac                                                    >  1:1833325/1‑92 (MQ=255)

                                     gcgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAg                                                    >  2:793182/1‑93 (MQ=255)

                                     gcgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAg                                                    >  2:111101/1‑93 (MQ=255)

                                      cgcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGt                                                   <  2:3008607/93‑1 (MQ=255)

                                      cccATTAAACGCTTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTCCCCGCCAGt                                                   <  1:6661643/91‑1 (MQ=255)

                                       gcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTc                                                  >  1:4726099/1‑93 (MQ=255)

                                       gcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAATTTCAGGTGCCCGCCAGTc                                                  <  1:6541311/93‑1 (MQ=255)

                                       ccATTGAGCGCGTGCGGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTc                                                  <  1:428758/92‑1 (MQ=255)

                                        cATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCa                                                 <  1:512493/93‑1 (MQ=255)

                                         aTTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAg                                                >  1:5312508/1‑93 (MQ=255)

                                           tGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGa                                              >  1:159124/1‑93 (MQ=255)

                                            gAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAt                                             <  2:5669912/93‑1 (MQ=255)

                                            gAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAt                                             <  2:4820456/93‑1 (MQ=255)

                                              aCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAttt                                           <  1:2291067/93‑1 (MQ=255)

                                              aCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAttt                                           >  1:6227796/1‑93 (MQ=255)

                                              aCGCGTGCTGGCCTCCGGTACTGGTGCGCGTAGTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAttt                                           <  1:6470371/93‑1 (MQ=255)

                                                     cTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTg                                    <  2:2039976/93‑1 (MQ=255)

                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   >  2:3595423/1‑93 (MQ=255)

                                                      ggcccTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   <  1:4594357/90‑1 (MQ=255)

                                                       ggCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgc                                  >  2:1701611/1‑93 (MQ=255)

                                                        gccgccGGTACTGGTGCGCGTATAGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgcg                                 <  2:5694941/93‑1 (MQ=255)

                                                        gCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGgcg                                 <  1:3959313/93‑1 (MQ=255)

                                                         ccTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGACAGTCAGGATTTCAAACTGgcgc                                >  1:2260921/1‑93 (MQ=255)

                                                          cTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTTGCGCa                               >  1:746986/1‑93 (MQ=255)

                                                          cTCCGGTACTGGTGCGCGTACTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCa                               >  2:3349629/1‑93 (MQ=255)

                                                              ggTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATaaa                           <  1:1731472/93‑1 (MQ=255)

                                                               gTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAag                          >  1:106037/1‑93 (MQ=255)

                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         <  1:3846951/93‑1 (MQ=255)

                                                                  cTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAgaga                       >  2:1709096/1‑93 (MQ=255)

                                                                   tGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      <  1:4344130/93‑1 (MQ=255)

                                                                   tGGTGCGCGTATTCTCACCAGCCGCGATGATTTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGAATTTCAAACTGGCGCATAAAGAGAt                      <  1:4025925/93‑1 (MQ=255)

                                                                    ggTGCGCTTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATc                     <  1:6147941/93‑1 (MQ=255)

                                                                    ggTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATc                     <  1:3272075/93‑1 (MQ=255)

                                                                    ggTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTAGCGCATAAAGAGATc                     >  1:4896354/1‑93 (MQ=255)

                                                                    ggTGCGCGTATTATCTCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGCTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATc                     <  1:4094671/93‑1 (MQ=255)

                                                                      tGCGCGTATTGTCCCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGa                   <  1:5468803/93‑1 (MQ=255)

                                                                      tGCGCGTATTGCCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGa                   <  1:1252445/93‑1 (MQ=255)

                                                                      ggcgcGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGa                   <  1:2100525/92‑1 (MQ=255)

                                                                       gcgcGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGAc                  >  2:6656652/1‑93 (MQ=255)

                                                                       gcgcGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGAc                  <  1:3724415/93‑1 (MQ=255)

                                                                        cgcgTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGAcc                 <  1:3136188/93‑1 (MQ=255)

                                                                           gTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCaaa              <  2:278984/93‑1 (MQ=255)

                                                                           gTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCaaa              >  1:5954778/1‑93 (MQ=255)

                                                                           gTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCaaa              >  1:1064875/1‑93 (MQ=255)

                                                                            tATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             >  2:1344895/1‑93 (MQ=255)

                                                                              ttGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTGAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATcc           <  1:145847/93‑1 (MQ=255)

                                                                              ttGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATcc           >  2:4047482/1‑93 (MQ=255)

                                                                                gTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTg         >  2:5192963/1‑93 (MQ=255)

                                                                                  ccccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGaa       <  1:3743895/91‑1 (MQ=255)

                                                                                  cACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCGTAAAGAGATCGACCAAATCCTGaa       >  2:1228342/1‑93 (MQ=255)

                                                                                  cACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGaa       >  2:3888447/1‑93 (MQ=255)

                                                                                    ccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc     >  1:1292752/1‑93 (MQ=255)

                                                                                     cAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcg    <  1:1162812/93‑1 (MQ=255)

                                                                                       gCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAACTCCTGAAAcgcg  >  2:1955253/1‑93 (MQ=255)

                                                                                       gCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg  <  1:276247/93‑1 (MQ=255)

AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444

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