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	98.8%	1037.2 / ‑0.0	340	H314R (CAC→CGC)	metL	aspartokinase II
Reads supporting (aligned to +/- strand): ref base A (1/3); new base G (180/155); total (181/159)
Fisher's exact test for biased strand distribution p-value = 3.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.76e-01

Only 100 of 363 total aligned reads displayed.

AACTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAACGCG  >  NC_000913/4128265‑4128444

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        >  2:4014130/1‑93 (MQ=255)

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        >  1:1840856/1‑93 (MQ=255)

 aCTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCCGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCga                                                                                        >  1:1511977/1‑93 (MQ=255)

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

  cTGCGCTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgat                                                                                       <  2:2142688/93‑1 (MQ=255)

    gcgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatga                                                                                     >  2:1914473/1‑93 (MQ=255)

    gcgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatga                                                                                     <  1:3946833/93‑1 (MQ=255)

     cgcTGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCgatgat                                                                                    <  2:7250327/93‑1 (MQ=255)

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

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

        tGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGGCTCCGGTACTGGTGCGCGTATTGTCACCAGCGGCGATGATGTc                                                                                 >  1:4666301/1‑93 (MQ=255)

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

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

        tGTAGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTc                                                                                 <  1:1293397/93‑1 (MQ=255)

           aGCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGt                                                                              <  1:7474354/93‑1 (MQ=255)

            gCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGtt                                                                             >  2:346760/1‑93 (MQ=255)

            gCTACACGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGtt                                                                             <  1:4409075/93‑1 (MQ=255)

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

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

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

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

               acacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCCGCCGCGATGATGGCTGTttgt                                                                          >  2:5659762/1‑92 (MQ=255)

               acacGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttga                                                                          >  1:2592950/1‑93 (MQ=255)

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

                 acGCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgatt                                                                        <  2:5113793/93‑1 (MQ=255)

                   gCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattga                                                                      >  2:3342156/1‑93 (MQ=255)

                   gCCGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTttgattga                                                                      >  2:1406308/1‑93 (MQ=255)

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

                     cGGATCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGt                                                                    >  1:699626/1‑93 (MQ=255)

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

                       gATCAGGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGttt                                                                  <  2:4241397/93‑1 (MQ=255)

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

                        aTCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTc                                                                 >  1:2126473/1‑93 (MQ=255)

                         tCAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCa                                                                >  2:4210044/1‑93 (MQ=255)

                          cAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAg                                                               <  2:7542836/93‑1 (MQ=255)

                          cAAGGTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAg                                                               <  2:4105999/93‑1 (MQ=255)

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

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

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

                            aGGTGCCACGCGCATTGAACGCGTCCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGt                                                             <  1:92127/93‑1 (MQ=255)

                             ggTTCCACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTg                                                            >  2:5517445/1‑93 (MQ=255)

                                 ccACGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCg                                                        >  1:5595937/1‑93 (MQ=255)

                                  ttcgcgCATTTAACGCGTGTTGGCCTCCGTTATTGGTGCGCGTATTGTCACTAGCCGCGATGATGTCTTTTTGATTGAGTTTCAGGTGCCCGc                                                       >  1:7036991/3‑93 (MQ=255)

                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      >  2:257384/1‑93 (MQ=255)

                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      >  1:6005047/1‑93 (MQ=255)

                                   aCGCGCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGcc                                                      >  1:3998165/1‑93 (MQ=255)

                                    cgggcATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     <  1:302074/90‑1 (MQ=255)

                                    cgcgCATTGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCa                                                     <  2:3205024/93‑1 (MQ=255)

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

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

                                          tttaaCGCGTGCTGGCCTCCGGTACTGGTGGGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAgg                                               <  1:5394492/90‑1 (MQ=255)

                                          ttGAACGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAgg                                               >  1:764157/1‑93 (MQ=255)

                                             accgcgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAtt                                            <  2:7106427/91‑1 (MQ=255)

                                             aaCGCGTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGAtt                                            >  1:5152544/1‑93 (MQ=255)

                                               cgcgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCGGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTc                                          <  2:5236437/93‑1 (MQ=255)

                                               cgcgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTc                                          <  2:1563053/93‑1 (MQ=255)

                                                gcgTGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCa                                         >  2:871865/1‑93 (MQ=255)

                                                   tGCTGGCCTCCGGTTCTGGTGCGCGTATTTTCACCAGCCGCGATGATGTGTGTTTGATTTAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAAc                                      <  1:4078042/93‑1 (MQ=255)

                                                   tGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAAc                                      <  2:5718742/93‑1 (MQ=255)

                                                   tGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAAc                                      >  1:5925727/1‑93 (MQ=255)

                                                   tGCTGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAAc                                      >  2:3828667/1‑93 (MQ=255)

                                                   ggCTGGCCTCCGAAACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGTTTGAGTTTCAGGGGCCTGCCAGTCAGGATTTCAAAc                                      >  2:7898294/2‑93 (MQ=255)

                                                   ggCTGGCCCCCGGTCCTGGTGCGGGTATTGTCCCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAAc                                      <  2:1782223/92‑1 (MQ=255)

                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   <  2:2441101/93‑1 (MQ=255)

                                                      tGGCCTCCGGTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTgg                                   <  1:452784/93‑1 (MQ=255)

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

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

                                                            cctgTACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTTCCCGCCAGTCAGGATTTCAAACTGGCGCATa                             <  2:4797566/90‑1 (MQ=255)

                                                               gTACTGGTGCGCGCATTGTCACCAGCCGCGATGATGTCCGTTTGATTGAGTTTCAGGTGCCCGCCAATCAAACTTTCTAACTCGCGCATAAag                          >  2:6938531/1‑93 (MQ=255)

                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         >  2:6686132/1‑93 (MQ=255)

                                                                tACTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAaga                         >  2:7397654/1‑93 (MQ=255)

                                                                 aCTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAagag                        <  1:6202212/93‑1 (MQ=255)

                                                                  cgggTGCGCGTATTGTCCCCAGCTGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAgaga                       <  2:4317471/91‑1 (MQ=255)

                                                                  cTGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAgaga                       <  2:433576/93‑1 (MQ=255)

                                                                   tGGTGCGCGTATTGTCCCCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      <  2:2373071/93‑1 (MQ=255)

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

                                                                   tGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      >  2:5809359/1‑93 (MQ=255)

                                                                   tGGTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGAt                      >  2:4577992/1‑93 (MQ=255)

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

                                                                    ggTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATc                     >  2:3419544/1‑93 (MQ=255)

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

                                                                     gTGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCg                    <  1:3465932/93‑1 (MQ=255)

                                                                      tGCGCGTATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGa                   >  1:7160987/1‑93 (MQ=255)

                                                                      tGCGCGTATTGTCACCAGCCGCGATGATGTCGGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGa                   <  2:1877616/93‑1 (MQ=255)

                                                                         gcgTATTGTCACCAGCCACGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCa                >  1:7220483/1‑93 (MQ=255)

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

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

                                                                            tATTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAAt             <  1:1716598/93‑1 (MQ=255)

                                                                             aTTGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATAGAACAAATc            >  2:6915215/1‑93 (MQ=255)

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

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

                                                                               tGTCACCAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCt          >  1:4811472/1‑93 (MQ=255)

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

                                                                                   cccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGaaa      <  1:1044725/92‑1 (MQ=255)

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

                                                                                    ccAGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAc     <  1:5300976/93‑1 (MQ=255)

                                                                                      aGCCGCGATGATGTCTGTTTGATTGATTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   <  1:68979/93‑1 (MQ=255)

                                                                                      aGCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgc   >  2:7876065/1‑93 (MQ=255)

                                                                                       gCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCATGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg  >  2:4670750/1‑93 (MQ=255)

                                                                                       gCCGCGATGATGTCTGTTTGATTGAGTTTCAGGTGCCCGCCAGTCAGGATTTCAAACTGGCGCATAAAGAGATCGACCAAATCCTGAAAcgcg  >  2:3158795/1‑93 (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: `‑`